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Abstract 



In this thesis, we describe some recent results obtained in the analysis of two-dimensional 
quantum field theories by means of semiclassical techniques. These achievements represent a 
natural development of the non-perturbative studies performed in the past years for conformally 
invariant and integrable theories, which have lead to analytical predictions for several measur- 
able quantities in the universality classes of statistical systems. Here we propose a semiclassical 
method to control analytically the spectrum and the finite-size effects in both integrable and 
non-integrable theories. The techniques used are appropriate generalization of the ones intro- 
duced in seminal works during the Seventies by Dashen, Hasllacher and Neveu and by Goldstone 
and Jackiw. Their approaches, which do not require integrability and therefore can be applied 
to a large class of system, are best suited to deal with those quantum field theories characterized 
by a non-linear interaction potential with different degenerate minima. In fact, these systems 
display kink excitations which generally have a large mass in the small coupling regime. Under 
these circumstances, although the results obtained are based on a small coupling assumption, 
they are nevertheless non-perturbative, since the kink backgrounds around which the semiclas- 
sical expansion is performed are non-perturbative too. 
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Introduction 



Non-perturbative methods in quantum field theory (QFT) play a central role in theoretical 
physics, with applications in many areas, from string theory to condensed matter. In the last 
years, a considerable progress has been registered in the study of two-dimensional systems. For 
these models, in fact, exact results have been obtained in the particular situations when the 
systems are conformally invariant or integrable. 

Conformal invariance does not merely represent a restriction on the solvable problems. In 
fact, it provides a fruitful link between QFT and statistical mechanics, realized through the 
ideas of scaling and renormalization, which explain the basic features of phase transitions and 
are naturally encoded in the QFT framework. In fact, physical systems at a second order phase 
transition (called critical) are characterized by the divergence of the correlation length, which 
makes unimportant the fine details of their microscopic structure and leads to an organization 
of the critical behaviours in universality classes. As a consequence, the order parameter fluctu- 
ations can be described with the language of conformal field theories (CFT). The corresponding 
dynamics can be exactly solved in two-dimensions, since in this case conformal invariance per- 
mits to determine universal amplitudes and to extract as well the entire spectrum of the theory. 

Furthermore, dynamics in the vicinity of a critical point can be described by perturbed CFT, 
obtained by adding to the action operators which break the conformal symmetry and introduce 
a mass scale in the system, inducing a renormalization group flow. Suitable choices of the 
perturbing operator make the off-critical massive field theory integrable, with consequent elas- 
ticity and factorization of the scattering. In two dimensions, this fact together with a simplified 
kinematics leads to the exact computation of the scattering amplitudes, which encode in their 
analytical structure the complete information about the spectrum. Moreover, the knowledge 
of the S-matrix permits to implement the so-called form factors approach, which makes the 
analysis of off-shell correlators possible. 

Within this program, during the the last years measurable universal quantities for many 
statistical models have been computed. The studied systems include, for instance, the Ising 
model in a magnetic field, the tricritical Ising model, the g-state Potts model, percolation and 
self-avoiding walks. In all other cases, the understanding of two dimensional QFT has been 
reached up to now either by conformal perturbation theory or numerical methods. 

A complete understanding of QFT also requires the control of finite-size effects, both for 
practical and theoretical reasons. In fact, the finite-volume energy spectrum contains a lot of 
information about the properties of the theory in infinite volume. On the one hand, this permits 
to control the systematic error induced by the finiteness of the samples in the extrapolation 
procedure of numerical simulations. On the other hand, scattering data can be extracted from 
the numerical analysis through this correspondence. Moreover, defining a theory in finite volume 
one can explicitly follow the renormalization group flow between its ultraviolet (UV) conformal 
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limit and the infrared (IR) massive behaviour, since these two limits correspond, respectively, 
to volumes much smaller or larger than the correlation length. Finally, QFT in finite volume 
are intimately related to another important subject, i.e. QFT at finite temperature. 

At present, finite-size effects have been studied non-perturbatively only in the above men- 
tioned cases of conformal or integrable QFT. At criticality, in fact, powerful analytical techniques 
are available to extract the entire spectrum of the transfer matrix. Furthermore, in integrable 
massive theories, the knowledge of the exact S-matrix permits to study the system at finite tem- 
perature, with the so-called Thermodynamic Bethe Ansatz. This technique provides integral 
equations, mostly solved numerically, for the ground state energy of the theory in finite volume. 
In all other cases, the control of finite-size effects has been reached either by perturbative or 
numerical methods. 

A natural development of the above mentioned studies of two-dimensional QFT consists in 
looking for some techniques to control analytically non-integrable systems and finite-size effects, 
both for theoretical reasons and their application to several condensed-matter or statistical sys- 
tems. My original contribution to this research stream has been presented in the papers [1-4] 
and will be described in this thesis. The basic tool used is an appropriate generalization and ex- 
tension of semiclassical methods, which have proved to be efficient in analysing non-perturbative 
effects in QFT since their introduction in the seminal works |2.'-}| I24j . The semiclassical approach 
does not require integrability, therefore it can be applied on a large class of systems. At the 
same time, it permits to face problems which are not fully understood even in the integrable 
cases, such as the analytic study of QFT in finite volume. In particular, it has led to new 
non-perturbative results on form factors at a finite volume PP, spectra of non-integrable models 
[2] and energy levels of QFT on finite geometries 11] . 

The semiclassical method is best suited to deal with those quantum field theories charac- 
terized by a non-linear interaction potential with different degenerate minima. In fact, these 
systems display kink excitations, associated to static classical backgrounds which interpolate be- 
tween the degenerate vacua and generally have a large mass in the small coupling regime. Under 
these circumstances, although the results obtained are based on a small coupling assumption, 
they are nevertheless non-perturbative, since the kink backgrounds around which the semiclas- 
sical expansion is performed are non-perturbative too. In any case, the restriction on the variety 
of examinable theories imposed by the above requirements is rather mild, since non-linearity is 
the main feature of a wealth of relevant physical problems. 

In the study of non-integrable spectra and finite-size effects, we have used two basic results. 
The first is the well-known semiclassical quantization technique, introduced by Dashen, Hass- 
lacher and Neveu (DHN) in [23]. It consists in solving a Schrodinger-like problem, associated to 
the "stability equation" for the small quantum fluctuations around the classical backgrounds, 
and in building the energy levels of the system in terms of the classical energy and the stability 
frequencies. The second tool is a covariant refinement of a result due to Goldstone and Jackiw 
|24| . They have shown that the classical kink backgrounds can be interpreted at quantum level 
as the Fourier transform of the form factors between kink states. In pQ, we have overcome the 
drawbacks of the original result, which was expressed non-covariantly in terms of the particles 
space momentum, using the so-called rapidity variable, which is Lorentz invariant. This refine- 
ment opened the way to the use of semiclassical form factors in the analysis of particle spectra 
and correlation functions, because it permits to go in the crossed channel and write the form 
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factor between the vacuum and a kink-antikink state. Since the analytical structure of these 
form factors encodes the information about the masses of all the kink-antikink bound states, in 
this way it is possible to explore non-perturbatively the complete spectrum of non-integrable 
theories, a purpose for which there are no other known analytical techniques. In jl], a predic- 
tion about the spectrum of the cj) 4 field theory in the broken Z2 symmetry phase was made, in 
agreement with previous approximate results. In the following work |2j, a detailed analysis of 
the spectrum was performed in the double Sine-Gordon model, which turns out to be a rele- 
vant description of many concrete physical systems and displays appealing features such as false 
vacuum decay and phase transition phenomena. 

The study of finite-size effects can be tackled after having identified suitable classical so- 
lutions to describe the kinks in finite volume. With this respect, the key result of is the 
construction of such backgrounds for the Sine-Gordon and the broken (j) 4 theories on a cylindri- 
cal geometry with antiperiodic boundary conditions. The form factors are then expressed as a 
Fourier series expansion of the classical solutions, and they are used to obtain an estimate of the 
spectral representation of correlation functions in finite volume. This result adds to few others, 
which however strictly rely on specific integrable structures of the considered models, while it 
can be in principle extended to any theory displaying topologically non-trivial backgrounds. The 
complete DHN quantization has been then performed in [3] for static classical solutions on a 
finite geometry. In particular, the example of the Sine-Gordon model with periodic boundary 
conditions has been explicitly treated, reconstructing the scaling function and the excited energy 
levels in terms of the size of the system. Finally, the same theory, defined this time on a strip 
with Dirichlet boundary conditions, has been studied in [I]. The semiclassical achievements pro- 
vide an explicit analytical control on the interpolation between the massive field theory and its 
UV limit, and can be compared with the numerical results obtained in the sine-Gordon model 
with different techniques. However, their application is more general, and permits to estimate 
the scaling functions also in non-integrable theories. 

The thesis is organized as follows. In Chapter^ we review the non-perturbative results 
available for conformal and integrable QFT in two dimensions, with particular emphasis on their 
applications to statistical mechanics. In Chapter[2] we introduce the semiclassical techniques 
established in the past for QFT in infinite volume, and we test the efficiency of the semiclassical 
approximation by comparing the results obtained in the sine-Gordon model with the exact ones 
provided by the integrability of the theory. Chapter^] is devoted to the analysis of the spectra 
of excitations in the non-integrable systems defined by the (f) 4 field theory in the broken Z2 
symmetry phase and by the Double Sine-Gordon model. Finally, in Chapter0]we present the 
study of finite-size effects. In particular, we describe the form factors in the sine-Gordon and 
broken (j) 4 theories defined on a twisted cylinder, and the energy levels in the sine-Gordon model 
defined both on a periodic cylinder and on a strip with Dirichlet boundary conditions. 
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Chapter 1 

Quantum field theories in 2D and 
statistical physics 

One of the most fruitful applications that QFT has found in recent years is the analysis of 
the universality classes of two-dimensional statistical models near their critical points, which 
correspond to second order phase transitions. Critical systems, in fact, fall into universality 
classes which can be classified by conformally invariant field theories (CFT). The corresponding 
dynamics can be exactly solved in two-dimensions through a systematic computation of cor- 
relators, allowed by conformal invariance. Furthermore, dynamics in the vicinity of a critical 
point can be described by perturbed CFT, obtained by adding to the action operators which 
break the conformal symmetry and introduce a mass scale in the system. Suitable choices of 
the perturbing operator make the off-critical massive field theory integrable, with consequent 
elasticity and factorization of the scattering. In two dimensions, this fact together with the sim- 
plified kinematics leads to the possibility of computing exactly the scattering amplitudes, which 
encode in their analytical structure the complete information about the spectrum. Moreover, 
the knowledge of the 5-matrix permits to implement the so-called form factors approach, which 
makes the analysis of off-shell correlators possible. Within this program, during the the last 
years measurable universal quantities for many statistical models have computed. The studied 
systems include, for instance, the Ising model in a magnetic field, the tricritical Ising model, the 
g-state Potts model, percolation and self-avoiding walks. 

In this Chapter, we will introduce the basic concepts related to conformal invariance and 
integrability. Section lTTl presents a brief overview of conformal field theories, with few examples 
which will be referred to in the following Chapters. In Section fL2l we discuss the main features 
of massive scattering in integrable quantum field theories, exploiting the powerful consequences 
of the simplified two-dimensional kinematics. Finally, Section fL3l describes a useful effective 
Lagrangian description of CFT. 

1.1 Critical systems and conformal field theories 

A statistical mechanical system is said to be critical when its correlation length £, defined as the 
typical distance over which the order parameters are statistically correlated, increases infinitely 
(£ — > oo). Correspondingly, length scales lose their relevance, and scale invariance emerges. 
This is peculiar of the continuous ( or second order) phase transitions, which consist in a sudden 
change of the macroscopic properties of the system as some parameters ( e.g. the temperature) 
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are varied, without finite jumps in the energy (characteristic of first order transitions). A 
remarkable property of these systems is that the fine details of their microscopic structure 
become unimportant, and the various possible critical behaviours are organized in universality 
classes, which depend only on the space dimensionality and on the underlying symmetry. This 
allows a description of the order parameter fluctuations in the language of a field theory, which 
is invariant under the global scale transformations 

provided that the fields transform as 

where A is called the scaling dimension of the field <&. 

The use of conformal invariance to describe statistical mechanical systems at criticality 
is motivated by a theorem, due to Polyakov, which states that local field theories which are 
scaling invariant are also conformally invariant [Sj. Therefore, every universality class of critical 
behaviour can be identified with a conformal field theory (CFT), i.e. a quantum field theory 
that is invariant under conformal symmetry. This way of studying critical systems started with a 
pioneering paper by Belavin, Polyakov and Zamolodchikov [BJ, and is systematically presented in 
many review articles and text books ( see for instance |HJ |§] ) • We will now concisely summarize 
the main properties of CFT. 

An infinitesimal coordinate transformation — > x^ + is called conformal if it leaves 

the metric tensor g^ u invariant up to a local scale factor, i.e. 

9fiu(x) -> g(x)g^ u (x). 

These transformations, which include rotations, translations and dilatations, preserve the angle 
between two vectors and satisfy the condition 

d^ v + d v ^ = ^(d-0, (1.1-1) 

where d is the dimension of space-time. In two dimensions, since the conformal group enlarges to 
an infinite set of transformations, it is possible to solve exactly the dynamics of a critical system, 
assuming conformal invariance and a short-distance operator product expansion (OPE) for the 
fluctuating fields. In fact, if we describe euclidean two-dimensional space-time with complex 
coordinates 

Z — OC I 1/X - Z — X 1/3C - 

ea. (|l.l.l|) specializes to the Cauchy-Riemann equations for holomorphic functions. Therefore 
the solutions are holomorphic or anti-holomorphic transformations, z — ► f(z) and z — > f(z), 
such that dzf = d z f = 0. These functions admit the Laurent expansion 

oo oo 
n=— oo n=— oo 

which has an infinite number of parameters. In this way, the conformal group enlarges to an 
infinite set of transformations. 



Defining now a two-dimensional quantum field theory invariant under conformal transfor- 
mation, we can associate to each field an holomorphic and an antiholomorphic comformal di- 
mensions h and h, defined in terms of the scaling dimension A and of the spin s as 

h = ±(A + s), h = ±(A-s). (1.1.2) 

A field is called primary if it transforms under a local conformal transformation z — > w = f(z) 
as 

(10"* (§) "*«"•*>■ (lx3) 

Conformal invariance fixes the form of the correlators of two and three primary fields up to a 
multiplicative constant: 

2jr=2h if h\ = hi = h and h\ = hi = h 
(Mzi,zi)Mz2,z 2 )) = { Z ^ Z ^ , (1.1.4) 

otherwise 

hizuz^M^^M^zs)) = C l2 3 hl+h2 _ ha h2+ l 3 . hl h 3+hl - h2 (1.1.5) 

z 12 z 23 z 13 

1 



hi+h2 — h 3 -h 2 +hs—hi -hz+hi—h^ ' 



Z 12 

where zij = Z{ — Zj and Zij = Zi — Zj ■ 

It is typical of correlation functions to have singularities when the positions of two or more 
fields coincide. The operator product expansion (OPE) is the representation of a product of 
operators (at positions x and y) by a sum of terms involving single operators multiplied by 
functions of x and y, possibly diverging as x — > y. This expansion has a weak sense, being valid 
within correlation functions, and leads to the construction of an algebra of scaling fields defined 
by 

<lH(x)Mv) = ^2C^ j (x,y)cj )k (y), (1-1.6) 
k 

where C^(x,y) are the structure constants. Translation and scaling invariance forces these 
functions to have the following form: 



C*j(x,y)- 



X - yjAi+Aj-Afc ' 

where CL are exactly the undetermined multiplicative constants of the tree-point correlators 

1EE5J. 

A particularly important operator is the stress-energy tensor T^ 11 , which expresses the vari- 
ation of the action under a transformation of coordinates x^ — ► x^ + £ M (x): 



Conformal invariance is equivalent to the vanishing of 5S under the condition ()1.1.1|) . and it 
is guaranteed by the tracelessness of the stress-energy tensor. Together with translation and 
rotation invariance (d^T^ = 0), the condition Tjf = is expressed in complex coordinates as 



d z T = and d z T = 0, 



where T(z) = T\\ — T22 + 2zTi2 and T(z) = T\\ — T22 — 2iT±2- Therefore the stress-energy 
tensor splits into a holomorphic and an antiholomorphic part. In two dimensions, it is possible 
to deduce the following OPE for the stress-energy tensor and a primary field of dimension (h, h): 

T(z) 4>(w, w) = 7 ttt 4>(w, w) H d w 6(w, w) + regular terms , (1.1.7) 

(z — w) z z — w 

T(z) d>(w, w) = — —k 4>(w, w) + dw4>(w, w) + regular terms . 

(z — w) z z — w 

Furthermore, it is possible to show that the OPE of the stress-energy tensor with itself has the 
form: 

T(z)T(w) = . C//2 . . + 7 — ~ ttt T(w) + dT(w) + regular terms , (1.1.8) 

(z — wp (z — w) 1 z — w 

where the constant c, called central charge, depends on the specific model. A similar expression 
holds for the antiholomorphic component. The holomorphic and antiholomorphic components 
of the stress-energy tensor can be expanded in Laurent series respectively on modes L n and L n , 
which are the quantum generators of the local conformal transformations 

n*) = E . m = £ J£j , (1.1.9) 

n=— 00 n=— 00 

and obey the Virasoro algebra 

[L n , L m ] = (n - m) L n+m + — n(n 2 - 1) 5 n+rnfi , 

[L n , L m ] = (n - m) L n+m + ^ n(n 2 - 1) 5 n+mfl , (1.1.10) 

[L n , L m ] = . 

In virtue of the decomposition of (jl.l.lflj) in the direct sum of two algebras, one in the 
holomorphic and the other in the antiholomorphic sector, the general properties of CFT have 
the same form in the two sectors, and from now on we will only restrict to the holomorphic part. 

Comparing definition (|1.1.9j) with the OPE Ql.1.7]) . we can deduce the action of some gener- 
ators on a primary field: 

(Lo^iz) = h<j>{z) 
{L^){z) = d<j>{z) (1.1.11) 
{L n <f>) (z) = if n > 1 

The relation [Lo,L n ] = —nL n leads to the interpretation of generators L n with n > as 
destruction operators and with n < as creation operators. Hence primary fields define highest 
weight representations of the Virasoro algebra, being annihilated by all destruction operators. 
The action of creation operators on these fields is encoded in the regular part of the OPE (|1. 1.7(1 . 
and defines the so-called descendant fields 

^i^--^) = (L_ ni L_ n2 ...L_ nfe )0, 

which are again eigenvectors of Lq: 



The number N = Y^i=i rii is called level of the descendant. As an example, the stress-energy 
tensor is a level two descendant of the identity (T = L_2l). The set [(f)] constituted by all 
the descendant fields of a primary operator (f> is called conformal family. It is possible to show 
that every correlation function involving descendant fields can be computed by applying a linear 
differential operator to the correlation function of the corresponding primary fields. 

The Hilbert space of states of a CFT is built by acting on the vacuum with the operators 
evaluated at z = 0. Therefore, the primary states are given by 

\h) = 0(O)| 0) , 

and the descendent states can be obtained from them as L_ rtl L_.„ 2 ...L_ n J h). 

In concluding this section, it is worth mentioning how the central charge c has the physical 
meaning of measuring the response of the system to the introduction of a macroscopic scale 
|10| I12j . In fact, the complex plane can be conformally mapped to an infinite cylinder of 
circumference R by the transformation (see Fig . 11 . 1 j) 

R 

z _> w( z ) = — Inz . (1.1.12) 
Implementing the above transformation on the stress-energy tensor one gets 

T cyL (w) = (j?J [t pL (z)z 2 - ^ . 

If we assume that the vacuum energy density (T p i) vanishes on the plane, we see that it is 
non-zero on the cylinder: 

(Tcyl.) = -g^2" 

The central charge is then proportional to the Casimir energy, which naturally goes to zero as 
the macroscopic scale R goes to infinity. In particular, the hamiltonian and the momentum are 
expressed on the cylinder in terms of the Virasoro generators as 

H = ^{ L o + Lo-^) P= 2 -^-(L -L ) (1.1.13) 



z — > w(z) 



X 



R 



Figure 1.1: Conformal map (|1.1.12j) from plane to cylinder. 



1.1.1 Examples 

We will now describe few examples of CFT, which will be further discussed in the thesis. We 
limit ourselves to the statement of the main results, whose proofs can be found in the literature. 



Gaussian CFT 

An interesting CFT, only apparently trivial, is given by the gaussian action which describes a 
free bosonic field 



A G = I d 2 x 8^8^, (1.1.14) 



where g is a normalization parameter. The two-point function of the field (ft is expressed as 

( (ft(z)(ft(w) ) = log(;z — w) + const. . 

Airg 

Using Wick's theorem, one can easily show that the OPE between the energy-momentum tensor 
T(z) = —2i:g : deft deft : and the operator deft is given by 

m/ \r>,/ \ d(ft(w) d 2 (ft(w) . 

1 [z)o(ft{w) = 7 + 7 ^ + regular terms , 

(z — w) z (z — w) 

from which it follows that d(ft is a primary field with conformal dimensions h = h = 1. Further- 
more, the OPE of T with itself 

T(z)T(w) = — rj + + 7 + regular terms 

(z — w;) 4 (z — u;) 2 (z — w) 

shows that the theory is characterized by a central charge c = 1. Besides d(ft, other primary 
operators can be built by normal ordering the exponentials of (ft: 



V a (z,z) 



■ p iouj>{z,z) . 



_ 2 

These fields, called vertex operators, have conformal dimensions h = h = 

A generalization of the above theory is obtained compactifying the bosonic field, i.e. making 
it an angular variable on a circle of radius 1Z. This can be done by defining the theory on a 
cylinder of circumference R with boundary conditions 

<f>(x + R,t) = (ft(x,t) + 2-TrnK . (1.1.15) 

The index n, called winding number, labels the various sectors of this CFT, together with an 
index s related to the discrete eigenvalues of the momentum operator on the circle. In each 
sector, the states | s, n) with lowest anomalous dimension are created by the vertex operators 

V 8 , n (z,z) =: exp [ia+ n tp(z) + ia~ n (p{z)] :, (1.1.16) 

i.e. 

| s,n) = V s>n (Q,0) | vac) , 

where 

a tn = ± 27T 9 nil ; 

(ft{x,t) = ip(z)+(p{z) . 



Their conformal dimensions are given by 
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h sn = 2itg[- — + -nTZ) , h sn = 2itg[- — --nil) . (1.1.17) 

\4ngK 2 J ' ' y \^gTZ 2 



Minimal models 



Virasoro minimal models are particular CFT characterized by a finite set of conformal families, 
in virtue of a truncation of the operator algebra. These theories can be labelled as A4(p,p') 
with two integers p and p' , in terms of which the central charge and the conformal dimensions 
of primary fields are expressed as 



1-6 



(p-p) 



'\2 



VP 

(pr - p' s) 2 - (p -p') 2 
App 1 



(1.1.18) 
(1.1.19) 



with 



1 < r < p' 



and 



1 < s < p 



The conformal dimensions are organized in a rectangle in the (r, s) plane, called Kac table. The 
number of distinct fields is (p — l)(p' — l)/2, since there is a symmetry h TfS = hp/_ r)P _ s which 
makes half of the Kac rectangle redundant. It can be shown that minimal models are unitary 
only if \p — p'\ = 1, and in this case they are usually labelled with p 1 = m and p = m + 1. 

The simplest unitary minimal model is -M(4, 3), which has central charge c= \ and the Kac 
table shown in Fig. ll.2l 



I 



a 



a 



12 3 
Figure 1.2: Kac table of the minimal model M(&, 3) 

This field theory is in the same universality class as the lattice Ising model defined by 
the usual configuration energy 

E[a] = - J^^UiUj - h^ai , ate {-1,1}. 

Besides the identity operator 0x,i = I) the theory contains the operator <j)\^ = a with conformal 
dimension h a = jq, which is the continuum version of the lattice spin crj, and 4>\ : s = e, with 
h e = 2 j which corresponds to the interaction energy <7j<7j + i. The algebra defined by the OPE 
(|1.1.6jl can be schematically represented in this model by the fusion rules 



a x a ~ I + 
a x e ~ a , 

£ X £ ~ I . 



This notation means, for instance, that the OPE of a with a (or fields belonging to their families) 
may contain terms belonging only to the conformal families of I and e. 
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Figure 1.3: Kac table of the minimal model A^(5,4) 

The next unitary model, .M(5,4), displays a richer structure. The Kac table of this model, 
which has central charge c = jq, is shown in Fig. 11.31 

It was recognized in |13j that the lattice model associated with this conformal field theory 
is the dilute Ising model at its tricritical fixed point (TIM), defined by 

E[a,t] = -J^aidjUtj - n^2(ti - 1) , Oi G {-1,1}, U G {0,1} , 

where [i is the chemical potential and U is the vacancy variable. The corresponding phase 
diagrams is drawn in Figure 11.41 where I and II denote respectively a first and second order 
phase transition, and the point (J/, 0) represents the Ising model, with all lattice's site occupied. 
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Figure 1.4: Phase diagram of the TIM. 



The field — with h £ = j^, corresponds to the energy density, while 4>i^ = e', with h e > = 
is the vacancy (or subleading energy) operator. The leading and subleading magnetization 
fields are respectively 02,2 = c arid 4>2,i = cr', with h a = and h a > = j^. The remaining field 
0x,4 = e" has conformal dimension h E n = |. Dividing the operators in even and odd with respect 
to the Z2 symmetry of the model under Oi — > —0%, we can list the fusion rules in the following 
way: 
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1.2 Integrable quantum field theories 

The scaling region in the vicinity of second order phase transitions can be described by a given 
CFT perturbed by its relevant operators (characterized by an anomalous dimension A, < 2), 
with the action 

A = Acft + ^\ d 2 x^(x), (1.2.1) 

i ^ 

where the couplings have mass dimension Aj ~ [m] 2_Ai . The relevant operators, being super- 
renormalizable with respect to UV divergencies, do not affect the behaviour of the system at 
short distances, but they change it at large scales. Any combination of the relevant fields defines 
a Renormalization Group (RG) trajectory which starts from the given CFT and can reach an- 
other critical point (defined by a different CFT) or a non-critical fixed point, corresponding to 
a massive QFT. From now on, we will only consider this second case. It was shown in [Tlj that, 
depending on the choice of the perturbing operator, the off-critical massive field theory can be 
integrable, with consequent elasticity and factorization of the scattering. Obviously, this kind 
of theories covers only a class of the statistical systems of interest, however this class includes 
relevant physical problems. For instance, integrable QFT correspond to the Ising model with 
thermal or magnetic perturbation separately, described respectively by 




or to the tricritical Ising model perturbed by the leading energy operator 




We will now present a brief overview of integrable massive quantum field theories, underlying 
their basic features (for an exhaustive review, see JH])- This will also give us the opportunity 
of introducing the most important kinematical quantities used in the following. 

An integrable QFT is characterized by the presence of an infinite set of conserved charges, 
which make the corresponding scattering theory purely elastic and factorized |16| . This implies 
that an arbitrary n-particle collision process can be described by the product of n(n — l)/2 
elastic pair collisions. Hence the determination of the complete S matrix reduces to that of the 
two-particle amplitudes, which are defined as 

\Ai (pi) A, (p 2 )) in = Sfj \A k (pg) A l ( Pi )) out , (1.2.2) 

where Ai (p\) and Aj (p 2 ) denote the incoming particles (with 2-momenta p± andpg), and A\. (p^) 
and Ai (p^) the outgoing states (see Fig. ll.5jl . 



Figure 1.5: Two-particle S'-matrix 



Lorentz invariance fixes the two body S-matrix to be a function of the Mandelstam variables 
s = (Pi + P2) , t = (pi — ps) and u = (pi — pi) , which satisfy the relation s + t + u = J2i=i m i- 
Since in (1+1) dimensions and for elastic scattering only one of these variables is independent, it 
is convenient to introduce a parameterization of the momenta in terms of the so-called rapidity 
variable 0: 

p\ = rrii cosh 0j , p\ = mi sinh B% , (1.2.3) 

which corresponds to the following expression for the Mandelstam variable s: 

s = (pi + p 2 )^ (pi + P2Y = m 2 +m 2 + 2mimj cosh 0y , (1-2.4) 

with Oij = Oi — Oj. The functions S!?J will then depend only on the rapidity difference of the 
involved particles: 

\A i (e 1 )A j (e 2 )) in = s%(e 12 ) \A k {e 2 )A l {e 1 )) (mt . 

Elasticity of the scattering processes implies a drastic simplification in the analytic structure 
of the S-matrix, which can be extended to be an analytic function in the complex s-plane 
|17j . In fact, contrary to the generic case, where many branch cuts are present, in the elastic 
case the two-particle S-matrix only displays two square root branching points at the two- 
particle thresholds (m; — rrij) 2 and (rrij + mj) 2 , and is real valued on the interval of the real 
axis between them. From (|1.2.4|) it follows that the functions Sfj (9) are meromorphic in 8, and 
real at Re(#) = 0. The two cuts in the s variable, in fact, are unfolded by the transformation 
(jl. 2.4(1 : for instance, the upper side of the cut along [(m, + m j) 2 j 00 ] is mapped into the positive 
semiaxis < < 00, while the lower side is mapped into the negative semiaxis — 00 < 8 < 0. The 
physical sheet of the s-plane goes into the strip < Im(#) < ir, while the second Riemann sheet 
is mapped into — ir < Im(8) < 0. The structure in the plane repeats then with periodicity 2iri. 
See Fig. ll.6l for a representation of the analytic structure of the S'-matrix in the two variables s 
and 0. 

The two-particle S-matrices satisfy the usual requirements of unitarity, expressed as 

E^TWSL = (1-2-5) 

n,m 

and crossing symmetry, given by 

flg(fl) = flg(m-e), (1.2.6) 

since the analytic continuation s — > t from the s-channel to the t-channel corresponds to the 
change of variable — > in — 0. Furthermore, the amplitudes are restricted by the star-triangle 
(or Yang-Baxter) equations 

(012) 5gj (9 l3 ) fig (0 23 ) = (0 l3 ) S{% (0 12 ) fiS ^) > ( L2 - ? ) 
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Figure 1.6: Analytic structure of the elastic S-matrix in the variables s and 9. 



where a sum over the intermediate indices is understood. The system of equations (jl.2.5j) . Q1.2.6j) 
and (|1.2.7j) is in many cases sufficient to solve the kinematics of the problem, determining 
a consistent solution for the two-particle S'-matrix, up to a so-called CDD ambiguity, which 
consists in multiplying the solution by factors that alone satisfy the same equations. 

Since the bound states of a theory correspond to singularities of the S'-matrix, its analytic 
structure encodes the spectrum of the system. Stable bound states are usually associated to 
simple poles in the s variable located in the real interval (m; — rrij) 2 < s < (m; + mj) 2 . These 
are mapped by Q1.2.4|) onto the imaginary #-axis of the physical strip (see Fig. ll.6|) . If a two- 
particle amplitude with initial states Ai and Aj has a simple pole in the s-channel at 9 = iu^j 
(A n is the associated intermediate bound state), in the vicinity of this singularity we have 

iR n - 

S$(0) ~ 7 tJL T > with K = r tf r » (1.2.8) 



where Tfj are the on-shell three-particle coupling constants of the underlying quantum field 
theory (see Fig. ll.7|) . Remembering that the corresponding singularity in the s variable is of the 
form (s — m^) and using relation (|1.2.4[) . we get the following expression for the mass of the 
bound state: 

,2 2 , 2 



m„ = m 4 



+ rrij + 2 mi rrij cos ufj. (1.2.9) 




Figure 1.7: First-order pole in the S-matrix 

The dynamics of the system can be determined by implementing the so-called "bootstrap 
principle", which consists in identifying the bound states with some of the particles appearing 
as asymptotic states. This leads to further equations which permit to fix the CDD ambiguities 
mentioned above and to identify the particle content of the theory. 



Finally, it is worth recalling that unstable bound states (resonances) are associated to s- 
variable poles in the second Riemann sheet at s = (jn^ — , where is the inverse life-time 
of the particle. These correspond to poles in 9 located in the strip —ir < Im(8) < at positions 
9 = — iufj + a& satisfying 

r 2 

m\ — = mf + to 2 + 2 rrti rrij cos ufj cosh cqj , 

mk Tfc = luiirrij sin it^- sinhajj- . (1.2.10) 



The knowledge of the exact S-matrix further permits to extract non-perturbative informa- 
tion on the correlation functions of the theory. In fact, the spectral function p^' associated to 
a given operator defined as 



(0|$(x)$(0)|0) 



^p(*V)e^ 



(27T) 



can be expanded as a sum over complete sets of particle states 



P 



*V) = 27rE^/ ^i-^n S(p -p 1 ...-p n )5(p 1 -p\...-p 1 n ) X 



where d£l 



dp _ dB_ 
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and 



a\...a n \" Li ■ ■ ■ i 



9 n ) = (0\<S>(0)\A ai (6 1 )...A an (e n )} 



(1.2.11) 



(1.2.12) 



The matrix elements (|1.'2.12|) . called "form factors" and pictorially depicted in Fig. 11.81 are 
subject to the Watson equations, which relate them to the scattering amplitudes. In the case 
of integrable theories these equations take a simplified form |18( I19[ I20| . which permit the 
determination of form factors once the S'-matrix is known. Furthermore, it has been shown in 
a series of works (see, for instance |2j) that the spectral representation has a fast convergent 
behaviour, therefore accurate estimates of correlators and other related physical quantities can 
be obtained by just using few exact terms in the series (|1.2.11|) . having consequently a great 
simplification of the problem. 




Figure 1.8: Form factor (|1.2.12|) 



Here we limit ourselves to the description of the main equations ruling the simplest non- 
trivial form factors, i.e. the two-particle ones, that we will use and discuss in the following. The 



discontinuity of the matrix elements across the unitarity cut leads to the following relation with 
the two-particle scattering amplitudes: 

F%{6) = S$(6)F? d {-6) , (1.2.13) 

where 9 = 9\ — 9i- The crossing symmetry of the form factor is expressed as 

F*(0 + 2*vr) = e- 2 ^Ft(-9) , (1.2.14) 

where the phase factor e - 2j7r 7*,a j s inserted to take into account a possible semi-locality of the 
operator which interpolates the particle a (i.e. any operator ip a such that (0|</? a |a) 7^ 0) with 
respect to the operator $(x). When 7$^ = 0, there is no crossing symmetric counterpart to 
the unitarity cut but when 7$^ 7^ 0, there is instead a non-locality discontinuity in the plane 
of the Mandelstam variable s defined in 1)1.2.4(1 . with s = as branch point. In the rapidity 
parameterization there is however no cut because the different Riemann sheets of the s-plane 
are mapped onto different sections of the #-plane; the branch point s = is mapped onto the 
points 9 = ±in which become therefore the locations of simple annihilation poles. The residues 
at these poles are given by 

Res e=±iw F*(9) = iS- ab (l - e^ 2i ^)(0\m . (1.2.15) 

Finally, the two-particle form factors inherit the s-channel bound state poles of the S-matrix, 
and the corresponding residua are given by 

Res e=iu o b F* b (e)=ir c ab F? . (1.2.16) 

where the couplings T c ab coincide with the ones in 1)1.2.8(1 . This relation, pictorially represented 
in Fig. 11.91 is the property of form factors that we will mostly exploit in the following. 




1.3 Landau— Ginzburg theory 

As we have seen in Sect. ll.ll statistical models at the critical point can be described by conformal 
field theories (CFT), and many systems of physical relevance have been identified with the 
Virasoro minimal models. The operators of these theories can be organized in a finite number 
of families, and this simplifies the dynamics allowing in principle for a complete solution. The 
only disadvantage of this kind of theories is that they have no Lagrangian formulation, therefore 
they cannot be studied by a path-integral approach and the underlying physics is not always 



transparent. However, for a class of minimal theories, the unitary ones M{m + l,m), there is a 
simple effective Lagrangian description, suggested in [22], which is realized by a self-interacting 
field (f> subjected to a power-like potential. The field <f> stands for the order parameter of the 
statistical system, and the potential V((f>), whose extrema correspond to the various critical 
phases of the system, is usually chosen to be invariant under the reflection <j> — > —<j). For a 
potential of degree 2(m — 1), this ensures the existence of m — 1 minima separated by m — 2 
maxima. Several critical phases of the system can coexist if the corresponding extrema coincide. 
The most critical potential is therefore a monomial of the form 

V m ((f>) = 2(m ~ 1} • 

Starting from the field <f>, one can construct composite fields : <p k : by normal ordering its powers. 
These have been shown in |22j to display the same fusion properties as the operators present in 
the minimal model A4(m + 1, to), supporting in this way the correspondence. In particular, the 
field (f) is always associated to the primary operator <^>2,2- 

One of the nicest features of the Landau-Ginzburg description is that it provides a very 
intuitive picture of the perturbation of CFT away from the critical points, since this simply 
corresponds to adding opportune powers of <fi to V m . 

For instance, the universality class of the Ising model is described by to = 3, the spin operator 
corresponds to a ~ 4> and the energy operator to e ~: 4> 2 :. Therefore the thermal perturbation 
of the critical point is described by the Landau-Ginzburg theory 

V{4>) = A(f> 4 + B<f, 2 + C , 

where the sign of B refers to high or low temperature, respectively (see Fig. ll.lU]) . 



(a) V{<f>) (b) V{4>) (c) V{<t>) 




Figure 1.10: Landau-Ginzburg potential for the Ising model (a) at the critical point, (b) at high 
temperature, (c) at low temperature. 



Another very interesting example is given by the case m = 4, which corresponds to the 
tricritical Ising model. The operator content of the theory is in this case identified as 

a ~ (f> 
e ~: <j) 2 : 
< a' ~ : (J) 3 : 
e' ~: ^ : 



Hence, for instance, the perturbation of the tricritical point by leading and subleading energy 
densities is described as 

V(<f>) = A</> 6 + B<t> 4 + C(t> 2 +D , 
and some of the resulting potentials are shown in Fig. ll.lll 

(a) V(4>) (b) V(4>) (c) V(</>) 




Figure 1.11: Landau-Ginzburg potential for the tricritical Ising model (a) at the tricritical 
point, (b) with positive leading and subleading energy perturbations, (c) with positive leading 
and negative subleading energy perturbations. 

In concluding this Section it is worth to remark how the above mentioned correspondence 
between Lagrangians and minimal models is based on euristic arguments which are not rigorously 
proven. Furthermore, quantum corrections can be very strong, hence a quantitative analysis of 
statistical system is generally not possible starting from the Lagrangian formulation. However, 
the qualitative picture offered by the Landau-Ginzburg theory, being essentially based on the 
symmetries of the model, remains valid also at quantum level. In particular, this description 
nicely illustrates the underlying physics in cases when the potential displays degenerate minima, 
since the classical solutions interpolating between the minima have a direct meaning of quantum 
particles, as we are going to show in the next Chapter. 



Chapter 2 

Semiclassical Quantization 



In this Chapter we will describe the two basic technical tools used in the following to investi- 
gate non-integrable spectra and finite-size effects. The first is represented by the semiclassical 
quantization technique introduced for relativistic field theories in a series of papers by Dashen, 
Hasslacher and Neveu (DHN) [22] by using an appropriate generalization of the WKB approxi- 
mation in quantum mechanics. The second is a result due to Goldstone and Jackiw [23], which 
relates the form factors of the basic field between kink states to the Fourier transform of the 
classical solution describing the kink. 

Here we only focus on the description of these two results, that we will fully exploit in the 
following. However, it should be mentioned that the semiclassical study of QFT was developed 
during the Seventies by few independent groups and with complementary techniques, such as the 
path-integral approach, the inclusion of fermions and the analysis of higher order corrections in 
the semiclassical expansion. It is worth mentioning, for instance, the work of Callan and Gross 
[25] , Gervais, Jevicki and Sakita [22], Faddeev and Korepin [22]. A complete review of these 
beautiful achievements can be found in |28j . 

We will restrict to the study of two-dimensional theories, in virtue of their simplified kine- 
matics that allows for a powerful applications of the semiclassical techniques. However, it is 
well known that the semiclassical methods are naturally formulated for QFT in any dimension 
(d+ 1), and have provided, for instance, relevant insight in the study of four-dimensional gauge 
theories. In fact, the main feature that makes interesting their use is the presence of non-linear 
interaction, which gives rise to topologically non-trivial classical solutions. Classical non-linear 
equations have been always recognized to play a crucial role in the description of physical phe- 
nomena, in virtue of their intriguing complex features (see for instance |3(Jl 131 j ) . The expectation 
of finding interesting properties also in the associated quantum field theories, first acquired by 
Skyrme [H2] , opened the way to the non-perturbative semiclassical studies that we are going to 
describe. 

The Chapter is organized as follows. In Section l2~T1 we introduce the DHN quantization 
technique for static non-perturbative backgrounds. In Section l2~2l we describe the Goldstone 
and Jackiw's result and its relativistic formulation, proposed in pQ, which is the basic tool used 
in the following to study the spectrum of non-integrable QFT. Finally, in Section l2~3l we apply 
the above techniques to the sine-Gordon model, where, in virtue of the integrability of the 
theory, exact result to be compared with the semiclassical ones are available. 
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2.1 DHN method 



The semiclassical quantization of a field theory defined by a Lagrangian 



£ = 7i(d^)(d^)-V(cf>) (2.1.1) 



is based on the identification of a classical background (j)d( x , t) which satisfies the Euler-Lagrange 
equation of motion 

d IM d^ d + V'(cf >d ) = 0. (2.1.2) 

The procedure is particularly simple and interesting if one considers finite-energy static classical 
solutions (pci( x ) hi 1+1 dimensions, usually called "kinks" or "solitons". Their presence is 
one of the main features of a large class of 1+1-dimensional field theories defined by a non- 
linear interaction V{4>) displaying discrete degenerate minima (pi, which are constant solutions 
of the equation of motion and are called "vacua" . The (anti)kinks interpolate between two next 
neighbouring minima of the potential, and consequently they carry topological charges Q = ±1, 
where 

oo 

Q = - — l —— I dx^±. (2.1.3) 

ax 



•H+l 



The conservation of this quantity, which can be easily deduced from the associated conserved 
current 

oo 

f = e^d v <j>, Q = 1 t [dxf, 

<Pi+l ~ Pi J 
— oo 

will play a crucial role in the following. 

Being static solutions of the equation of motion, i.e. time independent in their rest frame, 
the kinks can be simply obtained by integrating the first order differential equation related to 

(En 

U9 t) 2 = V(<f>d)+A, (2.1.4) 

further imposing that 4> c i(x) reaches two different minima of the potential at x — > ±oo. These 
boundary conditions, which describe the infinite volume case, require the vanishing of the in- 
tegration constant A. As we will see in the following, the kink solutions in a finite volume 
correspond instead to a non-zero value of A, related to the size of the system. 

For definiteness in the illustration of the method, we will focus on the example of the c/> 4 
theory in the broken Z2 symmetry phase, defined by the potential 

^) = ^ 4 "^ 2 + ^- (2-1.5) 

This theory displays two degenerate minima at <f>± = an d a static (anti)kink interpolating 

between them 

vn vn 

(f>cl(x) = (±) -j= tanh— (x-x ) , (2.1.6) 

where the arbitrariness of the center of mass position xq is due to the translational invariance 
of the theory. The corresponding classical energy is obtained by integrating the energy density 



oo 

£d = J dx e d (x) 



2 v / 2m 3 

~3 r 



(2.1.7) 



Fig. 12. II shows the potential, the classical kink and its energy density. 
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Figure 2.1: Potential V(tp), kink 4> c i(x) and energy density e c i{x). 



The kinks exhibit certain particle properties. For instance, they are localized and topo- 
logically stable objects, since the conservation of topological charge forbids their decay into 
mesons with Q = 0. Moreover, in integrable theories their scattering is dispersionless and, in 
the collision processes, they preserve their form simply passing through each other. All the 
above properties are an indication that the kinks can survive the quantization, giving rise to 
the quantum states in one-particle sector of the corresponding QFT. As it is well known, one 
cannot apply directly to them the standard perturbative methods of quantization around the 
free field theory since the kinks are entirely non-perturbative solutions of the interacting theory. 
Their classical mass, for instance, is usually inversely proportional to the coupling constant (see 
((2, 1.7(1 in the broken (j) 4 example). In infinite volume, an effective method for the semiclassical 
quantization of such solutions (as well as of the vacua ones) has been developed in a series of 
papers by Dashen, Hasslacher and Neveu (DHN) [23] by using an appropriate generalization of 
the WKB approximation in quantum mechanics. 

The DHN method consists in initially splitting the field 4>(x, t) in terms of the static classical 
solution (which can be either one of the vacua or the kink configuration) and its quantum 
fluctuations, i.e. 



<j>(x,t) = <t>d(x) +7](x,t) , rj(x,t) = ^\ 



Vk(x) 



and in further expanding the action of the theory in powers of t], obtaining for instance, in the 
example (|2. 1.5(1 . the expansion in A 



St 



dx dt C(4>ci)+ 



dx dt — r)(x, t) 



dt 2 



dx 2 



m z + 3A<^ i](x,t)+\ 



i J dx dt 



>d rf + 



(2.1.8) 

The semiclassical approximation is realized by keeping only the quadratic terms, and as a result 
of this procedure, %(x) satisfies the so called "stability equation" 



d 2 



<h:l 



ri k (x) = ulr) k (x) , 



(2.1.9) 



together with certain boundary conditions. The semiclassical energy levels in each sector are 
then built in terms of the energy of the corresponding classical solution and the eigenvalues Ui 
of the Scrodinger-like equation 1)2.1. 9JI . i.e. 



E {n . } = E d + h (n k + ^ oj k + 0(h 2 ) , (2.1.10) 

where n/% are non-negative integers. In particular the ground state energy in each sector is 
obtained by choosing all = and it is therefore given by 1 

E = + | J> fc + 0(fi 2 ) . (2.1.11) 
k 

The semiclassical quantization technique was applied in to the kink background 1)2.1.6)1 . 
in order to compute the first quantum corrections to its mass, given at leading order by the 
classical energy. The stability equation 1)2.1.9)1 can be cast in the hyper geometric form in the 
variable z = |(1 + tanh tt|), and the solution is 



,2 



v ( x ) = Z V 1_ 5^(1 - zyV 1 -^? F f 3, -2, 1 + 2y 1 - ^ ; zj . 

The corresponding spectrum is given by the two discrete eigenvalues 

w 2 = 0, with 7ft) (a) = 2 mx , (2.1.12) 

cosrr ^ 

and 

k>l = 7T"^, witn = TT^ ' (2.1.13) 



3 sinh ^ 

- m 2 , with = n 

2 /1V 7 cosh 2 ^ 

plus the continuous part, labelled by q € E, 



2 2 

uj = m 



2+ig 2 ) , with rto(jc) =e^ (^3tanh 2 ^|- 1 - g 2 - 3i 9 tanh 



(2.1.14) 

The presence of the zero mode ujq is due to the arbitrary position of the center of mass xo in 
1)2.1.6)1 . while lo\ and uj q represent, respectively, an internal excitation of the kink particle and the 
scattering of the kink with mesons 2 of mass y/2m and momentum mq/y/2. This interpretation 
will find a clear explanation in Sect. 12.21 

The semiclassical correction to the kink mass can be now computed as the difference between 
the ground state energy in the kink sector and the one of the vacuum sector, plus a mass 
counterterm due to normal ordering (see Fig. 12. 2)1 : 
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Figure 2.2: Tadpole diagram giving rise to 5m 2 . 

The discrete values q n and k n are obtained putting the system in a big finite volume of size 
R with periodic boundary conditions: 

„ mR 
2mr = k n R = Qn—^ + 0{q n ) , 

where the phase shift 5(q) is extracted from rf q (x) in 1)2. 1.141) as 



7} q (x) el^vl^M , %) = -2arctan' 39 



x— >±oo wy V2 — (/ 2 

Sending i? — ► oo and computing the integrals we finally have 

2y/2 m 3 ( 1 /3 3 \ . 

In concluding this section, it is worth mentioning that the construction of the complete 
Hilbert space requires to consider time-dependent multi-kink solutions with finite energy. Their 
semiclassical quantization can be performed with an appropriate modification of the DHN 
method [2B] to include the time-dependence. Moreover, the semiclassical computations can 
be extended at higher order in the small-coupling expansion, keeping cubic (and higher) powers 
of r] in (|2.1.8[) . This procedure is highly non-trivial, especially due to the presence of the zero 
mode wo = 0, which causes some divergencies at higher order in A. However, effective methods 
for solving problem have been developed and applied in |25U26| . in the context of a path integral 
technique. 

2.2 Classical solutions and form factors 
2.2.1 Goldstone and Jackiw's result 

A direct relation between the kink states and the corresponding classical solutions has been 
established by Goldstone and Jackiw [23; w h° have shown that the matrix element of the 
field (p between kink states is given, at leading order in the semiclassical limit, by the Fourier 
transform of the kink background. A detailed explanation of this result can be found in |28j . 
and a similar achievement in a non-relativistic context is presented in |29| . 

From now on we will fix h = 1, since the semiclassical expansion in h is equivalent to the expansion in the 
interaction coupling A. 

2 The mesons represent the excitations over the vacua, i.e. the constant backgrounds 4>± — ±77^, therefore 
their square mass is given by V"(<f>±) = 2m 2 . 



The technique to derive this result relies on the following assumptions, which consolidate the 
basic ideas already exposed in Sect. l2.ll and that we illustrate here in the case of the broken 4 
theory (f2~T3|) : 



1. The Hilbert space of the theory contains, in addition to the 'vacuum sector' (i.e. the 
vacuum and multi-meson states), the so-called kink sector. This is spanned by the states 
| p ) and | p*) , representing the kink particle of momentum p in its ground state and 
excited state, respectively, and the states | p, k±, k m } and | p*,kx, ...,k m ), representing 
the scattering states of the kink particle and m mesons. 

2. The kink sector is orthogonal to the vacuum sector, hence the kink is stable against decay 
into mesons. In fact, although the states in the kink sector have higher energy than 
the vacuum sector ones, they are prevented from decaying by the conservation of the 
topological charge, which holds also in the quantum theory. 

3. The mass of the quantum kink behaves as M ~ ^ in the weak coupling limit (see (J2.1.7J) ) . 

4. The matrix elements of 4>(x, t) between the above states behave in the weak coupling limit 

as 

(p'\Hp) ~ o(i/Va) 

{p | (j) | p, k } and (p' \ <j) \ p*) ~ O (1) 

( P ',k' 1 ,...,k' l \0\p,k l ,...,k m ) ~ o (a^™- 1 )/ 2 ) 

This assumption, which will find confirmation a posteriori, relies on the fact that the kink 
classical background itself is of order 1 / \/A, and that the emission or absorption of every 
additional meson and the internal excitation of the kink carry a factor This can be 
intuitively understood by noticing that, in the expansion Q2.1.8JI of the interaction V(4>), 
the leading perturbative term is of order \<f) c i, i.e. of order y/X. 

Let's now define the matrix element of the basic field 4>{x,t) between in and out one-kink 
states, with momenta p\ and p2, as the Fourier transform of a function /(a), to be determined: 

(P2 | HO) | pi ) = j da e ^-^ a f(a) . (2.2.1) 
Next we consider the Heisenberg equation of motion for the quantum field 4>(x, t) 

{3 2 t - dl) <l>(x, t) = m 2 <j>(x, t) - X <p 3 (x, t) , (2.2.2) 
we take the matrix elements of both sides 3 

h(pi " P2)M - Pif] e- i ^-P^( P2 | 0(0) | pi ) = 

= e -i(pi- M )^ {m 2 {p 2 \m\Pi) - A(p 2 |0 3 (O)bi)} , (2.2.3) 
3 Lorentz invariance imposes the relation (p2 \ <j)(x, t) \ pi ) = e~ ll " pl ~' P2 ^ lX ^ (p2 \ 0(0) |pi ) 



and we equate them to leading order 1/y/X. Assumption 3 implies that the kink momentum is 
very small compared to its mass, therefore in the left hand side of (j2.2.3j) we can neglect the 
energy difference 



(£i-£ 2 ) 2 =(^+...y=o(A 2 ) 



Hence the left hand side gives, at leading order, 



In the last term of the right hand side of (|2.2.3[) . the cubic power of (j> can be expanded over a 
complete set of states in the kink sector; in virtue of assumption 4, the leading term is obtained 
when the intermediate states are one-kink states: 

-\Y,(P2\m\p)(p\m\Q)(q\m\pi) = -a/ d ae ^-^ a [/>)] 3 . (2.2.4) 

p,q 

Hence, at leading order in A, the function /(a) obeys the same differential equation satisfied by 
the kink solution, i.e. 

^f(a) = X[f(a)f-m 2 f(a). (2.2.5) 

This means that we can take f(a) = (p c i(a), adjusting its boundary conditions by an appropriate 
choice for the value of the constant A in eq. Q2.1.4JI . 
Therefore, we finally obtain 

( p 2 \ (f)(0) \ Pl ) = J da e i{pi - p2)a (f) d (a) + higher order terms . (2.2.6) 

Along the same lines, it is easy to prove that the form factor of an operator expressible as a 
function of <j> is given by the Fourier transform of the same function of <p c i. For instance, the 
form factor of the energy density operator e can be computed performing the Fourier transform 

of£ ci(^) = |(^) 2 + V{cf> d ). 

Eq. (j2.2.6fl is the crucial achievement that opens the way to our study of non-integrable 
theories. However, it is worth to mention a further result obtained in [21], which confirms the 
interpretation of the stability frequencies w\ and uj q described in the previous Section. Let's 
consider the form factor of (j> between two states containing one kink and one kink plus one 
meson, respectively, and define a function fk(a) as 

( P 2,k\m\Pi) = fdae^-^ a f k (a). (2.2.7) 



Calling ivk the energy of the meson with momentum k, and equating again at leading order the 
two sides of the Heisenberg equation of motion evaluated between the above asymptotic states, 
we obtain the equation 

(4>ci) A (a) = <4f k (a) , (2.2.8) 



and consequently we can choose fk{°) = Vkip) I \F^ki where %(a) is one of the eigenfunctions 
of (|2.1.9f) and the constant \/yJ2ujk is chosen to realize a familiar boson normalization. This 
completely justifies the identification of the eigenvalues of (j2.1.9|) with the mesons energies. 



Similarly, if we define the form factor of <p between a kink in the ground state and one in the 
internal excited state as 

(P*2 I 0(0) | Pi ) = J da e^-v,)a ) (2 2 Q) 

and we call oj\ the energy of the internal excitation, we obtain 

f(a) = u\f*{a) , (2.2.10) 
which leads to f*(a) = rji{a) / ^J2oj\. 

2.2.2 Relativistic formulation of Goldstone and Jackiw's result 

The remarkable Goldstone and Jackiw's result has, unfortunately, two serious drawbacks: ex- 
pressing the form factor as a function of the difference of momenta, Lorentz covariance is lost 
and moreover, the antisymmetry under the interchange of momenta makes problematic any at- 
tempt to go in the crossed channel and obtain the matrix element between the vacuum and a 
kink-antikink state. 

In order to overcome these problems, in [Q we have refined the method proposed in |24) . 
This was done by using, instead of the space-momenta of the kinks, their rapidity variable 0, 
defined in (|1.2.3j) . The approximation of large kink mass used by Goldstone and Jackiw can be 
realized considering the rapidity as very small. For example, in the 4 theory (|2.1.5|) . where the 
kink mass M is of order 1/A, we work under the hypothesis that 8 is of order A. In this way we 
get consistently 

E = M cosh 6> ~ M , p = Msinh(9 ~ M6 <C M . (2.2.11) 

It is easy to see that the proof of (|2.2.6f) outlined in Sect. l2.2TT1 still holds, if we define the form 
factor between kink states as the Fourier transform with respect to the rapidity difference: 

( Pi I 0(0) |p 2 ) = f(0) =M J dae iM9a <p d (a) , (2.2.12) 

with the inverse Fourier transform defined as 

Ma) = J ^e- tMea f(6). (2.2.13) 

In fact, the analysis of the r.h.s. of (|2.2..3|) remains unchanged, while the one of the l.h.s. is even 
improved, since we can exploit the approximation (|2.2. 1 lj) directly on the covariant expression 

-(Pi -P2)m(pi -P2Y = -2M 2 (1 - cosh0) ~ M 2 9 2 . 

Our formalism only displays a small inessential difference with respect to |24j . in the orders in 
A of the various form factors (see assumption 4 in Sect. l2.2T|) , In fact, the M factor in front 
of (|2.2.12|) . which is a natural consequence of considering the rapidity as the basic variable, 
increases by 1/A the order of all form factors, without however altering any conclusion. 

With the above considerations, it is now possible to express the crossed channel form factor 
through the variable transformation 9 — > in — 8. We then have 



F 2 (9) = (O|0(O)|pi,p2> = /(i7r-0) 



(2.2.14) 
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Figure 2.3: Pictorial representation of (a) the form factor /(#), (b) the crossed channel form 
factor F2(8), (c) the form factor ^(0*) at the dynamical pole 9*. 



The analysis of this quantity allows us, in particular, to get information about the spectrum 
of the theory. As we have shown in (|1.2.16j) . in fact, its dynamical poles located at 9* = i(ir — u), 
with < u < 7T, correspond to kink-antikink bound states with masses given by 1)1.2. 9JI . which 
in this case specializes to 

u 

m {b) = 2Msin- . (2.2.15) 

(see Fig. l2.l-j)) . It is worth noticing that this procedure for extracting the semiclassical bound 
states masses is remarkably simpler than the standard DHN method of quantizing the corre- 
sponding classical backgrounds, because in general these solutions depend also on time and have 
a much more complicated structure than the kink ones. Moreover, in non-integrable theories 
these backgrounds could even not exist as exact solutions of the field equations: this happens for 
example in the (ft 4 theory, where the DHN quantization has been performed on some approximate 
backgrounds |23| . 

Once the matrix elements ([2.2.14)1 are known, one can estimate the leading behaviour in A of 
the spectral function in a regime of the momenta dictated by our assumption of small kink 
rapidity. The delta functions in (jl.2.11j) make meaningful the use of our form factors, derived 
in the small 9 approximation, only if we consider a regime in which p° ~ M and p 1 <C M and, 
from now on, we will always understand this restriction. The leading 0(1/ A) contribution to 
the spectral function, denoted in the following by p^\p 2 ), is given by the trivial vacuum term 
plus the kink-antikink contribution: 



5(t 



pW(p 2 ) = 2vr^W)l(0|<A(0)|0)| 2 + " 



4 M 2 / 2vr 



2dl + i7T ~M 



2 



(2.2.16) 



where the range of integration of the above quantity is of order p 1 /M (note that, being p l /M <C 
1, the integral can be roughly estimated by evaluating I-F2I 2 at 9± = 0: this is what we will do 
in the following). 



2.3 The Sine— Gordon model in infinite volume 

The Sine-Gordon model, defined by the potential 

2 

TTl 

= ^2(1- cos ^) , (2.3.1) 



gives us the chance of testing the validity of the semiclassical approximation. In fact, the 
integrability of this model has led to the exact computation of its spectrum and form factors 
|16 | I18[ 1191 12f)j. which can be evaluated in the small-/? limit to verify the semiclassical results. 

This model has been object of intensive research efforts, since it plays a relevant role in 
several physical contexts, either as a classical non-linear system or as a quantum field theory. For 
instance, its non-linear equation of motion appears in the study of propagation of dislocations 
in crystals, self-induced transparency effects in non-linear optics and spin wave propagation 
in liquid 3 He (see [SJ for a review). As a quantum field theory, it plays a crucial role in the 
description of quantum spin chains, in virtue of the bosonization procedure, which permits to 
express fermionic fields as exponentials of bosonic ones. In particular, an exact correspondence 
has been established in between the sine-Gordon and the Thirring model, which describes 
massive fermions interacting through a four-fermion coupling. 

The exact solution of the model shows that the spectrum of fundamental excitations is 
constituted by a soliton, an antisoliton and a given number of neutral particles, called "breathers" 
and obtained as soliton-antisoliton bound states ^0]. The number of breathers depend on (3 as 
the integer part of — , where 

-< = t4w* (2 - 3 ' 2) 

is the renormalized coupling constant. In particular, f3 2 = Air represent a free point, above which 
no bound states exist and the theory is in a repulsive regime, while the regime of small (3, that 
we will explore semiclassically, is deeply attractive. 

The potential l|2.3.1J) displays infinite degenerate minima, and the soliton solutions of the 
equations of motion are given by configurations which interpolate between two neighbouring 
vacua. Focusing on the two minima at <j> = and <j> = these are connected by the classical 
soliton 

^(x) = - arctan ( e m ( x - x °A , (2.3.3) 



P 

which is solution of eq. (|2. 1.4|) with A = and has classical mass M c \ = 8 -p. The semiclassical 
quantization of this background has been performed in |23| : eq. ()2.1.9|) can be cast in the 
hypergeometric form by using the variable z = \{1 + tanhmx), and its solution is expressed as 



1 /i W 2 1 /i ul 2 

v ( x ) = -z*V (I- z)~*V F 

\ 

The corresponding spectrum is given by the discrete zero mode 

ujq = , with Vo( x ) 




2 coshmx 
and by the continuous part 

.2 2f-, , „2\ .„ ;+l , „ f m \ _ Jmqx ^1 ~ tann mx 



u) q = m (1 + q ) , with Vq( x ) = e * 



1 + iq 



For this integrable theory we can easily verify the interpretation of the continuous spectrum as 
representing the scattering states of the kink with neutral particles. In fact, the semiclassical 
phase shift 5(q), defined as 

r]Jx) — ► e i[9mx^5{g)] with j/ n = 2 arctan I = _ i j 9 + » (2 _ 3 _ 4) 

x^ioo q q — i 



correctly reproduces the small-/? limit of the exact one 



£exact(0) = -2 log S kb (9) 

where the kink-breather scattering amplitude ^Hj 



Skb 



sinh 9 + i cos ^ 
sinh 9 — i cos ^ 

lb 



is evaluated in the kink rest frame, in which the breather momentum can be expressed as 
mq = to sinh #. 

Similarly to the broken (j) 4 case discussed in Sect. 12. ll the semiclassical correction to the kink 
mass is given by 

M-M d = \Ys [ m ^ 1 + In ~ \An + m 2 ] - I dx[l- cos^aix)} , (2.3.5) 



with 

oo 

2 = to 2 /? 2 r dk ^ 3 g ^ 

87r y V^ 2 + m 2 

— oo 

Determining the discrete values q n and A; n in a large finite volume of size R with periodic bound- 
ary conditions and sending R — > oo, one finally obtains the semiclassical quantum correction to 
the mass of the kink 

8m m 

M = T - - . (2.3.7) 

As we have already mentioned, it is known that the coupling constant renormalises as j3 2 — ► 7, 
and the exact quantum mass of the soliton is given by M = — , which coincides with the above 
expression (|2.3.7ft . The equality of the semiclassical and the exact results for the masses is a 
remarkable property of the SG model. 

In pQ we have derived and discussed the semiclassical form factor 1)2.2. 12JI relative to the 
kink background ()2.3.3)) . Since the first quantum corrections to M c i = 8^7 are of higher order 
in (3 2 , we consistently approximate the kink mass M with this value and we assume that the 
(anti) soliton rapidity is of order j3 2 . The semi-classical form factor 1)2.2.12)1 is explicitly given 
by 

2™ 1 2vr 2 / 6» \ , 

P 9 cosh j£9 P \P J 

In order to compare the corresponding F-2,{6) given by 1)2.2.14)1 with the semiclassical limit of 
the exact one, the only thing to take into account is that, in the definition of the exact form 
factor of the fundamental field <j)(x), the asymptotic two-particle state is actually given by the 
antisymmetric combination of soliton and antisoliton. Since at our level the form factor between 
antisoliton states is simply 

00 

< Pi I 0(O)|p 2 >= ^j- J dae^ M9 ^ a arctan {e~ ma ) = f(-9) , 
—00 

we finally obtain 

_ ... 4iri 1 



In the regime m — 9 ~ 0(/3 2 ), this function indeed agrees with the exact result JH] 



m = 2 4 
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(2.3.10) 



where £ = 7/87r and 



G(0) = exp 



eft sinhf(l-£) sin 2 |^ 
t sinh | £ cosh § sinh i 
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(2.3.11) 



Furthermore, we can also check that the dynamical poles of (|2,3.9|) . located at 
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2 <n< "2 + ^' 



consistently reproduce the odd part of the well-known semiclassical breathers spectrum (2 

(2n+ P 2 



m (2n+i) = 2Msin 



/3 2 

— (2n + l) 
16 v y 



(2n + l)m 



1 



3 x 8 3 



(2.3.12) 



Since in the vacuum sector ( | 4> | ) 
spectral function takes the form: 

p{p 2 ) = 4tt 3 5 



0, in this model the 1//3 2 leading contribution to the 



M 



P 2 (p 1 ) 2 cosh 2 



7T F_ 

2 m 



(2.3.13) 



Furthermore, in [H] we have presented a more quantitative comparison which permits to con- 
clude that formula (|2.2.12j) , though proven under the semiclassical assumption of small coupling 
and small rapidities, remarkably extends its validity to finite values of the coupling and to a 
large range of the rapidities. Consider, for instance, the form factor of the energy operator (up 
to a normalization N) g{9) = N(02 \ e(0) \6\), whose semiclassical and exact expressions are 
given, respectively, by 



5scmicl. (0) 



ffexact(#) = sinh- 



sinh 



Itt 

w 



sinh ^ 



G{6) 



(2.3.14) 
(2.3.15) 



with G(6) defined in ()2.3.11j) . Fig. l2.4l shows how, for small values of the coupling, the agreement 
between the two functions is very precise for the whole range of the rapidity. Furthermore, the 
discrepancy between exact and semiclassical formulas at larger values of /3 can be simply cured, 
in our example, by substituting the bare coupling (3 2 with its renormalized expression 7 into the 
semiclassical result ([2.3. 14(1 . as shown in Fig. 12.51 Hence we can conclude that the monodromy 
factor ([2.3.11[) . which is the relevant quantity missing in our approximation, does not play a 
significant role in the quantitative evaluation of the form factor even for certain finite values of 
the coupling 4 . 



4 It is easy to understand the reason of this conclusion in the above example: at small values of 9 we have 
G{9) ~ 1, whereas for 9 — > 00, when G(9) may contribute, the whole form factor goes anyway to zero. Similar 
conclusion can be reached for all other form factors which vanish at 9 — * 00. 
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Figure 2.4: Comparison between the exact function F given by 1)2.3.15(1 (continuous line) and 
its semiclassical approximation 1)2.3.14(1 (dotted line), at /? = 0.1 and (3 = 0.5. 




Figure 2.5: Comparison, at f3 = 1, between (a) the exact function g given by 1)2.3.15(1 (continuous 
line) and its semiclassical approximation 1)2.3.14)1 (dotted line), (b) the exact function g given 
by 1)2.3.15(1 (continuous line) and its semiclassical approximation 1)2.3.14(1 with the substitution 
(3 2 -> 7 (dotted line). 



Chapter 3 



Non-integrable quantum field 
theories 

As a natural development of the studies on integrable quantum field theories, there has been 
recently an increasing interest in studying the properties of non-integrable quantum field theories 
in (1 + 1) dimensions, both for theoretical reasons and their application to several condensed- 
matter or statistical systems. However, contrary to the integrable models, many features of 
these quantum field theories are still poorly understood: in most of the cases, in fact, their 
analysis is only qualitative and even some of their basic data, such as the mass spectrum, are 
often not easily available. Although one could always rely on numerical methods to shed some 
light on their properties, it is obviously important to develop and apply some theoretical tools 
to control them analytically. In this respect, there has been recently some progress, thanks to 
two different and complementary approaches. 

The first approach, called the Form Factor Perturbation Theory (FFPT) [31], is best suited 
to deal with those non-integrable theories close to the integrable ones. The second approach, 
given by the semiclassical method introduced in Chapter is on the other hand best suited to 
deal with those quantum field theories having kink excitations of large mass in their semiclassical 
limit. Although this method is restricted to work in a semiclassical regime, it permits however 
to analyse non-integrable theories in the whole coupling-constants space, even far from the 
integrable points. 

In this Chapter, after illustrating the FFPT technique (Sect . 1X1]) . we describe the original 
results obtained with the semiclassical method for the </> 4 field theory in the broken symmetry 
phase (Sect . HO]! and for the Double Sine-Gordon model (Sect. EO|l . The Appendices are devoted 
to some technical results obtained in the study of the Double Sine-Gordon theory. 

3.1 Form Factor Perturbation Theory 

The Form Factor Perturbation Theory technique (FFPT) , introduced in [HI] , permits to obtain 
quantitative predictions on the mass spectrum, scattering amplitudes and other physical quan- 
tities in non-integrable theories which can be seen as a perturbation of integrable ones. As any 
other perturbation scheme, it works finely as far as the non-integrable theory is an adiabatic 
deformation of the original integrable model, i.e. when the two theories are isospectral. This 
happens when the field which breaks the integrability is local with respect to the operator which 
creates the particles. If, on the contrary, the field which moves the theory away from integra- 
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bility is non-local with respect to the particles, the resulting non-integrable model generally 
displays confinement phenomena and, in this case, some caution has to be taken in interpreting 
these perturbative results. 

The method can be applied when the action A of the theory in exam is represented by a 
deformation of an integrable one .Ao through a given operator ^f: 

A = A + gj d 2 x^{x) . (3.1.1) 

One of the first consequences of moving away from integrability is a change in the spectrum of 
the theory: the first order corrections to the mass of the particle a belonging to the spectrum 
of the unperturbed theory is in fact given by 

8ml = 2gF^(i7r) + 0(g 2 ), (3.1.2) 

where the particle-antiparticle form factor of the operator ^(x), defined by the matrix element 

F%(0 1 -0 2 ) = (0 | *(0) | a(0i)a(0 2 )> , 

is introduced as function of the difference of the rapidity variables, defined in (|1.2.3|) , The mass 
correction (|3.1,2j) may be finite or divergent, depending on the locality properties of the operator 
ty(x) with respect to the particle a. The situation was clarified in [201 and it is worth recalling 
the main conclusion of that analysis. 

As we have seen in (|1,2.14|) . (|1.2,15|) . if the operator which interpolate particle a is non-local 
with respect to the perturbing operator vl/, the form factor F^ a (9) displays an annihilation pole 
at 9 = in, which is exactly the value at which the form factor is evaluated in (|3.1.2|) . Therefore, 
in this case the mass correction received by particle a under the perturbation ^ is formally 
infinite, unless the residue in (|1.2.15|) vanishes. 

Let's consider as an example the case in which the unperturbed theory is the sine-Gordon 
model (|2.3.1|) with frequency (3, and the perturbing operator is the exponential ^ a = e taip . 
This vertex operator is semi-local with respect to the soliton particle s, with the index 7 QS 
entering (jl.2.14|) given by -y a ^ s = a//3, whereas it is local with respect to the breather particles 
b (i.e. 7 Qi 6 = 0). This implies that formula (|3. 1.2(1 can be safely applied to compute the 
first order correction to the mass of the breathers, whereas a divergence may appear in an 
analogous computation of the mass correction of the solitons. This divergence has to be seen as 
the mathematical signal that the solitons of the original integrable model no longer survive as 
asymptotic particles of the perturbed theory, i.e. they are confined. 

3.2 Broken 4 theory 

In this Section, we will describe the original results obtained in and [2] about the spectrum 
of the (ft 4 field theory in the Z 2 broken symmetry phase. This non-integrable theory, defined by 
the potential (|2.1.5j) 

A 4 m 2 2 m 4 

V ^ = -^ 4 -ir^ + 4A' 

is invariably referred to as a paradigm for a wealth of physical phenomena. Just to mention 
a relevant example, in the Landau-Ginzburg language discussed in Sect. 11.31 it describes the 
universality class of the Ising model at low temperature. In spite of this deep interest, however, 
its non-perturbative features are still poorly understood. 



3.2.1 Semiclassical spectrum 



The main properties of the potential (J2.1.5j) and its kink background (|2.1.fi|) have been already 
discussed in Sect. 12. ll where the semiclassical quantization of the kink has been described. Here 
we directly move to the computation of the kink-antikink form factors, in order to extract from 
their analytical structure information about the complete spectrum of the theory. 

With our formulation in terms of the rapidity, we can write an unambiguous 1 Lorentz co- 
variant expression for the form factor (|2.2.12|) 
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(3.2.1) 



where the kink mass is approximated at leading order by the classical energy M = ^p^y-- 

It is of great interest to analyse in this case the dynamical poles of i^(^) f° r extracting 
information about the spectrum. They are located at 
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and the corresponding bound states masses are given by 
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(3.2.2) 



(3.2.3) 



Note that the leading term is consistently given by multiples of y/2m, which is the known mass 
of the elementary boson of this theory 2 . Contrary to the Sine-Gordon model, we now have all 
integer multiples of this mass and not only the odd ones: this is because we are in the broken 
phase of the theory, where the invariance under <f> — > —(ft is lost. Furthermore, this spectrum 
exactly coincides with the one derived in [2S] by building approximate classical solutions to 
represent the "breathers". 

Another important information can be extracted from the residue of F2{0) on the pole 
corresponding to the lightest bound state . This quantity, indeed, has to be proportional to 
the one-particle form factor (0| <j>\ ) through the semiclassical 3-particle on-shell coupling of 
kink, antikink and elementary boson T fc ^ fe , as shown in (|1.2.16|) . Since in our normalization the 
one-particle form factor takes the constant value l/-\/2, at leading order in the coupling we get 



kkb 



(3.2.4) 



a quantity so far unknown in the non-integrable </> 4 theory. 

Finally, the 1/A leading contribution of the spectral function is given in this case by 



P(P 2 ) = y^ {P°/™) $ {P 1 /™) + ^ 
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(3.2.5) 



1 This has to be contrasted with the tentative covariant extrapolation discussed in 12 II. 
2 The elementary bosons represent the excitations over the vacua, i.e. the constant backgrounds <j>± 
therefore their square mass is given by V"(cf>±) = 2 m 2 . 



3.2.2 Resonances 



The appearance of resonances in the classical kink-antikink scattering has been studied for this 
theory with numerical techniques in In this work, the key ingredient for the presence of 
resonances was identified in the presence of the so-called "shape mode", which is the discrete 
eigenvalue uj\ of the small oscillations around the kink background, given in ()2.1.13|) . corre- 
sponding to the internal excitation of the kink. We have shown in |5j how this mechanism 
can be analytically explained in our formalism. In fact, Goldstone and Jackiw's result on the 
form factor of the field <j) between asymptotic states containing a simple and an excited kink, 
described by formulas l|2,2.9j) and (|2.2.10j) . can be also refined in terms of the rapidity variable. 
In our case, where the fluctuation iji is given in (|2.1.13j) . the covariant form factor, analytically 
continued in the crossed channel, is expressed as 
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The dynamical poles of this object correspond to bound states with masses 
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have masses in the range 



2M < m$ < 2M + ui 



(3.2.8) 
(3.2.9) 



and, therefore, they can be seen as resonances in the kink-antikink scattering. 

Since the numerical analysis done in is independent of the coupling constant 3 , a quan- 
titative comparison with our semiclassical result is rather difficult, due to the dependence on A 
of (|3.2.8|) . However, the presence of many resonance states seen at classical level is qualitatively 
confirmed to persist also in the quantum field theory at small A, i.e. in its semiclassical regime, 
according to (|3.2.8j) . 



3.3 Double Sine— Gordon model 

An interesting non-integrable model where both FFPT and Semiclassical Method can be used 
is the so-called Double Sine-Gordon Model (DSG). It is defined by the potential 

V(<p) = -^cos/3^- ± cos (£<p + 6\ +C, (3.3.1) 

where C is a constant that has be chosen such that to have a vanishing potential energy of the 
vacuum state. The classical dynamics of this model has been extensively studied in the past by 
means of both analytical and numerical techniques (see |SZj for a complete list of the results), 
while its thermodynamics has been studied in [HE] by using the transfer integral method |39j . 
and with the path integral technique (see for instance |40| ) . 

3 Classically, in fact, one can always rescale the field and eliminate the coupling constant A. 



With A or \i equal to zero, the DSG reduces to the ordinary integrable Sine-Gordon (SG) 
model with frequency (3 or (3/2 respectively. Hence the DSG model with a small value of one 
of the couplings can be regarded as a deformation of the corresponding SG model and studied, 
therefore, by means of the FFPT On the other hand, for (3 — > 0, irrespectively of the value 
of the coupling constants A and (i, the DSG model reduces to its semiclassical limit. Despite the 
non-integrable nature of the DSG model, its classical kink solutions are - remarkably enough 
- explicitly known [371 138] and therefore the Semiclassical Method can be successfully applied 
to recover the (semi-classical) spectrum of the theory. As we will see in the following, the two 
approaches turn out to be complementary in certain regions of the coupling constants, i.e. both 
are needed in order to get the whole mass spectrum of the theory, whereas in other regions they 
provide the same picture about the spectrum of the excitations. 

Apart from the theoretical interest in testing the efficiency of the two methods on this specific 
model where both are applicable, the study of the DSG is particularly important since this model 
plays a relevant role in several physical contexts, either as a classical non-linear system or as a 
quantum field theory. At the classical level, its non-linear equation of motion can be used in 
fact to study ultra-short optical pulses in resonant degenerate medium or texture dynamics in 
He 3 (see, for instance, |1J and references therein). As a quantum field theory, depending on 
the values of the parameters A, //, (3, S in its Lagrangian, it displays a variety of physical effects, 
such as the decay of a false vacuum or the occurrence of a phase transition, the confinement of 
the kinks or the presence of resonances due to unstable bound states of excited kink-antikink 
states. Moreover, it finds interesting applications in the study of several systems, such as the 
massive Schwinger field theory or the Ashkin-Teller model [SHI , as well as in the analysis of 
the 0(3) non-linear sigma model with 6 term |42| . i.e. the quantum field theory relevant for 
understanding the dynamics of quantum spin chains [43] 144] . The DSG model also matters in the 
investigation of other interesting condensed matter phenomena, such as the soliton confinement 
of spin-Peierls antiferromagnets [45] . the dynamics of the spin chains in a staggered external 
field or the electron interaction in a staggered potential [4fij . 

Motivated by the above combined theoretical and physical interests, we have performed in [2] 
a thorough study of the spectrum of the DSG model. In some regions of the coupling constants, 
the obtained original semiclassical results have revealed to be necessary to complete the spectrum 
analysis, while in other regions they have been compared with the analogous ones produced by 
FFPT. Here we present the main result of the semiclassical analysis of the spectrum of the DSG 
model and its comparison with the results coming from FFPT fSect. l3~3~Tj) . Furthermore, the 
phenomenon of false vacuum decay will be also discussed fSect. l3~3.2|) . Finally, this Chapter's 
appendices are devoted to some technical results obtained in this study. In particular, Appendix 
!3IAI presents the computation of the kink mass corrections by using the FFPT, Appendix l3IBI lists 
the relevant expressions of the semiclassical form factors in DSG model, Appendix I3K 'I discusses 
the analysis of neutral states in comparison with the Sine-Gordon model, and Appendix I3IDI 
describes the basic results in a closely related model, i.e. the Double Sinh-Gordon model. 

3.3.1 Semiclassical spectrum 

We will study the theory defined by ([3.3.1]) in a regime of small f3, where the semiclassical results 
are expected to give a valuable approximation of the spectrum 4 . At the quantum level, the 

4 By applying the stability conditions found in 1301 to this model, they reduce to the condition /3 2 < 87r. Hence, 
for these values of /3 and, in particular in the semiclassical limit j3 — » 0, the potential 13.3.11 is stable under 



different Renormalization Group trajectories originating from the gaussian fixed point described 
by the kinetic term ^(d^(p) 2 of the lagrangian are labelled by the dimensionless scaling variable 
rj = A/i~( 8?r ~^ / 4 )/( 87r_ / 3 ) which simply reduces to the ratio rj = ^ in the semiclassical limit. 
When A or fx are equal to zero, the DSG model coincides with an ordinary Sine-Gordon model 
with coupling (3 or /3/2, and mass scale y/Ji or \J A/4, respectively. 

Since for general values of the couplings the potential (|3.3.1j) presents a -^-periodicity, it 
was noticed in |17j that one has an adiabatic perturbation of an integrable model only if the 
A = theory is regarded as a two-folded Sine-Gordon model. This theory is a modification of 
the standard Sine-Gordon model, where the period of the field (f) is defined to be -£-, instead of 
^ [IS]- As a consequence of this new periodicity assignment, such a theory has two different 
degenerate vacua \k ), with k = 0, 1 and \k + 2 ) = \k), which are defined by (k\<p\k) = jf k. 
Hence it has two different kinks, related to the classical backgrounds by the formula 

Ikiv A 

Kt\ k+l {x) = — + ~p ^ctan e mx k = 0, 1 , (3.3.2) 
and two corresponding antikinks, related to the classical solutions by the expression 
K c k l + i tk (x) = ^ + |arctan e~ mx 

= ^±1)1 _ l arctan e mx k = 0, 1 . (3.3.3) 

Finally, in the spectrum there are also two sets of kink-antikink bound states b$ , with I = 0, 1 
and n = 1,..., ^ . 

The flow between the two limiting Sine-Gordon models (with frequency (5 or (3/2, respec- 
tively) displays a variety of different qualitative features, including confinement and phase tran- 
sition phenomena, depending on the signs of A and /j,, and on the value of the relative phase 5. 
However, it was observed in |36| that the only values of 8 which lead to inequivalent theories are 
those given by \S\ < \. Furthermore, in virtue of the relations 

Vs(^+^,X,fl) = V 5+7r / 2 (0, A, -/i) , 

Vs(-(j),X,n) = V-s(4>, A,//) , 

we can describe all the inequivalent possibilities keeping \x positive and the relative phase in 
the range < 5 < 5. The sign of the coupling A, instead, simply corresponds to a shift or a 
reflection of the potential, without changing its qualitative features. As we are going to show in 
the following, the case 6 = 5 displays peculiar features, while a common description is possible 
for any other value of 6 in the range < 5 < ? . 

It is worth mentioning that the possibility of writing exact classical solutions for all the 
different kinds of topological objects in this model finds a deep explanation in the relation 
between the trigonometric potential (|3.3.1|) and power-like potentials. In fact, defining 

7171 4 

if = — ± -arctany , n = 0,1,2,3, (3.3.4) 

P P 



one can easily see that the first order equation which determines the kink solution 
renormalization and no countertems have to be added. 



is mapped into the equation for Y 

i(f)W), (3.3.5) 

where U(Y) describes various kinds of algebraic potentials, depending on the values of n, 5 and 

C. The 5 = case was analysed in (32] and its classical solutions are very simple because U(Y) 

only contains quartic and quadratic powers of Y. It is easy to see that a similar situation also 

i ( .. , A 2 



occurs in the 5 = ^ case; for instance, choosing n = 1 and C = — ^ + ) ) one obtains the 
quartic potential 

which has the well known classical background 



For generic 5, instead, also cubic and linear powers of Y appear, making more complicated the 
analysis of the classical solutions. 



5 = case 

It is convenient to start our discussion with the case 5 = 0. This case, in fact, displays those 
topological features which are common to all other models with < 5 < f, but it admits a 
simpler technical analysis, due to the fact that parity invariance survives the deformation of 
the original SG model. As we will see explicitly, in this case the results of the FFPT and 
the Semiclassical Method are complementary, since they describe different kinds of excitations 
present in the theory. 




Figure 3.1: DSG potential in the case 5 = 0. 

Fig l3.1l shows the shape of this DSG potential in the two different regimes, i.e. (i) < A < 4/x 
and (ii) A > The absolute minimum persists in the position (mod^) for any values of 
the couplings, while the other minimum at jj- (mod -£) becomes relative and disappears at the 
point A = 4fi. The breaking of the degeneration between the two initial vacua in the two-folded 
SG causes the confinement of the original SG solitons, as it can be explicitly checked by applying 
the FFPT. The linearly rising potential, responsible for the confinement of the SG solitons, gives 



rise then to a discrete spectrum of bound states whose mass is beyond 2M$g, where Msg is the 
mass of the SG solitons |36[ I45j . 

The disappearing of the initial solitons represents, of course, a drastic change in the topo- 
logical features of the spectrum. At the same time, however, a stable new static kink solution 
appears for A ^ 0, interpolating between the new vacua at and The existence of this new 
topological solution is at the origin of the complementarity between the FFPT and the Semi- 
classical Method. By the first technique, in fact, one can follow adiabatically the deformation 
of the SG breathers masses: these are neutral objects that persist in the theory although the 
confinement of the original kinks has taken place. It is of course impossible to see these particle 
states by using the Semiclassical Method, since the corresponding solitons, which originate these 
breathers as their bound states, have disappeared. Semiclassical Method can instead estimate 
the masses of other neutral particles, i.e. those which appear as bound states of the new stable 
kink present in the deformed theory. 

This new kink solution, interpolating between and ^r, is given explicitly by 
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(3.3.8) 



where 
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(3.3.9) 



is the curvature of the absolute minimum. Interestingly enough this background admits 
an equivalent expression in terms of the superposition of two solitons of the unperturbed Sine- 
Gordon model, centered at the fixed points ±R 

cp K (x) = (psc(x + R) + ¥sg(x ~ R) , 

where <psg( x ) = Jf arctan [e mx ] are the usual Sine-Gordon solitons with the deformed mass 
parameter (|3.3.9|) whereas their distance 2R is expressed in terms of the couplings by 



R = — arccosfH/-^- + 1 
m V A 



By looking at Fig l3.21 it is clear that this background, in the small A limit, describes the two 
confined solitons of SG, which become free in the A = point, i.e. where R — > oo. 
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Figure 3.2: Kink solution ()3.3.8j) 



The classical energy of this kink is given by 
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and in the A — > limit it tends to twice the classical energy of the Sine-Gordon soliton, i.e 



(3.3.10) 
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therefore confirming the above picture. In the fj, — > limit, the asymptotic value of the above 
expression is instead the mass of the soliton in the Sine-Gordon model with coupling (3/2. The 
expansion for small [i 



8VA/4 fi 32 2 , 



(3.3.11) 
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gives the first order correction which is in agreement with the result of the FFPT in the semi- 
classical limit (see Appendix 131 A j) . 

The bound states created by the kink 1)3.3.8)) and its antikink can be obtained by looking at 
the poles of the semiclassical form factors of the fields <p(x) and e(x), reported in Appendix I3IB1 
and their mass are given by 5 
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For small \i we easily recognize the perturbation of the standard breathers in Sine-Gordon with 
P/2: 
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while the expansion of the bound states masses for small A 
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deserves further comments: in fact, although the above masses have well-defined asymptotic 
values, they do not correspond however to any state of the unperturbed SG theory. The reason 
is that the classical background (|3.3.8|) in the A — ► limit does not describe any longer a localized 
single particle. This implies that its Fourier transform cannot be interpreted as the two-particle 
form factor and, consequently, its poles cannot be associated to any bound states. 

A technical signal of the disappearing of the above mentioned bound states in the A — > 
limit can be found by computing the three particle coupling among the kink, the antikink and 
the lightest bound state, given by 1)1.2.16(1 . At leading order in (3 we get 



KKb 



l§y/2 m 



5 Due to parity invariance, the dynamical poles of the form factor of ip between kink states only give the bound 
states with n odd. The even states can be obtained from the form factor of the energy operator e(x). 



The divergence of the coupling as A — ► indicates that the considered scattering processes cannot 
be seen anymore as a bound state creation, i.e. the corresponding bound state disappears from 
the theory. A general discussion of the same qualitative phenomenon for the ordinary Sine- 
Gordon model can be found in [^0], where the disappearing from the theory of a heavy breather 
at specific values of (3 is explicitly related to the divergence or to the imaginary nature of the 
three particle coupling among this breather and two lightest ones. 

Summarizing, in this model we have three kinds of neutral objects, i.e. meson particles. The 
first kind (a) is given by the bound states originating from the confinement potential of the 
original solitons. These discrete states have masses above the threshold 2Msg-> where Msg is 
the mass of the SG solitons, and merge in the continuum spectrum of the non-confined solitons 
in the A — > limit )361 145j . The second kind (b) is represented by the deformations of SG 
breathers, that can be followed by means of the FFPT and have masses, for small A, in the 
range [0, 2M$g]- Finally, the third kind (c) is given by the bound states 1)3.3. 12jl of the stable 
kink of the DSG theory and they have masses in the range [0,4Msg]- All these mass spectra 
are drawn in Fig. 13.31 
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Figure 3.3: Neutral states coming from: a) solitons confinement, b) deformations of SG 
breathers, c) bound states of the kink 1)3.3.8)) 

Obviously, due to the non-integrable nature of this quantum field theory not all these par- 
ticles belong to the stable part of its spectrum. Apart from a selection rule coming from the 
conservation of parity, decay processes are expected to be simply controlled by phase-space 
considerations, i.e. a heavier particle with mass Mh will decay in lighter particles of masses rrii 
satisfying the condition 

M h > rm . 

i 

Hence, to determine the stable particles of the theory, one has initially to identify the lightest 
mesons of odd and even parity with mass m* and m\ (m* < mi), respectively. Then, the 
stable particles of even parity are those with mass m below the threshold 2m* whereas the 
stable particles of odd parity are those with mass m < m*_ + m\. For instance, in the fj, — »• 
limit we know that the only stable mesons are those given by the particles (c), as confirmed 
by the expansion 1)3.3.13)1 . Hence, in this limit no one of the other neutral particles is present 
as asymptotic states. For the mesons of type (a), this can be easily understood since they are 
all above the threshold dictated by the lightest neutral particle. The situation is more subtle, 
instead, for the states (b). However, their absence in the theory with jjl — ► clearly indicates 
that at some particular value of A even the lightest of these objects acquires a mass above the 



threshold 2m|^, with given by 1)3.3.12)1 . Analogous analysis can be done for other values 

of the couplings so that the general conclusion is that most of the above neutral states are 
nothing else but resonances of the DSG model. As it is well known, unstable particles should 
have a non-zero imaginary part in their masses, as we have discussed in p. 2. 10(1 . but expression 
()3.3,12)) is real. This is due to the fact that, at semiclassical level, the fingerprint of instability is 
the imaginary nature of some of the frequencies Ui, eigenvalues of the stability equation (|2. 1.9(1 . 
which are not considered in the leading order expression ()3.3.12)1 for the bound states masses. 
Hence, although at the leading order we are missing the imaginary contributions to the masses, 
we always have to keep in mind that they come from some of the Wj. 

Another kind of resonances has been numerically observed in the classical scattering of kink 
and antikinks of the type 1)3.3.8)1 . The corresponding study, performed in |37j . relies on the ideas 
developed in for the broken (j) 4 theory, and it is essentially based on the presence of the 
"shape mode". Along the lines of the discussion presented in Sect. 13. 2) this mechanism can be 
easily interpreted in our formalism. Unfortunately, in the case of DSG we were not able to solve 
analytically the stability equation around the kink background (|3.3.8|) . However, the conclusion 
that this kind of resonances is described by the bound states of a normal and an excited kink 
remains unchanged. 

In addition to the above scenario of kink states and bound state thereof, in the region A < 4/i 
there is another non-trivial static solution of the theory, defined over the false vacuum placed 
at ip = It interpolates between the two values ^ and ^ — j| arccos(l — A/2/x), and then it 
comes back (see Fig. l3,4[) . Its explicit expression is given by 
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(3.3.14) 



m} = f i-' 1 (3.3.15) 

is the curvature of the relative minimum. Similarly to the kink (|3.3.8|) . it admits an expression 
in terms of a soliton and an antisoliton of the unperturbed SG model: 

<Pb(x) = <p SG (x + R)+ Vsg(-(x - R)) , (3.3.16) 

where now <psg( x ) = -* arctan [e mf x ] are the Sine-Gordon solitons with the deformed mass 
parameter 1)3.3.15)1 whereas their distance 2R is now given by 



R = arcsinhA/-^- — 1 . (3.3.17) 

rrif V A 

In the small A limit, it is clear that this background describes the confined soliton and 
antisoliton of the SG model, which become free in the A = point, i.e. where R — ► oo. 

The classical background 1)3.3. 14|) is not related to any stable particle in the quantum theory. 
This can be directly seen from equation 1)2.1.9)1 : in fact, Lorentz invariance always implies the 
presence of the eigenvalue oJq = 0, with corresponding eigenfunction rjo(x) = ^fci(x)- However, 
in the case of the solution ()3.3.14)) the eigenfunction 770 clearly displays a node, which indicates 
that the corresponding eigenvalue is not the smallest in the spectrum. Hence, there must 
be a lower eigenvalue to 2 _ 1 < 0, with a corresponding imaginary part of the mass relative to 
this particle state. Furthermore, the instability of 1)3.3.14)1 can be related to the theory of false 
vacuum decay j^El]: due to the deep physical interest of this topic, we will discuss it separately 
in Section 13.3.21 





Comments on generic 5 case 

We have already anticipated that the qualitative features of the theory relative to 5 = case are 
common to all other theories associated to the values of 5 in the range < 5 < \. This can be 
clearly understood by looking at the shape of the potential, which is shown in Fig. 13. 51 for the 
case 5 = ^ . 

In contrast to the 5 = case, parity invariance is now lost in these models, and the minima 
move to values depending on the couplings. Furthermore, in addition to the change in the nature 
of the original vacuum at which becomes a relative minimum by switching on A, there is also 
a lowering of one of the two maxima. These features make much more complicated the explicit 
derivation of the classical solutions, as we have mentioned at the beginning of the Section. 

However, it is clear from Fig. lH.ftl that the excitations of these theories share the same nature 
of the ones in the 5 = case. In fact, the original SG solitons undergo a confinement, while a 
new stable topological kink appears, interpolating between the new degenerate minima. Hence, 
the analysis performed for 5 = still holds in its general aspects, i.e. also in these cases the 
spectrum consists of a kink, antikink, and three different kinds of neutral particles. 



5 = 5 case 

The value <5 = § describes the peculiar case in which no confinement phenomenon takes place, 
since the two different vacua of the original two-folded SG remain degenerate also in the per- 
turbed theory. As a consequence, the original SG solitons are also asymptotic states in the 



perturbed theory. By means of the semiclassical method we can then compute their bound 
states, which represent the deformations of the two sets of breathers in the original two-folded 
SG. Hence, in this specific case FFPT and semiclassical method describe the same objects, and 
their results can be compared in a regime where both j3 and A are small. 

V(» = l,\ = l) V(ji = l, A = 4) 
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Figure 3.6: DSG potential in the <5 = \ case. 



Fig. 13.61 shows the behavior of this DSG potential. There are two regions, qualitatively 
different, in the space of parameters, the first given by < A < 4/i and the second given by 
A > 4/i. They are separated by the value A = 4/i which has been identified in as a phase 
transition point. We will explain how this identification is confirmed in our formalism. 

Let's start our analysis from the coupling constant region where A < 4/i. Switching on A, the 
original inequivalent minima of the two-folded Sine-Gordon, located at cf) m i n = 0, (mod -£-), 



remain degenerate and move to 4> m i n 
common curvature of these minima is 
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(3.3.18) 



Correspondingly there are two different types of kinks, one called "large kink" and interpolating 
through the higher barrier between — <f)Q and ^ + (f>o, the other called "small kink" and inter- 
polating through the lower barrier between + 4>o and — 4>q. Their classical expressions are 
explicitly given by 
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(mod 4vr) , (3.3.19) 
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(3.3.20) 



With the notation previously introduced, the vacuum structure of the corresponding quantum 
field theory consists of two sets of inequivalent minima, denoted by | ) and | 1 ) , identified 
modulo 2, i.e. | a + 2n) =| a). The spontaneous breaking of the symmetry T : (p — > 2tt — <p 
selects one of these minima as the vacuum. If we choose to quantize the theory around | ), the 
admitted quantum kink states are | L) =| Kq^ ) and | S) =\ Kq-i), with the corresponding 
antikink states | L) =| K\ o ) and | S) =\ i^-i,o ), and topological charges 



Multi-kink states of this theory satisfy the selection rule coming from the continuity of vacuum 
indices and are generically given by 
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The leading contributions to the masses of the large and small kink are given by their classical 
energies, which can be easily computed 

8m ( A / 7r A \ 1 , 



The expansion of this formula for small A is given by 



and the first order correction in A coincides with the result of FFPT in the semiclassical limit 
(see Appendix I3IA|) . 

Since two different types of kink | L) and | S) are present in this theory, one must be careful 
in applying eq. (|2.2.12j) to recover the form factors of each kink separately. In fact, one could 
expect that both types of kink contribute to the expansion over intermediate states (|2.2.4j) used 
in [21] to derive the result. For instance, starting from the vacuum | 0} located at (fi m i n = — 4>q 
there might be the intermediate matrix elements q(S \ O \ L)q and q(L \ O \ S)o. However, if O is 
a non-charged local operator, it easy to see that these off-diagonal elements have to vanish for 
the different topological charges of | L) and | S). Hence, the expansion over intermediate states 
diagonalizes and one recovers again eq. (|2.2.12|) . 

Therefore, from the dynamical poles of the form factor of if on the large and small kink- 
antikink states, reported in Appendix I3IB1 we can extract the semiclassical masses of two sets 
of bound states: 

m$ = 2M L sin (n L ^\ , < n L < ir^± , (3.3.23) 

= 2M S sin (n s , < n s < vr^ . (3.3.24) 

Expanding for small A, we can see that these states represent the perturbation of the two sets 
of breathers in the original two-folded Sine-Gordon model: 
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+ 0(A 2 ) (3.3.25) 



A discussion of these results, in comparison with previous studies of this model [17]) is reported 
in Appendix I3ICI 

Concerning the stability of the above spectrum, for A < 4/i the only stable bound states 
are the ones with ^ < 2m|^; for A close enough to 4^z, however, the small kink creates no 

bound states, hence the stability condition 6 becomes < 2m|^ . 



(i 



Specifically, this stability condition holds for A > A*, where A* is defined by 



A* /tt . A*\ „ B 

arcsm 



^/I6fi 2 - A* 2 V 2 8tt 



Resonances in the classical scattering of the small kinks and antikinks (|3.3.2U|) have also been 
numerically observed in Our analysis is in agreement with these results, and it adds in this 
case another possibility. In fact, these resonances can be related both to the bound states of 
small excited kink-antikink with masses m 



g£' in the range 



2M S < mf) < 2M S + wi , 
and to the large kink-antikink bound states with masses in the range 

2M S < < 2M L , 

where are given by 1)3.3.231) . 

In the limit A — > 4/x, </>o tends to 5, the two minima at + <^o an d — </>o coincide and 
the small kink disappears, becoming a constant solution with zero classical energy. All the large 
kink bound states masses collapse to zero, and in this limit all dynamical poles of the large kink 
form factor disappear. This is nothing else but the semiclassical manifestation of the occurrence 
of the phase transition present in the DSG model (see [SB])- 

In the second coupling constant region, parameterized by A > 4//, there is only one minimum 
at fixed position — ^ (mod %), with curvature 
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There is now only one type of kink, given by 
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(3.3.26) 



(3.3.27) 



Its classical mass, expanded for small /x, is again in agreement with FFPT (see Appendix 131 Aj) : 

16m | „ A ( 7r .A — 8/A I 



M K 



0" 



7T .A 

1 H = arcsm — 

4y^(A - 4/x) V 2 

^37x + 0(M) 



a-o (/3/2) 2 

The bound states of this kink (see Appendix I3IBI for the explicit form factors) have masses 



(3.3.28) 
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For small fx, these states are nothing else but the perturbed breathers of the Sine-Gordon model 
with coupling (3/2: 
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In closing the discussion of the 5 = ir/2 case, it is interesting to mention another model 
which presents a similar phase transition phenomenon, although in a reverse order. This is the 
Double Sinh-Gordon Model (DShG), discussed in Appendix I3IUI The similarity is due to the 
fact that also in this case a topological excitation of the theory becomes massless at the phase 
transition point, but the phenomenon is reversed, because in DSG the small kink disappears 
when A reaches the critical value, while in DShG a topological excitation appears at some value 
of the perturbing coupling. 



3.3.2 False vacuum decay 



The semiclassical study of false vacuum decay in quantum field theory has been performed by 
Callan and Coleman |51j . in close analogy with the work of Langer |52j . The phenomenon occurs 
when the field theoretical potential U((p) displays a relative minimum at <y?+: this classical point 
corresponds to the false vacuum in the quantum theory, which decays through tunnelling effects 
into the true vacuum, associated with the absolute minimum (see Fig. I3.7JI . 
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Figure 3.7: Generic potential for a theory with a false vacuum 

The main result of [SJ is the following expression for the decay width per unit time and unit 
volume: 
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(3.3.30) 



det[-d 2 + U"(<p+)] 

specialized here to the case of two-dimensional space-time. Omitting any discussion of the 
determinant, about which we refer to the original papers [H5> we wm present here an explicit 
analysis of the coefficient B. 

It has been shown that B coincides with the Euclidean action of the so-called "bounce" 
background ifB- 

2 



(3.3.31) 



B = Se = 27r / dp p 



This classical solution is the field-theoretical generalization of the path of least resistance in 
quantum mechanical tunnelling; it only depends on the Euclidean radius p = Vt 2 + x 2 and 
satisfies the equation 

d 2 ip B , 1 dip B 



dp 2 
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with boundary conditions 
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(3.3.32) 
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Although in general one does not know explicitly the bounce solution, it is possible to set 
up some approximation to extract a closed expression for the coefficient B. The so-called "thin 
wall" approximation consists in viewing the potential U ((f) as a perturbation of another potential 
U + (ip), which displays degenerate vacua at ip± and a kink <pk(x) interpolating between them. 
The small parameter for the approximation is the energy difference e = U(<p+) — U((p—). 



In this framework, one can qualitatively guess that the bounce has a value <p(0) very close 
to </?_, then it remains in this position until some vary large p = R and finally it moves quickly 
towards the final value <p+. For p near R, the first-derivative term in eq. (|3.3.32[) can be 
neglected; if in addition one also approximates U with U+, then one can express the unknown 
bounce solution as 

if- p < R 

<Pb(p) = I pk(p-R) P~R (3.3.33) 

^<p+ p > R . 

Since the bounce has to represent the path of least resistance, the parameter R, free up to this 
point, can be fixed by minimizing the action 

S E = -TTR 2 e + 2irRM K , 

which is given by the sum of a volume term and a surface term. Hence, the condition = 
is realized by the balance of these two different terms in competition, and it finally gives 

- B-**&. (3.3.34) 

In the DSG model, however, we know explicitly the bounce background in the thin wall 
regime (here we have e = 1?), without any approximation on the potential. This is given by the 
solution (|3.3.14j) with x replaced by p, that can be directly used to estimate the decay width. 
Unfortunately the integral in (|3.3.31|) does not admit a simple expression to be expanded for 
small A, but it is clear from eq. (|3.3.16|) and Fig. 13.41 that the leading contribution is given by 

R+Ar 2 

S E c2,R J cfcpgS (R-x)] -^log(^), (3-3.35) 

R-Ar 

with R given by (|3. 3.17(1 . This behavior in A does not agree with the general prediction (|3.3.34j) . 
The reason can be traced out in the fact that eq. (|3. 3.16(1 explicitly realizes the relation between 
the bounce and the kink of the unperturbed theory, but in a more sophisticated way than 
(|3.3.33|) . In fact, the mass parameter nif of the SG kink ipsG is dressed to be the one of the 
DSG theory, and the parameter R is not free, since ([3.3. 14[) is already the result of a minimization 
process, being a solution of the Euler-Lagrange equations. The thin wall approximation can be 
still consistently used because R is very big for small A, while the crucial difference is that the 
volume term is now missing from the action, since the value <Pb(Q) = fi is the so-called classical 
turning point (see Fig. 13. 7|) . degenerate with the false vacuum. It is worth noting that the path 
of least resistance in quantum mechanics precisely interpolates between the false vacuum and 
the turning point. 

Up to the determinant factor, our result for the leading term in the decay width is then 

r 4 / A \ 8 " //32 , /16/A 



It will be interesting to investigate whether the above mentioned difference with the prediction 
(|3.3.34() is a particular feature of the DSG model or it appears for a generic potential if one 
improves the approximate description of the bounce along the lines discussed here. 



3.4 Summary 



In this Chapter we have shown how the semiclassical method is an efficient tool for studying 
non-integrable QFT. In fact, applied to the broken t/> 4 theory, it has provided new simple results 
about mass spectrum, three-particle couplings and resonances. 

In cases where also the Form Factor Perturbation Theory can be applied, the semiclassical 
method may be complementary to this technique or it may provide results comparable with 
it. We have applied both methods for analysing the mass spectrum of the non-integrable QFT 
given by the Double Sine-Gordon model, for few qualitatively different regions of its coupling- 
constants space. This model appears to be an ideal theoretical playground for understanding 
some of the relevant features of non-integrable models. By moving its coupling constants, in 
fact, it shows different types of kink excitations and confinement phenomena, a rich spectrum of 
meson particles, resonance states, false vacuum decay and the occurrence of a phase transition. 
In light of the many applications it finds in condensed matter systems, it would be interesting 
to investigate further its properties. 

Finally, it is worth mentioning that the semiclassical method can be easily extended to the 
analysis of several power-like potentials, which have interesting physical interpretations along 
the lines of the Landau-Ginzburg theory discussed in Sect. ll.3l 



Appendix 3. A. Kink mass corrections in the FFPT 



In this Appendix we compute by means of the FFPT the corrections to the kink masses in the 
semiclassical limit, which is relevant for a comparison with our results. 

For small A, we have to consider the DSG model as a perturbation of the two-folded Sine- 
Gordon In the 5 = tt/2 case, the perturbing operator is ^ = sin^0. Its form factors 
between the vacuum and the two possible kink-antikink asymptotic states are obtained at the 
semiclassical level by performing the Fourier transform of sin 
given by eq. 1)3.3. 2jl . Hence we obtain 
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The first order correction in A to the kink masses is then 
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in agreement with the correction to the classical masses (J3.3.22 
the large kink, and K 12 with the small one. 

In the 5 = case, instead, we can explicitly see how the solitons disappear from the spectrum 



as soon as A is switched on. The form factor of the operator \I/ 
at = m 
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The other interesting regime to explore is the small fx limit. In the case 5 = 0, this can 
be seen as the perturbation of the SG model at coupling (3 = (3/2 by means of the operator 
= cos 2j3(p. The semiclassical form factor is 

^aaW = tS-^( 1 - 2 ^)' ( 3A - 4 ) 
S p z sin nxy 

where we have defined y = -Jz(iTV — 9). The corresponding mass correction is given by 

in agreement with 1)3.3.11)1 . 

The case <5 = f can be described by shifting the original SG field as ip — > ip + 5. In this 
way the perturbing operator becomes — \? and we finally obtain the same mass correction but 
with opposite sign, as in ()3.3.28j) . 



Appendix 3.B. Semiclassical form factors 



In this Appendix we explicitly present the expressions of the two-particle form factors, on the 
asymptotic states given by the different kinks appearing in the DSG theory, of the operators 
<p(x) and e(x), the last one defined by 

These matrix elements are obtained by performing the Fourier transforms of the corresponding 
classical backgrounds, as indicated in (|2.2.12|) and 1)2.2.14)1 . We use the notation: 

F* R (9) = (0\*{0)\ K{9 X ) K{9 2 )} , 

with e = e 1 - e 2 . 



For the kink 1)3.3. 8JI in the 5 = case we have 
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while to and Mk are given by ()3.3.9)) and ()3.3.1U|) . respectively, and 
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where c = 1 + % ■ 



For the large kink 1)3.3.19)1 in the 8 = 5 case (with A < 4/x) we have 
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where c = — -A- . 



For the small kink (|3.3.20|) in the 8 = \ case (with A < 4/x) we have 
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where 
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while m and Ms are given by (|3. 3.18(1 and (|3.3.21|) . respectively, and 
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Finally, for the kink (|3.3.27|) in the 5 = ? case (with A > 4//) we have 
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while m and Ml are given by ((3.3.26(1 and (|3.3.28|) . respectively, and 
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where c = 1 — ^ . 



Appendix 3.C. Neutral states in the S = | case 



The semiclassical results reported in the text, i.e. eqs. ()3,3.23|) . (|3.3.24|) and (|3.3.25j) . pose an 
interesting question about the nature of neutral states in the DSG model at <5 = \. It should be 
noticed, in fact, that the first order correction in A obtained by the Semiclassical Method does 
not match with the results reported in where, by using the FFPT and an extrapolation of 
numerical data, the authors concluded that this correction was instead identically zero 7 . It is 
worth discussing this problem in more detail. 

In the standard Sine-Gordon model, the breathers | b n ), with n odd (or even), are defined 
as the bound states of odd (or even) combinations of | K K) and | K K) , where K represents 
the soliton and K the antisoliton. The combinations | K K ± K K) are eigenstates of the parity 
operator P : (f) —> —(f), which commutes with the hamiltonian and acts on the soliton trans- 
forming it into the antisoliton. The above mentioned identification of the bound states relies on 
a very peculiar feature of the Sine-Gordon S"-matrix in the soliton sector JH], whose elements 
are defined as 

K{e 1 )K{6 2 ) = S T (9 12 )K(9 2 )K(9 1 ) + S R (9 12 )K(9 2 )K(9 1 ) , 

K(e 1 )K(e 2 ) = s{e 12 )K{e 2 )K{e 1 ) , (3.C.1) 

K(e 1 )K(6 2 ) = S(9x2)K{B 2 )K{Bx). 

In fact, both the transmission and the reflection amplitudes St{9) and S R {9) display poles at 
#* = i(ir — n£), with residua which are equal or opposite in sign depending whether n is odd or 
even. Hence, the diagonal elements 

S-{9) = \[S T {e) - S R (9)] , (3.C.2) 

S+{9) = \[S T {e) + S R {9)] (3.C.3) 

have only the poles with odd or even n, respectively, and for each n there is only one bound 
state with definite parity. 

However, this is a special feature of the Sine-Gordon model which finds no counterpart, for 
instance, in other problems with a similar structure. As an explicit example, one can consider 
the (RSOS)3 scattering theory, which displays a 3-fold degenerate vacuum and two types of kink 
and antikink with the same mass. The central vacuum is surrounded by two other minima, as 
in the Sine-Gordon case, and this gives the possibility to define both a kink-antikink state and 
an antikink- kink state around it. The minimal scattering matrix, given in |53j . can be dressed 
with a CDD factor to generate bound states. It is easy to check that the common poles in 
the transmission and reflection amplitudes have in this case different residua, giving rise to two 
distinct bound states, degenerate in mass, over the central vacuum. 

Hence, if we call the bound states of kink-antikink and \bn^) the bound states of 

antikink-kink, in general we have to consider them as two distinct excitations, and if they have 
the same mass we can build two other states from their linear combinations 

\h {0) \ + \h {1) \ 

|6 W )= \bn)±_\b n ) _ (3C4) 



7 It is worth stressing that the linear correction ffi.3.25p in A is very small even for finite values of f3 (it is easy 



to check, indeed, that the first term of its expansion is (-jcf ) ) ant ^ somehow compatible with the numerical 
data given in (471 . 
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The peculiarity of the Sine-Gordon model is the removal of this double multiplicity due to the 
fact that the states | fe^n+i) anc ^ \^2n) decouple from the theory. This feature is shared also 
by the two-folded version of the model, since the kink scattering amplitudes have the same 
analytical form as in SG |48j . 

In the two-folded SG there are two different kink states | K-\fl) and | i£o,l) (see Sect. 13.31 
and ref. [IS] f° r the notation), and the parity P, which is still an exact symmetry of the theory, 
acts on them transforming the kink of one type into the antikink of the other type: 

P : | K 0tl ) -> | tf ,-i> , | K- lt0 ) | K lt0 ) . (3.C.5) 

If we quantize the theory around the vacuum | 0), we can define | bn^) as the bound states of 
I #0,1 Ki,o), and | b}P) 
transform under P as 



Kq\K\q) : and \bn^) as the bound states of | Kq _i K~.\ q). These degenerate states, which 



can be still organized in parity eigenstates | bn ) 5 and the particular dynamics of the problem 
causes the decoupling of half of them from the theory. Furthermore, it is easy to see that the 
form factors of an odd operator between two of these states has to vanish in virtue of the relation 

(Olsin^l&W&W) = (0 | P-'P (sin P^P\b^b^) = 

= -(Olsin^l&W&W), 

leading to the FFPT result that the breathers receive a zero mass correction at first order in A, 
as it is claimed in |47j . 

However, FFPT can be applied by taking into account the nature of neutral states in the 
DSG model, where the addition to the Lagrangian of the term sin ^tp spoils the invariance 
under P. The kinks | K-\q) and | -Ko,i) are deformed into the small and large kinks | S) and 
| L), respectively, which are not anymore degenerate in mass and cannot be superposed in linear 
combinations. Hence, the neutral states present in the theory are | bn ) an d | b^}, deformations 
of | bn^} and | b n ^) respectively. In virtue of the general considerations presented above, one can 
see that this interpretation does not lead to any drastic change in the spectrum. In fact, the 
states | &2n+i) an d I &2n ) have no reason to decouple in the DSG theory, but they have to carry 
a coupling which is a function of A adiabatically going to zero in the two-folded SG limit. 

A proper use of the FFPT on | bn) and | bn) reproduces indeed the situation described by 
l)3.3.25J) . in which the two sets of breathers receive mass corrections including also odd terms 
in A, but with opposite signs. This is easily seen by considering the P transformations in the 
two-folded SG model: 



( | sin 



§*|&SM?)) = (0 | P-'P (sin f^P-^lW) 



( 1 sin 2 Y , 



which gives, at first order in A, 



$m(L) = -$m {s) , (3.C.7) 

in agreement with our semiclassical result (|3.3.25|) . It is worth noting that also with this in- 
terpretation the total spectrum of the DSG model remains unchanged under the action of P, 



which corresponds to the transformation A — > —A. In fact, the two types of kinks and breathers 
are mapped one into the other. This is consistent with the observation that P, although it is 
not anymore a symmetry of the perturbed theory, simply realizes a reflection of the potential, 
hence the total spectrum should be invariant under it. 

Presently the above symmetry considerations seem to us the correct criterion to define the 
neutral states, and find confirmation in our semiclassical result (|3.3,25|) . However, the available 
numerical data presented in [17] pose a challenge to this interpretation and further studies are 
needed to solve this interesting and delicate problem. In fact, although Smt^ and 6m^ are 
not forced to vanish by symmetry arguments, there is in principle the possibility that both of 
them are identically zero in the complete quantum computation. This could follow from the use 
of the exact kink masses entering eqs. (|3.3.23j) and (|3.3.24j) . together with a proper shift of the 
semiclassical pole in the form factors, due to higher order contributions. The exact cancellation 
of the linear corrections is a very strong requirement, in support of which we have presently no 
indication in the theory, but a careful analysis of this point is neverthless an interesting open 
problem. 



Appendix 3.D. Double Sinh— Gordon model 



Among the different qualitative features taking place in perturbing integrable models, a situation 
particularly interesting is the one in which the perturbation is adiabatic for small values of 
the parameters but nevertheless a qualitative change in the spectrum occurs by increasing its 
intensity. 

This is indeed the situation in the 8 = ? case of DSG model, where we have two types of 
kinks for small A, but at A = 4/x one of them disappears from the spectrum. This phenomenon 
is obviously unaccessible by means of FFPT, hence the semiclassical method is the best tool to 
describe it. 

Here we consider another interesting example of this kind, realized by the Double Sinh- 
Gordon Model (DShG). In this case the phenomenon is even more evident, because in the 
unperturbed Sinh-Gordon model there are no kinks at all, but just one scalar particle, while 
perturbing it, at some critical value of the coupling a kink and antikink appear, i.e. there is a 
deconfinement phase transition of these particles. 

The DShG potential, shown in Fig. 13.81 is expressed as 
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Figure 3.8: DShG potential 

In the regime A < 4// the qualitative features are the same as in the unperturbed Sinh-Gordon 
model. At A = 4//, however, the single minimum splits in two degenerate minima, which for 
A > 4/z are located at <p± = i^arccosh^. A study of the classical thermodynamical properties 
of the theory in this regime has been performed in [23] with the transfer integral method. 

The kink interpolating between the two degenerate vacua is 
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Its classical mass is given by 
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From the form factor of ip on the kink-antikink asymptotic state, expressed as 
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we derive the bound states spectrum 
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All the kink-antikink bound states disappear from the theory at a certain value A* > 4/i such 

that tt^^- = 1, and the kink becomes a constant solution with zero classical energy when 

m x* 

A — ► 4//. This is the semiclassical manifestation of a phase transition, analogous to the one 
observed in DSG with 5 = ^. As we have already anticipated, here the phenomenon occurs in 
a reverse order, since in this case a kink appears in the theory by increasing the coupling A. 



Chapter 4 

Finite— size effects 



Quantum field theory on a finite volume is a subject of both theoretical and practical interest. 
It almost invariably enters the extrapolation procedure of numerical simulations, limited in 
general to rather small samples, but it is also intimately related to quantum field theory at 
finite temperature. It is therefore important to increase our ability in treating finite size effects 
by developing efficient analytic means. In the last years, a considerable progress has been 
registered in particular on the study of finite size behaviour of two dimensional systems. Also 
for these models, however, an exact treatment of their finite size effects has been obtained only 
in particular situations, namely when the systems are at criticality or if they correspond to 
integrable field theories. At criticality, in fact, methods of finite size scaling and Conformal 
Field Theory permit to determine many universal amplitudes and to extract as well useful 
information on the entire spectrum of the transfer matrix. Away from criticality, exact results 
can be obtained only for integrable theories which, on a finite volume, can be further analysed 
by means of Thermodynamical Bethe Ansatz. This technique provides integral equations for the 
energy levels, mostly solved numerically. In all other cases, the control of finite size effects in 
two dimensional QFT has been reached up to now either by perturbative or numerical methods. 

The semiclassical method is capable of shedding new light on this topic ^|31I3]. In fact, 
once the proper classical solutions are identified to describe a given geometry, the spectral 
function in finite volume can be easily estimated by adapting the Goldstone and Jackiw's result. 
Furthermore, the finite-volume kinks can be quantized semiclassically by adapting the DHN 
technique, which permits to write in analytic form the discrete energy levels as functions of the 
size of the system. This procedure can be also extended to geometries with boundaries. 

The Chapter is organized as follows. A brief overview of the known non-perturbative tech- 
niques to study finite-size effects is presented in Section l4~Tl Section l4~21 describes the semi- 
classical results obtained about form factors and correlation functions in finite volume, for the 
Sine-Gordon and broken <^> 4 theories on a cylinder with antiperiodic boundary conditions. Fi- 
nally, the semiclassical quantization of the Sine-Gordon model on a periodic cylinder and on 
a strip with Dirichlet boundary conditions are presented in Sections l4.3l and 14. 4( respectively. 
The Chapter also contains few appendices, devoted to some technical discussions. Appendix 
I4IAI presents the quantization of a free bosonic theory in a finite volume and a comparison of 
finite-volume and finite-temperature computations of the simplest one-point correlation func- 
tion. Appendix 1411:51 collects relevant mathematical properties of the elliptic functions used in 
the text whereas Appendix I4ICI displays the main properties of the Lame equation, which plays 
a central role in the finite-size quantization. 
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4.1 Known facts 



4.1.1 General results 



The understanding of finite-size effects is an important tool in the study of QFT. The finite- 
volume energy spectrum, in fact, contains a lot of information about the properties of the 
theory in infinite volume. This was first pointed out by Liischer |55j . who showed how the 
mass of a particle in a large but finite volume approaches exponentially its asymptotic value in 
a way controlled by the scattering data of the infinite-volume theory. Liischer's results were 
obtained in the context of Monte Carlo simulation of lattice field theories, with the intent of 
correctly interpreting the numerical data, unavoidably affected by finite-size effects. However, 
these findings proved to be useful for purposes other than merely controlling a systematic error 
source. In fact, they can be applied with the opposite spirit, for extracting scattering data from 
the numerical analyses. It is worth mentioning that this also provides an independent tool for 
checking the S*-matrices of integrable theories, often obtained on the basis of conjectures about 
the conserved quantities and particle content of the considered theories. 

We illustrate here Liischer's result for the mass corrections in finite- volume, restricting to 
the case of (1 + l)-dimensional space-time, with the space variable compactified on a periodic 
cylinder of circumference R. For the general d + 1-dimensional result and its proof, see |55M56| . 
The deviation of the mass of a particle from its value in infinite volume is shown to originate 
from polarization effects. In fact, the self-energy of the particle receives contributions from 
processes in which virtual particles "travel around the finite-size world" . The phenomenon can 
take place at leading order in two different ways, represented in Fig. 14.11 The first, possible 
only in theories with appropriate cubic couplings and indicated by (1), is a virtual process in 
which the particle splits in two constituents, which travel around the world before recombining 
to give back the original particle. The infinite-volume analog of this process is also depicted in 
Fig, l4.1l (T), The second process, indicated as (2) in Fig. l4,H is the analog of a tadpole diagram 
in infinite volume, and involves the scattering amplitudes of the given particle with the others 
in the theory, which travel around the world before annihilating again. 
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Figure 4.1: The two finite-volume contributions to the self-energy, and their infinite-volume 
analogs. 



Therefore, the finite-size mass of a particle a is given by 



m a (R) = m a (oo) + Am ^ (R) + Am^(R) + O (e~ aaR ) , 



(4.1.1) 



where the constant a a indicates the order of the errors, and 



Am«(i?) = -Y,S{m 2 a -\ml-m 2 c \)^ abc R abc e^ R , (4.1.2) 

b,c 

oo 

Am {2 \R) = V f d8 e- mbR cosh 9 m b cosh 6 (s%(9 + m/2) - l) (4.1.3) 

b -oo 

(P denotes the principal part of the integral). is the two-particle scattering amplitude 
defined in (jl.2.2)) and 9 is the rapidity variable introduced in 1)1.2. 3)) . In correspondence of a 
bound state pole at 9* = iu c ab , the S-matrix has the residue R a bc introduced in p.2.8j) . Finally, 
the quantity \i a \, c is defined as \i a bc = mb srnu aft' an< ^ the s * e P function 9 (m 2 a — \m 2 — m 2 \) 
restricts the sum in 1)4. 1.2)1 to those processes in which the virtual particles b and c can effectively 
recombine into a in the forward direction. It is worth noticing that the term Am^(R) is of 
order O (^e~ fJ-abcR ^ , while the term Am^ (R) is of order O (^ e ~ mbR ^ . 

As we have already commented, Liischer result has been derived in the context of numerical 
simulations in lattice QFT. Another efficient numerical method for studying two-dimensional 
QFT in finite-size, which does not require a discretization of space, was suggested in |HZj, and is 
called Truncated Conformal Space Approach (TCSA). This method is based on the possibility 
of viewing a massive QFT as a perturbation of some CFT, and its description gives us the 
possibility of discussing some general properties of finite-size energy levels. 

We have shown in Sect. ll.ll that a conformally invariant theory can be mapped on a cylindrical 
geometry, and the finite-size Hamiltonian can be expressed in terms of the Virasoro generators 
and the central charge as 

„ 2ir f - c 

ticFT = ~PT I + Lq 



R V u u 12, 

Perturbing this CFT with a relevant operator $ with zero spin and scaling dimension A = h + h, 
the Hamiltonian becomes (see 1)1.2.1)1 ) 



R 



H X = H C FT + XV, with V = j <S>(x,t)dx (4.1.4) 

o 

(we are considering for simplicity a single perturbing operator, but the discussion can be straight- 
forwardly generalized to multiple perturbations). The induced Renormalization Group flow is 
parameterized by the dimensionless variable Ai? 2_A , and the energy eigenvalues Ei can be 
expressed in terms of scaling functions fa of this variable: 

E t (R,X) =^ fi (\R*-*) . 

These functions interpolate between the massless CFT, corresponding to the ultraviolet (UV) 
limit XR 2 ~ A -> 0, and the massive infrared (IR) theory, corresponding to XR 2 A — ► oo. In the 
two limits the energies are dominated by 

Ei(R,X) ~ ^(hi + hi-^-) for Aii 2 - A ^0, (4.1.5) 



R V 12, 

2 - ->2-A 



Ei{R,X) ~ £ X^R + mi for XR —* oo , (4.1.6) 



where hi, hi are the dimensions of the conformal state related to the iih level, £q is the dimen- 
sionless bulk vacuum energy and mi is the (multi)particle mass term of the ith level. 

The TCSA technique takes advantage of the fact that the eigenstates of He ft can be chosen 
as a basis for the Hilbert space of the perturbed problem. As a consequence, the matrix elements 
of V in this basis 

<&l^l&> = £(&l*(0,0)|^-> 

can be computed in terms of the structure constants which enter the OPE in CFT (see (|1.1.6|> ). 
The obtained Hamiltonian can be finally truncated to the desired level, putting an upper bound 
on the conformal dimensions of the states in the basis. As a result of its diagonalization, one 
typically obtains finite—size spectra as shown in Fier. l4.2l (usually, one sets A = 1 and plots the 
spectrum as a function of R). 




Figure 4.2: TCSA results obtained in [HI] for the discrete spectrum Ei{R) in the Ising model with 
magnetic perturbation (left picture), and in the Ising model with both magnetic and thermal 
perturbation (right picture). 

From the UV and IR asymptotic behaviours of the obtained energy spectra it is there- 
fore possible to extract several conformal and scattering data, according to (|4. 1.5(1 and (|4. 1.6(1 . 
Furthermore, also the elastic two-particle S-matrix can be measured |2E], by considering two- 
particle states | A a A b ) with zero total momentum, relative momentum k and energy 

E = y/ml + k 2 + ^m 2 b + k 2 . 

In fact, the momentum k can be extracted from the above equation after the numerical de- 
termination of E, m a and m&, and the 5"-matrix can be then deduced by the quantization 
condition 

e tkR S ab (k) = 1. 

TCSA can be applied to both integrable and non-integrable theories, and this permits to 
observe an interesting phenomenon. In fact, in integrable theories it is quite common to observe 
crossing of different energy levels at finite values of R, while these intersections are immediately 
removed as one adds a non-integrable perturbation (see for instance Fig. l4.2j) . In non-relativistic 
quantum mechanics it is well known that two energy levels can cross only if they belong to differ- 
ent irreducible representations of the symmetry group of the Hamiltonian, otherwise, in absence 
of any symmetry, they have to repel. This has a clear counterpart in QFT: the presence of 
conserved charges ensures the stability of particles above threshold, whose lines cross an infinite 



number of one-particle lines. As soon as a perturbation is switched on, however, integrability 
is lost and inelastic processes become allowed, causing the decay of the particle. As a conse- 
quence, the particles lines above threshold become broken lines, whose difference in slope with 
the threshold is proportional to the decay width of the unstable states |59j . Therefore, the study 
of finite-size effects also gives evidences about the integrability of the theory in exam. 



4.1.2 Integrable QFT 

We have seen so far how important are finite-size effects for the understanding of QFT. It is 
therefore useful to have the possibility of controlling them also analytically. Apart from the case 
of CFT, where conformal invariance permits a complete treatment of the problem, interesting 
results can be obtained for integrable QFT, passing through a thermodynamic interpretation of 
the finite size of the system. 

The cylindrical geometry can be physically interpreted in two different ways, illustrated 
in Fig. 14. 31 One is the finite-volume picture described above, in which the space variable is 
compactified on a circle of circumference R while the time evolution takes place along an un- 
bounded direction. The other picture is the finite-temperature one, in which the space variable 
is unbounded but the (euclidean) time lives on a circle of circumference R. In the Matsubara 
formalism, this corresponds to introducing a finite temperature in the system, related to the 
compactification size by the inverse proportionality R oc l/T. 



Rot l/T 
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Figure 4.3: Finite volume and finite temperature interpretations of the cylindrical geometry. 



In the case of integrable QFT, a lot of information about the thermodynamics can be obtained 
from the knowledge of the S'-matrix. In fact, the partition function 

Z = Tr (e- RH ) 

can be determined by means of the so-called Thermodynamic Bethe Ansatz (TBA) technique, 
proposed in [00] • Here we limit ourselves to illustrate the basic ideas underlying this method, 
referring to the original literature for a consistent description. If we put the theory on a large 
box < x < L, assuming for simplicity that the spectrum contains only one type of particle 
with mass m and scattering matrix S, the quantization conditions of the momenta are given by 
the Bethe Ansatz equations 

jmLsmhOi - 6j) = 1 , 



with the rapidity variable 6 defined in (|1.2.3j) . The thermodynamic limit of the equivalent 
condition 

mLsinhfl, + ^ 6(&i - 6j) = litm , G Z , (4.1.7) 

with 6(6) = —i log S(6), is given by 

m cosh 6 + 2vr ip(6) * p x (6) = 2irp(6) , 

where p(6) and p±(6) represent, respectively, the density of levels and of occupied states per unit 
volume, * indicates the convolution operator and f(6) = -^6(6). The minimization of the free 
energy with respect to the density of states leads to the integral equation 

e(6) = mR cosh 6 -<p(6)* log(l + e- £ ( e) ) (4.1.8) 

for the so-called pseudo-energy e(6), defined as 

Pi _ 1 
p 1 + e £ ' 

The partition function can be then expressed in terms of the pseudo-energy, which is numerically 
determined from ()4.1.8|) . as 



Z(R, L) = exp 



mL J ^ cosh# log (l + e^ e ^) 



If we now come back to the finite-volume picture, performing the L — > oo limit we realize that 
the partition function has to be dominated by 

Z(R,L) = Tr (e~ LH *) ~ e~ LE °^ , 

where Eq(R) is the ground state eigenvalue of the finite-volume Hamiltonian Hr. Parameter- 
izing it conveniently as Eq(R) = ^/o(mi?), we therefore obtain the scaling function |6U1 161j 

f (mR) = J ^ cosh^ log (l + e-W) , (4.1.9) 

which in the UV limit has to reduce to f(mR) ^ ho + ho — c/12, where ho and c are the lowest 
conformal dimension and the central charge of the corresponding CFT (see (f 1. 1 . 13fl 1 . 

In this framework, the determination of excited energy levels is not possible or extremely 
involved. To this respect, however, other techniques exploiting integrability have been developed 
and applied for instance in |62| I63[ I64j . 



Besides energy levels, the understanding of QFT obviously requires also the knowledge of 
correlation functions. In the finite-temperature picture, the matrix elements are computed by 
performing a thermodynamic averaging over the usual infinite-volume particle states: 

(Cq ... O n ) R = ^Tv(e- RH O l ... O n ) , 

where H is the infinite-volume Hamiltonian and Z is the partition function. Therefore, in the 
case of integrable QFT, a method has been proposed in jHS] for expressing the finite-temperature 



correlators in terms of the infinite volume form factors and some thermodynamical quantities 
available from TBA. For instance, the one-point functions are given as 

(0(x,t)) R = ^3/^-^ (flWi)e- <e A <*„,.. .,*|O(O,O)|0 ll ...A> > 

71=0 ' 71 \i=l J 

where the "filling fractions" are defined as f{6) = (l + e _e ^) , with the pseudo-energy e(8) 
given by the solutions of (|4,1.8|) . Higher point correlation functions can be in principle ex- 
pressed in a conceptually similar way, although there is some controversy about their correct 
determination 

Further insight on this topic could come from a comparison with the finite-volume picture, 
as it was done for the ground state energy. In this case, the computation of matrix elements 
involves the notion of finite-size ground state \0) R ■ 

{d ... O n ) R = R {0\O 1 ... O n \0) R . 

Unfortunately, at present the precise structure of this ground state is not clear, and moreover, 
the finite-size matrix elements of operators are still poorly understood. It is precisely from 
this problem that we will start, in the next Section, our discussion of the semiclassical results 
obtained in finite volume. 

4.2 Semiclassical form factors in finite volume 

A particularly interesting problem in the study of QFT in finite volume regards the possibility 
of defining a "form factor representation" for the correlation functions, in analogy to the infinite 
volume expression 2.11(1 . There are reasons to expect, in fact, a fast convergent behaviour of 
these series also for finite volume correlators, as it happens in infinite volume (see, for instance 
|21j). If this would be indeed the case, accurate estimates of finite volume correlators and other 
related physical quantities, could be obtained by just using few exact terms of their spectral 
representations, having consequently a great simplification of the problem. This observation 
makes clear that it is worth pursuing the research on finite- volume form factors, looking in 
particular for an efficient scheme of approximation. 

In the non-perturbative study of form factors at a finite volume, there are so far only 
semiclassical results relative to a conformal theory |67| as well as exact calculations relative to 
the Ising model [HE] (for a Bethe Ansatz approach see, however, 123). Although these findings 
are very interesting, the techniques employed in the above papers are however strictly related 
either to the specific integrable structure of the considered models or to the free nature of the 
Major ana fermion field of the Ising model. 

In pQ, we have proposed to face this problem with a semiclassical approach based on Gold- 
stone and Jackiw's result. This method, contrary to the ones previously mentioned, does not 
require integrability and it is then of more general applicability, of course within its range of 
approximation. It is worth to recall an important feature that has come out from the study 
of the semiclassical form factors in infinite volume. As we have seen in Sect. l2.IHl their accu- 
racy seems to extend, somehow, beyond the regime in which they were supposed to be valid. 
Together with the known fast convergency properties of the spectral series, the above result 
suggests that the semiclassical method may provide a rather precise estimate of finite volume 
correlation functions, an outcome which may be useful for many applications. 



The application of the procedure described in Sect. l2.2l to the finite volume case is straightfor- 
ward, thanks to the possibility of choosing /(a) as a solution of eq. (|2,1.4|) with any constant A. 
As explicitly shown by the examples discussed below, this is equivalent to define a kink solution 
configuration on a finite volume, with the constant A directly related to the size of the system. 
We have now to consider the matrix elements of </>(0) between two eigenstates \p ni ) and \pn 2 ) 
of the finite volume hamiltonian Hr. These states can be naturally labelled with the so-called 
" quasi-momentum" variable p n , which corresponds to the eigenvalues of the translation operator 
on the cylinder (multiples of vr/L), and appears in the space dependent part of ea. (|2.2.3f) in 
the case of finite volume. The Bethe ansatz equations (|4. 1.7(1 . valid for large R, are exactly 
a relation between this variable and the free momentum p°° of the infinite volume asymptotic 
states, through a phase shift 5(p co ) which encodes the information about the interaction: 

«) _2nvr 
R ~ R ~ 

Defining 9 n as the " quasi-rapidity" of the kink states by 

p n = M(R) sinh 9 n ~ M(R)9 n , 

we can now write the form factor at a finite volume by replacing the Fourier integral transform 
with a Fourier series expansion: 

R/2 

f(9 n ) = (p n2 \m\P ni ) = M(R) J dae iM ^ a (p cl (a) , (4.2.1) 

-R/2 
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where 



(2m -l)vr (2n 2 -l)vr 2nvr 

M{R)9 n ~ Pni - Pn2 = y- 



R R R 

Since the energy eigenvalues of the finite volume hamiltonian cannot be related to the quasi- 
rapidity as in eq. (jl.2.3|) . in principle we are not allowed to express the crossed channel form 
factor F2(9) via the change of variable 9 — > m — 9. However, it is easy to show that in our 
regime of approximations the deviations of the kink energy from (|1.2.3|) are of higher order in 
the coupling and can be neglected at this stage. 

On the cylinder, the spectral function can be expressed as a series expansion on the form 
factors: 



P^E kiPk )=2^^ £ EkiE l_ EK KE k - E hl ... - E hn ) S( Pk - Pkl ... - Pkn )x 



x| (0| 0(O)| n) | 2 , (4.2.3) 

where p ki are the quasi-momenta of the intermediate states, and E k . are the finite volume energy 
eigenvalues, that will be discussed in the following Sections. For the moment their knowledge is 
not necessary, because at leading semiclassical order the spectral function pW (E k ,p k ) is given 
by the trivial vacuum term plus the kink-antikink contribution, and the kink energies can be 



consistently approximated with their classical masses M, which can be exactly computed as 
functions of the volume. We then have 



pM(E k ,p k ) = 2v6(E k )6(p k )\ ( 0| 0(O)| ) | 2 + \ KM M2 > £ |*b (29 kl + m - g 



(4.2.4) 

As in the infinite volume case, the consistency of the semi-classical approximation selects as 
the relevant values of the above series those with 9 k ~ and therefore it can can be roughly 
estimated by simply evaluating I-F2I 2 at &ki = 0. 

We will now explicitly construct the form factors on the kink states in the sine-Gordon 
model (j2.3.1|) and in the broken 4 theory (|2.1.5|) . both defined on a cylindrical geometry with 
antiperiodic boundary conditions, respectively given by 

{Six + R) = -£■ — Mx) for sine-Gordon , 

V ; P K ' (4.2.5) 
<p(x + R) = —<j){x) for broken 4 . 



Sine Gordon model 

In order to identify a kink on the twisted cylinder, we have to look for a static finite energy 
solution of the SG model satisfying the anti-periodic boundary condition ()4.2.5|) . For the first 
order equation 

l{^t) 2 = % {1 - COSP4>d+A) (42 - 6) 
a solution with this property can be found for — 2 < A < 0, and it is expressed as 

2 

^d( x ) = t; arccos [k sn (m{x — xq), A: 2 )] , (4.2.7) 

y 

where sn(u, k 2 ) is the Jacobi elliptic function with modulus k 2 = ^J 2 -, and the presence of a free 
parameter x$, which represents the kink's center of mass position, is due to the translational 
invariance of the theory around the cylinder axis. The plot of (|4.2,7[) as a function of the real 
variable x is drawn in Fig l4.4l For a given value of A, this solution oscillates with a period 
4K(fc 2 ) between 4>q and ^ — 4>o, where <po is defined by the condition V(4>o) = — A, and 
K(/c 2 ) is the complete elliptic integral of the first kind. For the definitions and properties of 
elliptic integrals and Jacobi elliptic functions, see Appendix I4IB1 

The solution (|4.2.7|) has been proposed in [JO] as a model of a crystal of solitons and antisoli- 
tons in the sine-Gordon theory in infinite volume. In our finite volume case, the solution (|4.2.7|) 
has to be interpreted, instead, as a single (anti)soliton defined on a cylinder of circumference 

R = — 2K (A; 2 ) . (4.2.8) 
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Within this interpretation, the periodic oscillations of the solution represent the soliton circling 
around the cylinder. Eq. (|4.2,8|) is the explicit relation between the size of the system and the 
integration constant A; one can consistently recover the infinite volume limit for A — > 0: in 
this limit R goes to infinity and the function (|4.2,7|) goes to the standard (anti)soliton solution 

HSU- 




We can now write the finite volume form factor (|4.2.1j) in terms of the antikink background 

6221): 
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(4.2.9) 
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where k' = y/l - k 2 and K'(k 2 ) = K(k' 2 ). The kink mass M is represented at this order by the 
classical energy, given by 
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S d {R)= j dx 

-R/2 



d ~) + F (l-co S /3^) 



m 

~2 



B(k 2 ) - -(1 - k 2 )K(k 2 ) 



(4.2.10) 

where E(fc 2 ) is the complete elliptic integral of the second kind. In order to obtain the result 
()4.2.9|) one has to use the relation 

arccos [k sn(ma)] = - log [k sn(ma) + i dn(ma)] , 



and, after an integration by parts, finally compare the inverse Fourier transform ()4.2.2j) with the 
expansion [201 
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The form factor ()4.2.9j) has the correct IR limit 1 , and leads to the following expressions for 



x The function - — can be shown to tend to — iTv5(x) in the distributional sense for R — > oo, and in the 

same way one can show that cot, ( xR / 2 ) tends to zero. 



i*2(0) and for the spectral function: 
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(4.2.11) 
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(4.2.12) 

Note that the finite volume dependence of both the form factor 1)4.2. lljl and the spectral function 
1)4.2.12)1 is not restricted to the second term only. The M(R) K'(/c 2 ) factor in the first term carries 
the main independence, although it is not manifest but implicitly defined by eq. 1)4.2.8)) . 



Broken eft 4 theory 

The differential equation 
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has a solution given by 



•> d (x) = (±)<j2-<f% sn ( -j= (x - s ), k 2 



(4.2.13) 



with k 2 = -Jj — 1, V(4>o) = —A and 1 < (f>o < y/2, where we have rescaled the variables as 



m 



x = mx . 



This function oscillates with period 



2 R = 1 )H K (k 2 ) 

m (p 



(4.2.14) 



between the two values — y2 — </>q and y 2 — 0q, and it satisfies the anti-periodic boundary 
condition 1)4.2.5)1 . Moreover, it goes to the standard (anti)kink solution 1)2.1.6)1 for ^>o — ► 1, i-e. 
when R — > cxd. 

From the analytic knowledge of the background 1)4.2.13)) . we can immediately extract an 
important scattering data of the non-integrable cj) 4 theory. In fact, the leading term in the kink 
mass is given by the classical energy, expressed for generic R as 



£d(R) 
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(4.2.15) 



It is easy to see that for R — > oo this quantity indeed reproduces the infinite-volume energy 
1)2.1.7)1 . From its asymptotic expansion for large R, we can obtain the leading order of the kink 
mass correction in finite volume, and compare it with Liischer's result 1)4. 1.2114.1. 3)1 . Taking into 



account the k — > 1 (kf — ► 0) expansions of E and K (see Appendix l4IB|) and noting from (|4.2.14() 
that 

V2mR _ 



1 

256 1 



!- + ■•• , 

we derive the following asymptotic expansion of E c \ for large R: 

S d (R) =M 00 -8V2 Hf- e -^mR + Q ^-2V2mR 



(4.2.16) 



The counterpart of this behaviour in Liischer's theory is given by the process (1) in Fig. l4.H 
where particle a is the kink, particle b is the elementary meson, and particle c is another kink. In 
fact, in the broken </> 4 theory the elementary boson has semiclassical leading mass m& = \/2m, 
which is much lower than the kink one (|2.1.7|) . therefore the dynamical pole relative to this 
process is located at u kb ~ ^. From the comparison with (|4. 1.2(1 we finally extract the leading 
semiclassical expression for the residue of this 3-particle process: 



Rkkb — 8 



m 



It is easy to see that this result correctly reproduces, as required by (|1.2.8|) . the square of the 
value (|3.2.4() previously obtained for the 3-particle coupling T k j ib by looking at the residue of 
the kink-antikink form factor in infinite volume 2 . 

Moving then to the finite-volume form factors, these can be computed by comparing the 
inverse Fourier transform (|4.2.2|) with the expansion |90| 
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obtaining the expression 
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Therefore, the 1/A leading contribution to the spectral function is given by 
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(4.2.17) 



(4.2.18) 



Again, as in the Sine-Gordon case, the finite volume dependence of these quantities comes from 
the factor M(R)K'(k 2 ), where M(R) is the kink mass given by (|4.2.15j) . 



4.3 Sine— Gordon model on the cylinder 

As we have discussed, the knowledge of form factors permits to estimate the the spectral density 
representation of correlation functions at a finite volume in both integrable and non-integrable 

2 Crossing symmetry implies the equality R k ^. h = Rkkb- 



theories. However, correlation functions need another set of data for their complete determina- 
tion, precisely the energies of the intermediate states at a finite volume. 

The paper [3] has been mainly devoted to fill this gap, that is, to face the problem of a 
semiclassical computation of the energies Ei{R) of vacua and excited states as functions of 
the circumference R of a cylindrical geometry. The analytic form of the semiclassical scaling 
functions for two-dimensional QFT admitting static kink solutions can be achieved by suitably 
adapting the DHN method to the finite geometry. The example discussed in is the sine- 
Gordon model, which is particularly appealing for its simplified semiclassical results whereby 
the significant physical effects are not masked by other additional complications. Moreover, 
due to the integrable nature of this theory, its finite size effects have been previously studied 
by means of Thermodynamical Bethe Ansatz [621 I64j . and it would be interesting to perform 
a quantitative comparison between these results and the semiclassical ones, in order to directly 
control their range of validity. However, as already pointed out, semiclassical methods apply not 
only to integrable theories and this opens the way to describe analytically the finite size effects 
also in non-integrable models. 

The Sine-Gordon model ((2.3.1(1 on a cylindrical geometry admits quasi-periodic boundary 
conditions (b.c.) 

<f>(x + R,t) = 4>{x,t) + — 

(the arbitrary winding number n £ Z originates from the invariance of the potential 1(2.3. lj) 
under <f) — > <fi + ^p). In particular, since we are interested in the one-kink sector, which is 
defined by n = 1, we will impose the b.c. 

2ir 

<p(x + R,t) = if>(x,t) + — . (4.3.1) 

P 

The first step for applying the semiclassical method to this problem is to find the finite size analog 
of the kink solution, satisfying now the b.c.'s ((4.3.1J) . However, the success in constructing 
the scaling functions depends on whether one is able to solve the corresponding Schrodinger 
equation (|2. 1.9(1 and to derive an analytical expression for its frequencies uj}.. It turns out that 
the semiclassical finite size effects in SG model are intrinsically related to the simplest (N = 1) 
Lame equation, which admits a complete analytical study. 

The subject is organized as follows: after the introductory Section l4.3.1( where we discuss 
the simplest scaling function in a finite volume in order to clarify the nature of divergencies 
encountered in such computations, in Section l4.3.2l we present the complete semiclassical analysis 
of the energy levels in the one-kink sector. Finally, Section l4.3.3l is devoted to the finite-size 
form factors for the SG model on a periodic cylinder. 



4.3.1 Ground state energy regularization 

As shown by eq. (|2. 1.11(1 . quantum corrections to energy levels are given by the series on the 
frequencies u n . However, this series is generally divergent (this is the usual UV divergence in 
field theory) and a criterion is needed to regularize it. It is quite instructive to consider the 
simplest example where such divergence occur, i.e. in the calculation of the ground state energy 
£q 3C (R) of the vacuum sector of the SG theory on a cylindrical geometry of circumference R. 
This can be constructed by implementing the DHN procedure for one of the constant solutions, 
for instance <^f' c = 0, imposing periodic boundary conditions for the corresponding fluctuations 



7/ vac (x). Obviously, what comes out is nothing else but the Casimir energy of a free bosonic field 
(j>(x, t) with mass m. In this case the frequency eigenvalues are fixed to be 



with p n = 2-rrn/R and n = 0, ±1, ±2, .... 

The ground state energy has to be regularized by subtracting its infinite-volume continuous 
limit: this ensures in fact the proper normalization of this quantity, expressed by 

hm er(R) = o . 

R— >oo 

The ground state energy at a finite volume is therefore defined by 
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Isolating the zero mode, it can be conveniently rewritten as 
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where r = mR. Since the divergence of the series is due to the large n behaviour of the first two 
terms in the expansion 
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we begin our calculation by subtracting and adding these divergent terms to it: 
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The first series in the right hand side of the above expression is now convergent, whereas the last 
two terms should be coupled to the analogous ones coming from the integral, whose divergencies 
have to be handled in strict correspondence with those coming from the series. Hence, by 
subtracting and adding the leading divergence to the integral 
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n 2 + J -nj + J dun, (4.3.4) 

we can combine the last term in this expression with the one in (|4.3.3|) and implement the well 
known regularization 
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However, the first term in 1)4.3.41) still contains a subleading logarithmic divergence, as it can be 
seen by explicitly computing the integral by using a cut-off A, in the limit A — ► oo 



A , 

J dn { 



n z + 



n 



l/r\ 2 , , l/r\ 2 l/r\ 2 , r 
- I — I In 2A + - ( — ) - r ( — ) In — 



2 V2vr/ 



4 V2tt/ 



2 V2vr/ 2vr 



(4.3.6) 



This divergence can be cured by subtracting and adding the term ^ (^) 2 lnA. By combining 
this last term with its analogous in the series we have 
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where 7^ is the Euler-Mascheroni constant, while the remaining part of 1)4.3.6)) with the above 
subtraction is now finite. 

Collecting the above results, the finite expression of the ground state energy on a cylinder 
is then given by 
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(4.3.8) 

It is now easy to see that (|4.3.8|) nicely coincide with the analogous expression obtained in 
the finite-temperature picture, given by the TBA result (|4.1.9f) for the free case, in which the 
pseudo-energy is simply given by e{9) = mRcosh.0: 
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In fact, this integral formula can be expressed in terms of Bessel functions, which admit a series 
representation that directly leads to 1)4.3.8(1 (see Ref. |E2)- Moreover, one can also check that 
the above regularization scheme ensures the agreement between the finite-volume and finite- 
temperature calculations of the one-point functions {4> 2k ). The interested reader can find the 
simplest example of these calculations in Appendix I4IAI 

Finally, it is worth noting that the result 1)4.3.8(1 can also be obtained by using a simpler 
prescription which automatically includes the subtraction of the various divergencies, fastening 
the calculation. This consists in ignoring the divergent part of the integral, keeping only its 
finite part, and in regularizing the divergent series as 
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(4.3.9) 
(4.3.10) 



Formula 1)4.3.9)1 is the standard regularization of the Riemann zeta function £(— 1), where £(s) = 
Sn=i n 3 " ' an d usually corresponds to normal ordering with respect to the infinite volume vacuum 
(see, for instance, chapter 4). On the contrary, the regularization of the second series is 
a-priori ambiguous due to its logarithmic divergence, and its finite value 1)4.3.10)1 was chosen 
according to the above discussion. 



4.3.2 Scaling functions 



We will now develop a complete semiclassical scheme to analyse the energy of the quantum 
states in the one-kink sector of SG model on the cylinder. This can be achieved by applying 
the DHN method to an appropriate kink background. 

Properties of the periodic kink solution 

In order to identify a kink on the cylinder, we have to look for a static finite energy solution 
of the SG model satisfying the quasi-periodic boundary condition 1)4.3.1 J) . For the first order 
equation 1)4.2.6)) a solution which has this property can be found for A > 0. It can be expressed 
as 

7T 2 / mix — Xn) 9 \ , 9 2 . 

Mx) = ? + - 3 ,m (-^r^.^) . ^ = (4.3.11) 

provided the circumference R of the cylinder is identified with 

R = —kK(k 2 ) , (4.3.12) 

m 

where K.(k 2 ) denotes the complete elliptic integral of the first kind 3 . The parameter xo in 1)4.3.11)) 
represents the kink's center of mass position, and its arbitrariness is due to the translational 
invariance of the theory around the cylinder axis. The behaviour of 1)4.3.11)1 as a function of the 
real variable x is shown in Figure ll~o1 
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Figure 4.5: Solution of eq. 1)4.2.6)1 with A > and xq = 0. 

The function (|4.3.11|) has been first proposed in [70] and interpreted as a crystal of solitons 
in the sine-Gordon theory in infinite volume. In our finite volume case, instead, ()4.3.11)) has to 
be seen as a single soliton defined on a cylinder of circumference R (given by eq. 1)4.3.12)1 ). while 
its quasi-periodic oscillations represent winding around the cylinder. As shown in eq. ()4.3.12j) . 
there is an explicit relation between the size of the system and the integration constant A. It 

3 The definition and basic properties of K(fc 2 ) and the Jacobi elliptic amplitude am(u, k 2 ) can be found in 
Appendix I4IBI 



is easy to see that the infinite volume solution (|2.3.3|) is consistently recovered from (|4.3. llf) in 
the limit A — > 0, i.e. when R goes to infinity. 

The classical energy of the kink on the cylinder is given by 
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(4.3.13) 

where E(fc 2 ) is the complete elliptic integral of the second kind. In the R — > 00 limit (which 
corresponds to k' — > 0, with (A;') 2 = 1 — k 2 ), £ c i(R) approaches exponentially the correct value 
Moo = • This can be seen expanding E and K for small k' (see Appendix l4IB|) , and expressing 
the result in terms of mR, which can be itself expanded in k' in virtue of the relation (|4.3.12|) : 

-mR = 2_n t f)* + ... . 

16 v ; 

Hence the large R expansion of the classical energy is 
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£ cl (R) = M^ + ^m e~ mR + O {e~ 2mR ) . (4.3.14) 

We will comment more on the interpretation of this result in the following. 

Similarly, one can derive the behaviour of £ c i{R) for small r = mR, which corresponds to 
the limit A — > 00 (or k 2 — > 0) : 

£d{R) = -^^ + m^-m( — ) + (4.3.15) 
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This formula will be relevant in the later discussion of the UV properties of the scaling functions. 

Before moving to the quantization of the kink-background (|4.3.11|) , it is worth recalling that 
another simple kind of elliptic function, given in (|4.2.7|) . was also proposed in [70] and interpreted 
as a crystal of solitons and antisolitons in the infinite volume SG. This background corresponds 
as well to a kink on the cylinder geometry but satisfying the antiperiodic boundary conditions 

2-7T 

(j){x + R,t) = —-4>{x,t) . 

The associated form factors were obtained in [I] and have been described in Sect. l4.2l Although 
the quantization of this second kink solution is technically similar to the one of (|4.3.11f) presented 
here, it displays however some different interpretative features that justify its discussion in a 
separate future publication |72j . 



Semiclassical quantization in finite volume 

The application of the DHN method to the periodic kink (|4,3.11[) requires the solution of 
eq. (|2. 1.9(1 for the quantum fluctuations r/ w , which in this case takes the form 

\d% + k2 (° 2 + l ) -2k 2 sn 2 (x,k 2 )}^(x) =0, (4.3.16) 

where sn(x, k 2 ) is the Jacobi elliptic function defined in Appendix I4llj| and we have introduced 
the rescaled variables 

mx _ uj 

x = — , uj = — . 
k m 



Due to the periodic properties of (j) c i{x) expressed by eq. (|4.3.11j) . the boundary condition (|4.3.1I) 
translates in the requirement for rj^(x) 

rjo + = Vo(x) ■ (4.3.17) 

Eq. (|4. 3.16(1 can be cast in the so-called Lame form, which admits the two linearly independent 
solutions 

V±a(x) 
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a(x + iK') 

where the auxiliary parameter a is denned as a root of the equation 

2 - k 2 

V{a) = . 

The Weierstrass functions V(u), C( u ) an d cr(u) are defined in Appendix I4IC[ where the Lame 
equation and its relation with 1)4.3. 16JI are discussed in detail. 

As it is usually the case for a Schrodinger-like equation with periodic potential, the spectrum 
of eq. 1)4.3.16)) has a band structure, determined by the properties of the Floquet exponent 

F(a) = 2i[K£(a)-a((K)] , (4.3.18) 

which is defined as the phase acquired by r/± a in circling once the cylinder 

V±a (x + 2K) = e ±lF ^rj ±a (x) . 

We have two allowed bands for real F(a), i.e. 

-9 1 -■? 1 

< u K iu2 ~~ 1 and u > U2 > 

and two forbidden bands for F(a) complex, i.e. 

-9 1 1 

u < and - 1 < uj < . 

The band < < is described by a = K + iy, where y varies between and K' and, 
correspondingly, F(a) goes from to ir. The other allowed band lo 2 > corresponds instead 
to a = iy and, by varying y, F(a) goes from tt to infinity, as it is shown in Fig. 14.61 

By imposing the periodic boundary conditions (|4.3.17|) on the fluctuation rj(x), one selects the 
values of uj 2 for which the Floquet exponent is an even multiple of 7T, thus making the spectrum 
of eq. (|4.3.16|) discrete. These eigenvalues are (Dq = 0, which is the zero mode associated with 
translational invariance and has multiplicity one, and the infinite series of points 
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with multiplicity two, placed in the highest band uj 2 > p-, with y n determined by the equation 

F(iy n ) = 2K % C(ty n ) + 2y n C(K) = 2nvr , n = 1, 2, . . . (4.3.20) 

In the IR limit (A — > 0) the spectrum goes to the one related to the standard background (|2,3.3|) : 
the allowed band < uj 2 < p- — 1, in fact, shrinks to the eigenvalue <Dq = 0, while the other 
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Figure 4.6: Spectrum of eq. (j4.3.1(ij) 



allowed band uj 2 > merges in the continuous part of the spectrum Q 2 = 1 + q 2 . Therefore, 
by recalling the interpretation of the continuous spectrum in infinite volume as representing the 
scattering states of the kink with breathers, we can now give to the requirement (|4.3.20|) the 
physical meaning of a quantization condition for the momenta of the states constituted by a 
kink and a breather on the finite volume, analogous to the Bethe ansatz equation (|4.1.7j) but 
valid for any value of the size of the system. In particular, using the small-Zc' expansions of 
the elliptic integrals and Weierstrass functions presented in Appendices 141131 and I4IC Jl it easy to 
check that for large R the condition (j4.3,2()j) indeed reduces to the form 

mR q + 5(q) = 2nn , 

where mq represents the momentum of the lightest breather, and 5(q) is the semiclassical phase 

shift (Un- 
it is useful to note that, although the R dependence of the frequencies ()4.3.19j) is quite 

implicit, since it passes through the inversion of eq. (|4.3.12|) . nevertheless these are analytic 

functions of R and it is extremely simple to plot them. The corresponding curves, shown in 

Figure l4~7| provide an important piece of information, since they are nothing else but the energies 

of the excited states with respect to their ground state £o(R)- 

To complete the analysis, it remains then to determine the finite volume ground state energy 

£o(R) of the kink sector. In analogy with the infinite volume case (see eq. (|2.3.5|0 . this is defined 

by 

R/2 

£ (R) = £ d {R) + ]T u n {R) -JjL J dx[l- cos/fyc] - • (4-3.21) 

n=1 -R/2 

Before commenting in detail each of these terms, let's focus first on the main problem in de- 
riving a closed expression for £q(R), which consists in the evaluation of the infinite sum on the 
frequencies uj n (R) or, better, in isolating its finite part. We need therefore a method for solving 
the transcendental equation (j4.3.20f) for y n (k 2 ) in order to make the expression (|4.3.19|) for the 
frequencies oj n (k 2 ) explicit. As we have already seen for the classical energy, two kinds of expan- 
sion are possible, one in the elliptic modulus k and the other in the complementary modulus k' , 




which are efficient approximation schemes in the small and large r regimes, respectively. Here 
for simplicity we only present the small r expansion. By taking into account the series expansion 
in k for K, ((u) and V(u) (see Appendices I4IBI and I4I(!|) . we are led to look for a solution of 
eq. (j4,3.20j) in the form 
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Here we give the result for the first few coefficients y„ , s = 0, 1, 2: 
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which are those relevant in the later analysis of the UV properties of the scaling function. As a 
consequence, we obtain the following simple expression for the frequencies: 
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Comparing it order by order with the small-fc expansion of eq. (|4.3.12|) 
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we finally obtain the explicit i?-dependence 
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(4.3.22) 



It is worth noting that this series expansion in r, which can be easily extended up to desired 
accuracy, efficiently approximates the exact energy levels also for rather large values of the 
scaling variable. Fig. 14.81 shows a numerical comparison between the first energy level given by 




Figure 4.8: Comparison between the exact energy level oj\/m given by (|4.3.19j) (continuous line) 
and the approximate expression (j4,3.22[) (dotted line). 

()4.3.19|) and its approximate expression (|4,3.22j) . and for the higher levels it is possible to see 
that the agreement is even better. 

With the above analysis, the sum over frequencies in the ground state energy (j4.3.21j) takes 
the form 

71 = 1 71 = 1 71=1 



As we will see below, the subtraction of counterterm and vacuum energy in (|4.3.21l) leads to the 
cancellation of all the divergencies, producing a finite expression for the ground state energy in 
the kink sector. 

Moving now to the analysis of the remaining terms in (j4.3.21|) . a similar series expansion 
can be easily performed on each of them. The classical energy £ C [(R), given in eq. (j4.3.13|) . has 
already be treated in this way in eq. (j4.3.15|) . The finite volume counterterm (C.T.), where the 
one-loop mass renormalisation is given by 

c 2 m 2 (3 2 2vr ^ 1 
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and 4> c i is given by ()4.3.11j) . takes the explicit form 

K(k 2 )-B(k 2 ) 



C.T. = m 



kK (k 2 ) 



E 



a/(2?t,7t) 2 + r 2 



The first terms of its expansion in R are then 
C.T. 
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Finally, the vacuum energy <5g ac (i?) is the one precisely computed in Sec. 14.3.11 Since its role is 
to cancel certain divergencies present in the other terms of £q(R), in complete analogy with the 



infinite volume case (see eq. (|2.3.5|l ). we will now consider its "naive" formulation, given by 

£ ^ C (R) J_ l/2mr\ 2 ~_1 2vr ^ I(M 3 y^ — 
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(4.3.25) 

Hence, in the final expression for the ground state energy all the divergent series present in the 
sum over frequencies, in the counterterm and in the vacuum energy cancel out, and one obtains 
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where we have used |90j 



(4.3.26) 
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in order to evaluate explicitly the coefficient of the r 3 term. 

Repeating the above calculations, one can also easily write the finite expressions of the 
excited energy levels 1)2.1.10)1 . whose series expansion in r is given by 
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where is a set of integers defining a particular excited state of the kink. 
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UV- IR correspondence 



The semiclassical quantization of the periodic kink 1)4.3.111) provides us with analytic expressions, 
albeit implicit, for the scaling functions in the kink sector for arbitrary values of the scale 
r = mR. These quantities control analytically the interpolation between the Hilbert spaces 
of the ultraviolet (UV) and infrared (IR) limiting theories. It is worth noting that, although 
we obtain them in the framework of a particle-like description proper of the IR limit, the 
kink background l)4.3.11|) is intrinsically formulated on a finite size, leading to the possibility of 
extracting UV data. Hence, it is important to check whether our scaling functions reproduce 
both the expected results for the IR (r — ► oo) and UV (r — ► 0) limits. 

Concerning the IR behaviour, we have already seen that in the R — > oo limit all the quantities 
in exam, i.e. the classical solution, its classical energy and the stability frequencies, correctly 
reach their asymptotic values. In addition, we can now analyse the asymptotic approach of the 
kink mass to the infinite- volume limit, in order to compare it with the known scattering data 
of the SG theory, according to Liischer's theory described in Sect. l4.ll Restricting for simplicity 
our analysis to the leading term in (5 in the kink mass (which, in our approach, is simply given 
by the classical energy), we have then to compare the expansion presented in 1)4.3.14)1 with the 
term that dominates Liischer's result for small (3. This is given by formula 1)4.1.2)1 evaluated for 
the process in which a kink splits into the lightest breather and another kink, i.e. 
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Recalling the kink-breather S-matrix ^f)] 
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and selecting its s-channel pole 6** = iiij£ 6 = i (f + we find 4 
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Substituting in 1)4.3.28)) . for small 1 we have 
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AM(R) = m ^ e" mR + . . . , 

which therefore reproduces eq. ()4.3.14|) . It is a remarkable fact that the classical energy alone, 
being the leading term in the mass for f5 2 — ► 0, contains the IR scattering information which 
controls the large-distance behaviour of £q(R). 

As we have seen in Sect. l4.l| the UV behaviour for r — > of the ground state energy Eq(R) of 
a given off-critical theory is related instead to the Conformal Field Theory (CFT) data (A, A, c) 
of the corresponding critical theory and to the bulk energy term as 

27T / — C \ 

E (R) ~ —(h + h- — ) + BR + --- (4.3.29) 
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where c is the central charge, h + h'vs the lowest anomalous dimension in a given sector of the 
theory and B the bulk coefficient. For the Sine-Gordon model the bulk energy term is given by 

HUES] 

o 

an *y 

B = 16 tan — , (4.3.30) 
7 Z 16 

while its UV limit is described by the CFT given by the gaussian action (jl,1.14j) with the 
normalization constant fixed to the value 5 g = 1 

Ag = \ J d2xd ^ d ^i 

and the free bosonic field compactified on a circle of radius TZ = (see ea. ()1.1.15j) ). We have 
described the main properties of this CFT in Sect. ll.ll In particular, the central charge takes 
the value c = 1, and the theory is divided into sectors labelled by two integers s and n. In each 
sector, the state with lowest anomalous dimension (|1.1.17JI is created by the vertex primary 
operator V s>n , introduced in 1)1.1.16(1 . It is worth noticing that the SG action can be seen as 
the perturbed gaussian CFT of the form (I4T41 with A = ^ and $ = Vi, + V-i,o. These 

perturbing operators carry anomalous dimension A±i t o = /i±i,o + ^±1,0 = which is zero in 
the semiclassical limit, determining the form of the scaling variable as r = mR. 

The vacuum sector is described by s = n = 0, with h vac + h vac = 0, while the kink sector, 
defined by the boundary condition 1)4. 3.1(1 . naturally corresponds to s = 0,n = 1, in which the 
lowest anomalous dimension is 

A ,i = V + h ,i = ■ (4-3.31) 



The negative sign of the residue is due to the odd parity of the lightest breather, which implies F fc = — r fcb 

in turn 

5 Note that the usual normalization adopted in the CFT literature is instead g — j-. 



The conformal vertex operator Vo 1 has been put in exact correspondence with the soliton- 
creating operator of SG in Ref. |73] . 

The question to be addressed now is whether the small r expansion of £Q ac {R) and £$(R) 
given by eqs. (|4.3.8|) and (|4.3.26|) reproduces, in semiclassical approximation, the above data 
controlling the UV limit of SG model. For the vacuum sector, comparing (|4.3.8() with (|4.3,29j) . 
we correctly obtain c = 1 and h vac = h vac = 0. We do not expect, however, to obtain the bulk 
term B relative to SG model by looking at (|4.3,8[) . simply because the semiclassical expression 
of the ground state energy in the vacuum sector applies equally well to any theory which has 
a quadratic expansion near the vacuum state. Namely, apart from the value of the mass m, 
eq. (|4,3.8j) is a universal expression that does not refer then to SG model. The kink scaling 
function (|4.3.26|) has instead a richer structure. The obtained scaling dimension 

7 T ^ 
k + h = j2 

is the expected one, given by ()4.3.31|) . for the soliton-creating operator in Sine-Gordon, while the 
central charge contribution c = 1 is absent, simply because in (|4.3.26|) we have subtracted the 
vacuum ground state energy from the kink one 6 . Moreover, the bulk coefficient B = present 
in (|4.3.26|) correctly reproduces the semiclassical limit of the exact one, given in eq. (|4.3.30|) . In 
principle, this bulk term should be present in all the energy levels, included the ground state 
energy in the vacuum sector, but its non-perturbative nature makes impossible to see it in the 
semiclassical expansion around the vacuum solution, which is in fact purely perturbative. Hence 
it is not surprising that to extract the bulk energy term we have to look at the kink ground state 
energy, in virtue of the non-perturbative nature of the corresponding classical solution. Finally, 
the expression (|4,3.27|) for the excited energy levels explicitly show their correspondence with 
the conformal descendants of the kink ground state. In fact, their anomalous dimension is given 
by 

k {kn} + hkn} = ^2 + Yl ^ U ■ (4.3.32) 

The successful check with known UV and IR asymptotic behaviours confirms the ability 
of the semiclassical results to describe analytically the scaling functions of SG model in the 
one-kink sector. It would be interesting to further test them at arbitrary values of r through a 
numerical comparison with the results of |621 164j in an appropriate range of parameters. This 
was not pursued here because the results presently available in the literature were obtained 
for values of (3 which are beyond the semiclassical regime and moreover the energy levels were 
plotted as functions of a different scaling variable, i.e. the one defined in terms of the kink mass. 
We hope however to come back to this problem in the future. 

4.3.3 Form factors and correlation functions 

The semiclassical scaling functions, derived in Sect. l4.37fl provide an important information 
about the finite size effects in SG model. As in the infinite volume case, however, the complete 
description of the finite volume QFT requires to find, in addition to the energy eigenvalues 
(|4.3.27|) . the kink form factors and the correlation functions of local operators. This section 

6 The value c = 1, coming out from the regularization of the leading term of the series on the frequencies 
14.3.22H . is in fact exactly cancelled by the same term in the vacuum energy. 



is devoted to the analysis of this problem, i.e. to the determination of the finite volume form 
factors and the corresponding spectral functions. 

The form factors can be determined with the same procedure described in Sect. 14.21 and 
applied to the SG model and the broken </> 4 field theory, both defined on a cylindrical geometry 
with antiperiodic boundary conditions. In what follows we will apply it instead to the case of SG 
model with periodic boundary conditions. The corresponding finite volume form factor 1)4.2. lj) 
can be written in terms of the soliton background 1)4.3.11)) : 
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In order to obtain this result one has to compare the inverse Fourier transform Q4.2.2|) with the 
expansion [90] 
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The form factor (|4.3.33j) has the correct IR limit , and leads to the following expressions for 
F2(0) and for the spectral function 8 : 
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(4.3.34) 



(4.3.35) 



where 6 = in — 9. Note that the finite volume dependence of both the form factor ()4.3.34|) 
and the spectral function (|4,3.35|) is not restricted to the second term only. The /cK'(/c 2 ) M(R) 
factor in the first term carries the main independence, although it is not manifest but implicitly 
defined by eq. (I4.3.12J) . 

Another quantity of interest is the two-point function ( | e(x)e(0) \ ) of the energy density 
operator. One can calculate it by evaluating its spectral function 
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7 The functions cos ( xR / 2 ^ anc j sm [f a / 2 ) can be shown to tend to zero in the distributional sense for 7? — > oo. 
8 Here we are considering the matrix elements on the antisymmetric combinations of kink and antikink. 



in terms of the form factors of e(x), similarly to what we have done for the field <j) above. In 
order to find the semiclassical form factor 



fe(9n) = (Pn 2 |e(0) \Pni) , 

we need to compute the Fourier transform of 

, 2m 2 ,, , 9n 4m 2 9 fmx\ 

<^ = W* {1 + k) -JP m (x) 

This can be easily obtained from the following expansion 
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and we finally have 



Ue„)=M* . (4.3.37) 



The corresponding semiclassical spectral function is thus given by 



P (£) (^n,Pn) = ^ 5 (§ - 2) j\ N . (4.3.38) 
Z 3 V M / sinh 2 ( ^ 



It is worth mentioning that it is also possible to obtain the two-point functions of certain 
vertex operators V^(x, t) = e ±l P b <t>(xt) (f or 5 = i ; 1, | ; 2, ...), since the required Fourier expansion 
formulas of the type (|4.3.36|) are known in these cases [2^ • 



4.4 Strip geometry 

In our study of finite-size effects, we have focused so far only on the cylindrical geometry, since 
the relevant features we were interested to underline can be already observed in this simple case. 
However, an even richer phenomenology is produced by including non-trivial boundaries, and 
the semiclassical techniques are also suited to describe such geometries. In this Section, after a 
brief overview of boundary effects in CFT, we will describe how the semiclassical method has 
been applied in [1] to study a quantum field theory defined on a strip of width R, with certain 
boundary conditions at its edges. In particular, the example discussed is the sine-Gordon model 
subjected to Dirichlet boundary conditions at both edges of the strip. 



4.4.1 Boundary effects in CFT 

With the purpose of giving an intuitive idea of boundary effects in QFT, we focus here on the 
case of conformally invariant theories, due to the powerful analytical techniques available in 
this situation. We limit ourselves to a brief description of the main results, referring to the 
original literature for a complete discussion lllj. Furthermore, it should be mentioned that 
several exact results have been obtained also for integrable QFT with boundary, in virtue of the 
factorization property of the scattering in the bulk and off the boundary |75j . 

In a critical system with a boundary, conformal transformations must map the boundary onto 
itself and preserve the boundary conditions. As a consequence, holomorphic and antiholomorphic 



fields no longer decouple, and only half of the conformal generators remain. This can be easily 
seen in the prototype geometry for a two-dimensional system with boundary, i.e. the upper half 
plane. In fact, infinitesimal local conformal transformations of the form z — > z + e(z) map the 
real axis onto itself if and only if e(z) = e{z), i.e. e is real on the real axis, and this constraint 
eliminates half of the conformal generators. 

A powerful tool for describing the half-plane geometry is the so-called method of images. It 
consists in regarding the dependence of the correlators on antiholomorphic coordinates z% on the 
upper half-plane as a dependence on holomorphic coordinates z* = Z{ on the lower (unphysical) 
half-plane. Having introduced in this way a mirror image of the system, the definition of the 
stress-energy tensor T(z) can be extended into the lower half-plane as T(z*) = T(z). Such an 
extension is compatible with the boundary conditions, because 

T\ b = T\b (4-4.1) 

at the boundary, which in cartesian coordinates means T xy \s = 0, i.e. there is no energy or 
momentum flux across the surface. In this framework, it is possible to show that a n-point 
correlation function { 0i(^i) ••• (j> n ( z n) )hp on the half-plane can be identified with the 2n- 
point function on the entire plane ( 4>i(zi) 4>i(z\) ... 4> n (z n ) ^) n {z^) ) p , in which a mirror image 
is associated to each field. The role of the boundary is therefore simulated by the interaction 
between mirror images of the same holomorphic field. 

The most intriguing feature of boundary CFT is the existence of a relation between boundary 
conditions and the bulk operator content of the theory. In fact, condition (|4.4.1j) is expressed in 
terms of the Virasoro generators acting on a boundary state |a) as 

(L n - L_„) | a ) = , 

and the solutions are the so-called Ishibashi states 

\j)) = ^2\j; N)®UWW), 

N 

where | j; N) and \j; N) are the holomorphic and antiholomorphic states belonging to the 
conformal family of the bulk primary operator labelled by j, and U is an antiunitary operator 
introduced for technical reasons. The boundary states | a ) are therefore linear combinations of 
Ishibashi states associated with different primary operators 

I") = ^2 C aj \j)) , 

i 

and the coefficients C a j can be determined by exploiting the modular invariance of the the- 
ory. This procedure consists in first mapping the half plane on a strip through the conformal 
transformation 

z _> w {z) = —\az, (4.4.2) 

7T 

and then in interpreting the strip geometry in two different physical ways, similarly to what we 
have discussed in Sect. l4.1l for the cylinder (see Fig. l4.9j> . 

In the first physical picture the boundaries are placed at x = 0, R, and the partition function 
is obtained by evolving the Hamiltonian H a p along the time direction, while in the second 



1/3) 

A * 

^ H 
I a ) 



the half plane to the strip. 

scheme space is unbounded, and the infinite-volume Hamiltonian H evolves from a state | a ) 
to another state | j3 ) . The equality of the result in the two different pictures permits to express 
the coefficients C a j in terms of the so-called modular matrix of the system, establishing a 
1 : 1 correspondence between boundary states and primary operators of the bulk theory. For 
instance, in the Ising model (introduced in Sect. II. l|) there are three possible boundary states, 
symbolically indicated as \j) and given by 

T^ l)) + ^ £)) + T2 W)) 

72 |I)) + 72 |e)) "?2 |cj)) 




i) = 

e) = 
~a) = 



Finally, it is useful for the future discussions to illustrate the features of the Gaussian CFT, 
described in Sect. ll.ll defined here on a strip with Neumann or Dirichlet boundary conditions. 
The transformation (|4.4,2[) leads to the following expression for the Hamiltonian on the strip 

H = t; f L o - ■ (4.4.3) 



R V 24 
In the case of Neumann boundary conditions 

d x <p(o,t) = d x <p(R,t) = o , 

there is no winding, and the sectors of the theory are only labelled by the momentum index s. 
The lowest eigenvalue of Lq in each sector is given by 

1 s 2 

h s = ^- . (4.4.4) 

2irg TZ 2 v ' 

On the contrary, Dirichlet boundary conditions 

0(0, t) = <p , </>(R,t) = <j) R 

constrain the momentum index to zero, but leave the winding free. In this case, the lowest 
conformal dimension in each sector is expressed as 



h n = 7T [(<f>R ~ 0o) + 2-ktiR} 2 



(4.4.5) 



4.4.2 Sine— Gordon model on the strip 

The Sine-Gordon model (j2.l-t.ljl can be defined on a strip x G [0, R]. Among the several possible 
boundary conditions which preserve the integrability of the model, a particularly interesting 
example is given by the Dirichlet boundary conditions (D.b.c.) 

2ir 2ir 
0(0, t) = 4>o + yn , (f>(R,t) = <p R + —n R , Vt (4.4.6) 

with < 4>o,R < and tiq^r € Z. The topological charge of this model is conserved also in the 
presence of boundaries and it can be conveniently defined as 

Q = I J dx4> dx - (cf) R - (j) ) | = n R - riQ . 

Hence the space of states is split in topological sectors with Q = 0, ±1, ±2..., and within a given 
Q-sector the states are characterized by their energies only. 

It is worth mentioning that, in recent years, this problem (and variations thereof) has at- 
tracted the attention of several groups: the case of half-plane geometry, for instance, has been 
discussed by bootstrap methods in [7^1 E7J and by semiclassical ones in (751 ESI MH whereas 
the thermodynamics of different cases in a strip geometry has been studied in a series of publi- 
cations (see [8^ E21 EH1 EU E3 ESI EZj ) • In completing the work [3], a paper on a (semi) classical 
analysis of Sine-Gordon model on a strip [HE] also appeared, which partially overlaps with it. 

The semiclassical quantization presented in [I] adds new pieces of information on this subject 
and it may be seen as complementary to the aforementioned studies. For the Sine-Gordon 
model with periodic boundary conditions, alias in a cylindrical geometry, this program has been 
completed in [2| and presented in Sect. l4.3l Given the similarity of the outcoming formulas with 
the ones appearing in in the sequel we will often refer to Sect. l4.3l for the main mathematical 
definitions as well as for the discussion of some technical details. There is though a conceptual 
difference between the periodic example and the one studied here: in the periodic case, in fact, 
the vacuum sector is trivial at the semiclassical level (it simply corresponds to the constant 
classical solution) and therefore the semiclassical quantization provides non-perturbative results 
just starting from the one-kink sector. Contrarily, on the strip with Dirichlet b.c, the vacuum 
sector itself is represented by a non-trivial classical solution and its quantization is even slightly 
more elaborated than the one of the kink sectors. 

In order to describe the classical solutions of this problem, it is worth to preliminarily recall 
some results already obtained. As we have seen in the previous sections, the SG equation of 
motions admit three kinds of static solution, depending on the sign of the constant A in (j4.2.6j) . 
The simplest corresponds to A = and it describes the standard kink in infinite volume (j2.3.3j) . 
In the following, we will be concerned with the solutions relative to the case A ^ 0, which can 
be expressed in terms of Jacobi elliptic functions. In particular, for A > we have 

, 2 ( m(x - x ) r \ , 2 
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.-how ii in Fig. l4,lUl For —2 < A < 0, the solution is given instead by 



which has the monotonic and unbounded behaviour in terms of the real variable u -1 " — r 
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and it oscillates in the real variable u = m(x — xq) between the ^-dependent values (j) and 




u 



Figure 4.10: Solutions of eq. (|4.2.6jl . A > (left hand side), -2 < A < (right hand side). 

The SG model with the Dirichlet b.c. Q4.4.6|) can be classically described by using the two 
building functions <f)^(x) and (j)~,(x), thanks to their free parameters xq and k, which can be 
fixed in terms of <J)q, <f)R and R. However, in order to simplify the notation, in writing down our 
solutions we will rather use R and xq, both considered as functions of (fro, and k (as a matter 
of fact, k can be recovered by inverting the elliptic integrals which enter the corresponding 
expression of R). 

As shown below, both types of solutions <p^{x) and (^{x) are needed, in general, to define 
the classical background in the vacuum sector whereas only one of them, </>^(x), is employed for 
implementing the Dirichlet b.c. in the kink sector. 



The vacuum sector: Q = 



To discuss the vacuum sector, it is sufficient to restrict the attention to the case 9 no = ur = 0, 



(fro < <f>R and 



cos ^4>o 
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COS ^4>R 



It is also convenient to introduce the compact notation 



cq,r 



cos 



In order to write down explicitly the classical background corresponding to the vacuum state 
with Dirichlet b.c, it is necessary to introduce preliminarily two particular values R\ and R2 of 
the width R of the strip, which mark a change in the nature of the solution. They are given by 



mRi = arctanh (co) — arctanh (cr) , 
mR 2 = K(k)-F (arcsin £f , k 
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col 



With these definitions, the classical vacuum solution, as a function of a; £ [0, R], has the following 
behaviour in the three regimes of R: 
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(4.4.7) 



9 A11 other cases can be described in a similar way, denning properly xo and R, and by using antikinks when 
necessary. 
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with 



mx = -2K(fc) + F (arcsin ^ , fc) 

mi? = F (arcsin ^ , kj — F (arcsin ^fi- , k) 

k < k < 1 
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mxo = — F (arcsin , A;) 
4 ?ni? = 2K(fc) - F (arcsin ^ , k) — F (arcsin , jfe) 
< < 1 



In addition to the quantities already defined, we have introduced here the incomplete elliptic 
integrals 

F((p,k) = 

o ' " " - " o 



da 

yl — k 2 sin 2 a 



E(ip, k) = J da^/ 1 - k 2 sin 2 a 



which reduce to the complete ones at <p = tt/2. It is easy to check that at the particular values R± 
and R2, the different definitions of the background nicely coincide. Fig. 14. Ill shows the classical 
solution at some values of R, one for each of the three regimes 10 . 




Figure 4.11: Classical solution ()4.4.7|) at some value of R, in the case ficfio = 1 and f3(j)R = 2. 



The classical energy of the background Q4.4.7JI is expressed as 
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We have chosen for the plot the specific values fl(f>o — 1 and j3<j)R = 2, for which mRi = 0.76 and mR2 = 1.49. 
The same values will be considered in all other pictures since their qualitative features do not sensibly depend on 
these parameters, except for few particular values of <j>o,R discussed separately 
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and it is plotted in Fig. l4.12l As expected, the quantity Q4.4.8JI has a smooth behaviour at R\ 
and i?2i which correspond to the minimum and the point of zero curvature of this function, 
respectively. The non monotonic behaviour of the classical energy gives an intuitive motivation 
for the classical background being differently defined in the three regimes of R. 
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Figure 4.12: Classical energy ()4.4.8|) for 04>o = 1 and 04>r = 2. 

Furthermore, the classical energy can be easily expanded in the ultraviolet (UV) or infrared 
(IR) limit, i.e. for small or large values of mR, which correspond to k — > in the regime 
< R < R\ or to k — ► 1 in the regime R2 < R < 00, respectively. 

In fact, expanding the elliptic integrals in (|4.4.8|) (see [013] and Appendh d4lljl for the relative 
formulas), and comparing the result order by order with the small- A: expansion of mR defined 
in the first regime of (|4.4.7|) 
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one obtains the small-mi? behaviour 
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Later we will comment on the meaning of this result in the UV analysis of the ground state 
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energy. On the other hand, comparing the expansion for k — > 1 of £\ (R) in the third regime 
with 
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(4.4.11) 



one obtains the large-mi? behaviour 

4 3) (^) = |? (2-cos^«-cos^o)-^Ptan^otan^i J e- miJ + -- - . (4.4.12) 

The first term of this expression is the classical limit of the boundary energy of the vacuum 
sector jS^, since it is the term that needs to be subtracted by choosing to normalise the energy 
to zero at R — > oo. 

The classical description of the vacuum sector can be completed by mentioning the existence 
of two particular cases in which the three different regimes of R are not needed. The first is 
given by (fro = cj)R, for which the whole range of R is described by (fr cl (x) in (|4.4.7|) . since 
mi?2 = in this situation. The second case, defined by (fro arbitrary and (fru = 0, can be instead 
described by the antikink (fr^\x) = (fr^\—x) alone, since mRi = oo for these values of the 
boundary parameters (note that xo and R have to be defined as opposite to the ones in ()4.4.7|0 . 
As a consequence, these two cases display a monotonic behaviour of the classical energy, whose 
UV and IR asymptotics, respectively, require a separate derivation, which can be performed by 
simply adapting the above procedure. 

Finally, it is also worth discussing an interesting feature which emerges in the IR limit of the 
classical solution (|4, 4.7)1 , As it can be seen from Fig. 14. Ill by increasing R the static background 
is more and more localised closely to the constant value (fr{x) = and this guarantees the 
finiteness of the classical energy in the R — > oo limit, given by the first term in (|4,4.12l) n . 
However, if the IR limit is performed directly on the classical solution, we obtain one of the 
static backgrounds 12 studied in jSUj 

cfr Mx) — > — arccos [tanh m(x — Xq* 3 )] , with = — arctanh(co) . 

R^oo p 

The last expression tends to zero as x — > oo and consequently has classical energy E c i = 
^jpr ^ 1 — cos j ■ This phenomenon can be easily understood by noting that the minimum 

of (frfi(x) (which goes to zero in the IR limit), is placed at mx = tjixq + K(fe) (see Fig. l4.10l) 
and this point tends itself to infinity as k — > 1. Hence, the information about the specific value 
of <j)R is lost when R — > oo, i.e. only the states with cfru = survive in the IR limit. 

We will now perform the semiclassical quantization in the vacuum sector, around the back- 
ground (|4.4.7jl . Depending on the value of mR, the stability equation (|2.1.9|) takes the form 
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when R\ < R < R2 and R2 < R < 00. 

Equations ()4.4.1«-{|) and ()4.4.14)) can be cast in the Lame form with N = 1, which is described 
in Appendix I4ICJI and has been studied in detail for the periodic case in Sect. 14 .31 The only 

11 When I Co | < |cr|, the same qualitative phenomenon occurs, but the constant value is <f){x) = ^ in this case. 
12 Obviously, the same function is obtained as lim <jy\ (x), in the case 4>R = mentioned above. 
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differences with the periodic case are the presence of a non-trivial center of mass xq and the 
larger number of parameters entering the expression of the size R of the system: these make more 
complicated the so-called "quantization condition" that determines the discrete eigenvalues, 
although they do not alter the general procedure to derive it. 

The boundary conditions (|4.4.fi|) . which translate in the requirement 

ifc>(0) = VaiR) = , 
select in this case the following eigenvalues, all with multiplicity one, 
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and the y n 's are defined through the "quantization condition" 
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This equation comes from the consistency condition associated to the boundary values 

(D + Va (0) + D_ »7_ o (0) =0 , 
{D +Va (R) + D- V - a {R) = , 

where rj± a are the two linearly independent solutions of the Lame equation which are used to 
build the general solution rj{x) = D + r] a (x)+D^ r/_ a (x) (see Appendix l4ICI for details). Along the 
same lines discussed for (|4.3.2Uj) in the periodic case, the requirement (|4.4,16|) can be physically 
interpreted as a quantization condition for the momentum of a state containing a neutral particle 
between the two Dirichlet boundary states. The agreement between the large-i? limit of this 
expression and the Bethe ansatz equations involving the appropriate reflection matrices has been 
verified in [EE]- 

As it can be seen directly from 1)2.1.10(1 . the frequencies 1)4.4.15)) are nothing else but the 
energies of the excited states with respect to the ground state Eq &c (R). They can be easily 
determined from the above equations and their behaviour, as functions of R, is shown in Fig. 14. L-fl 

As in the periodic case, a more explicit expression for the energy levels 1)4.4.15)1 can be 
obtained by expanding them for small or large values of mR. The UV expansion, for instance, 
can be performed extracting from 1)4.4.16)1 a small-/c expansion for y n , inserting it in 1)4.4.15)1 . 
and finally comparing the result order by order with l)4.4.9)) . Exploiting the several properties 
of Weierstrass functions which follow from their relation with (9-functions (see for instance [HJ ) , 
one gets 
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Figure 4.13: The first few energy levels 1)4.4.15)1 for (3<fio = 1 and 04>r = 2. 
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where we have introduced the compact notation / = (3(4>r — cfto), s = (sin/3(/>^ — sin/?^o)- This 
leads to the UV expansion 
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(4.4.17) 



In order to complete the above analysis and obtain the reference value of the energy levels, 
i.e. the ground state energy Eq 3C (R) of the vacuum sector, we need the classical energy (|4.4.8|) 
and the sum on the stability frequencies given in (|4.4.15|) . i.e. 
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(4.4.18) 
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The above series is divergent and its regularization has to be performed by subtracting to it a 
mass counterterm and the divergent term coming from the infinite volume limit - a procedure 
that is conceptually analogous to the one discussed in Sect. l4.3l for the periodic case and therefore 
it is not repeated here. Furthermore, as already mentioned, equation 1)4.4. 18|) can be made 
more explicit by expanding it for small or large values of mR. Here, for simplicity, we limit 
ourselves to the discussion of the leading 1/R term in the UV expansion since it does not receive 
contributions from the counterterm and therefore it can be simply regularised by using the 
Riemann ("-function prescription (see Sect. l4.3T1 for a detailed discussion). The higher terms, 
instead, require a technically more complicated regularization, although equivalent to the one 
presented in [3]. 

The UV behaviour of the ground state energy is dominated by 
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where the coefficient —1/24 comes from the regularization of the leading term in the series of 
frequencies ()4.4.17|) , while the first term simply comes from the expansion of the classical energy 



(|4.4.1U|) . It is easy to see that the above expression correctly reproduces the expected ground 
state energy 1)4,4.3114.4.5)) for the gaussian Conformal Field Theory (CFT) on a strip of width 
R with Dirichlet boundary conditions. 

Finally, it is simple to check that also the excited energy levels display the correct UV 
behaviour, being expressed as 
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The kink sector: Q = 1 

In discussing the kink sector we can restrict to no = , tir = 1, since all other cases, as well 
as the antikink sector with Q = — 1, are described by straightforward generalizations of the 
following formulas. 

The classical solution can be now expressed only in terms of the function cft^ (x) as 
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(4.4.21) 

since in this case the whole range < mR < oo is spanned by varying k in [0, 1]. This can be 
intuitively understood by looking at the behaviour of 1)4.4.21)1 in Fig. l4.14l 
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Figure 4.14: Classical solution 1)4.4.21)1 at some values of R, in the case 



1 and (3(f>R = 2. 



As a consequence, the classical energy and the stability frequencies of this sector can be 
obtained from sjp and lo^ of the vacuum (given respectively in eq. 1)4.4.8)1 and 1)4.4.15)1 ). by 
simply replacing — ► + =jj-. The leading UV behaviour of the energy levels in this sector, 
given by 
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with Q = 1, correctly matches the CFT prediction. 

The only result which cannot be directly extracted from the vacuum sector analysis is the 
IR asymptotic behaviour of the classical energy, since now the k — > 1 limit has to be performed 



on sjp . We have in this case 
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which leads to 
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Analogously to the vacuum sector, the first term of this expression is related to the classical 
limit of the boundary energy in the one-kink sector. Notice that, differently from the vacuum 
case, where the asymptotic IR value of the classical energy was approached from below (see 
(|4.4.12|) ). the coefficient of the exponential correction has now positive sign, in agreement with 
the monotonic behaviour of the classical energy shown in Fig. l4.15l 
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Figure 4.15: Classical energy in the Q = 1 kink sector for (3<j)Q = 1 and (3(pR = 2. 

When R — > oo, a mechanism analogous to the one discussed for the vacuum also takes place 
here: the classical energy is finite for any value of (j)R, but since ^ mk (x) assumes the value ^ 
at mx = mxQ + kK(k) (see Fig. l4.lUll . a point which tends to infinity as k — ► 1, only the states 
with (ftn = survive in this limit. 

It is worth noticing that <j)^~' (x) can be also used to satisfy, at finite values of R, Dirichlet 
b.c. in sectors with arbitrary topological charge (see Fig. 14.16]) . giving rise to the correct UV 
behaviour (|4.4.22jl with Q = ur—uq. However, since <p£~' \x) always assumes the value ^(no + l) 
at mx = mxQ + kK(k), which is once again the point going to infinity when k — ► 1, in the IR 
limit it can only correspond to Q = 1. This result seems natural though, since in infinite 
volume, static classical solutions can only describe those sectors of the theory with Q = 0, ±1, 
while time-dependent ones are needed for higher values of Q. Hence, in the topological sectors 
with \Q\ > 1 the space of states will contain, at classical level, the time-dependent backgrounds, 
defined for any value of R (which are not discussed here), plus the static ones of the form (^'(x), 
which however disappear from the spectrum as R — > oo. 
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Figure 4.16: Classical solution in the Q = 3 sector (no = 0, ur = 3) at some values of R, in the 
case P4>q = 1 and fi<\>n = 2. 



4.5 Summary 

In this Chapter we have shown how the semiclassical methods can provide an analytic description 
of finite size effects in two-dimensional quantum field theories displaying degenerate vacua. In 
particular, we have applied the form factors technique to study the correlation functions in SG 
and broken (j) 4 theories on a cylindrical geometry. Furthermore, implementing on the cylinder 
the DHN quantization method, we have derived the scaling functions of the ground (and excited) 
states in the one-kink (with Q = 1) sector of the SG model. The semiclassical approach provides 
analytic and non-perturbative expressions for the energy levels, valid for arbitrary values of the 
size R of the system, which permit to link the IR data of the massive theory with the UV 
conformal data of CFT. It is particularly interesting the application to non-integrable models, 
where the large R behaviour of the obtained scaling functions can be compared with Liischer's 
theory to extract information about the unknown scattering properties of the model, as we have 
shown in the broken <^> 4 case at the level of classical energy. 

For the integrable case of sine-Gordon model, the next step in the semiclassical program is 
the extension of the DHN method to describe the multi-kink states ( Q = ±2, ±3, ...) as well as 
the non-vacua ("breather" -like) part of the Q = sector. These states are related to certain 
time-dependent solutions on the cylinder, i.e. to the finite volume analog of soliton-soliton, 
soliton-antisoliton and breather solutions. Although more complicated from the technical point 
of view, the determination of these classical solutions and the study of their scaling functions 
and form factors is a well stated open problem in the semiclassical framework, which deserves 
further attention. Another interesting problem to be studied is related to the the exactness 
(or very high accuracy) of the semiclassical results, which is a peculiar feature of SG model in 
infinite volume, as we have observed in Sect. 12.31 It would be very useful to understand whether 
similar phenomena take place for the semiclassical scaling functions and form factors in finite 
volume as well. An indication on this issue could be found by extending to finite volume the 
analysis of higher loop quantum corrections in the semiclassical expansion. 

The semiclassical method is suited also for the description of finite geometries with bound- 
aries, and the example of the SG model on a strip with Dirichlet b.c. has been considered. The 
resulting analytic expressions for the energy levels permit to link the IR data on the half-line 
with the UV conformal data of boundary CFT at c = 1. In comparison with a cylindrical 
geometry, an interesting new feature of the quantum field theory defined on a strip consists in 



a non-trivial (and non-perturbative) semiclassical description of its vacuum sector. Therefore, 
we have discussed in detail the classical solutions and energy levels in the Q = case, together 
with the Q = 1 that can also be described by static backgrounds. As we have already men- 
tioned, however, the semiclassical methods are not restricted to static backgrounds only, and 
a complete description of the theory in all sectors requires also the study of time-dependent 
solutions. Finally, it is worth noticing that the analysis performed here has natural and direct 
extension to other quantum field theories with various kinds of boundary conditions. 

One of the advantages of the semiclassical method is that it works equally well for both 
integrable and non-integrable models, if they admit kink-type solutions. In fact, we have 
chosen to test the efficiency of the semiclassical quantization on the example of SG model, 
mainly because it leads to the simplest N = 1 Lame equation. Static elliptic solutions for 
other models can be easily obtained by integrating equation (|2.1.4|) with ^4^0 and appropriate 
boundary conditions. This was done, for instance in pQ, where we have derived the form factors 
between kink states in the broken </> 4 model on the cylinder with twisted boundary conditions. 
In this case, the quantization of the finite volume kink involves a Lame equation with N = 2 and 
will be presented in Lame equations with N > 2 are also expected to enter the quantization 
of other theories. 



Appendix 4. A. Free theory quantization on a finite geometry 



Let us consider a free bosonic field (j)(x, t) of mass m defined on a cylinder of circumference 
R, i.e. satisfying the periodic boundary condition 

4>(x + R,t) = cf)(x,t) . (4.A.1) 

Imposing the equation of motion and the commutation relation 

[4>(x,t),U(y,t)} = i5 P {x-y) , 

where H(x,t) = ^(x,t) is the conjugate momentum of the field whereas 

1 oo 

S P (x) = - ^ X . S P (x + R) = 5 P (x) 

n=—oo 

is the periodic version of the Dirac delta function, we obtain the mode expansion of the field 
4>(x,t). This is given by 
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(4.A.3) 



Using the above expansion together with the commutation relation of A and A\ it is easy to 
compute the propagator of the field, given by 



oo ^ 

A F (x-x',t-t') = (0(x,t)<p{x',t')) = y- 



-i[u) n (t-t')-p n (x-x')} 



2u> n R 

The vacuum expectation value of the operator (p 2 (x,t) is then formally given by 

{<P 2 (x,t)) = A F (0) 



(4.A.4) 



(4.A.5) 



and, by translation invariance, is independent from x and t. However this expression is divergent 
and needs therefore to be regularized. Analogously to what has been done in the text for the 
ground state energy <?Q ac (it!), we need to subtract the corresponding expression in the infinite 
volume, so that the finite quantity, simply denoted by (Pq(R), satisfies the usual normalization 
condition 



lim <g(R) = . 

H— >oo 



Hence we define 



Isolating its zero mode, the series needs just one subtraction, i.e. 
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(r = mR). In the above expression, the first series is now convergent whereas the second series, 
which is divergent, has to be combined with a divergence coming from the integral. Indeed we 
have 

CO 

T(r) = / dn — — = lim jln2A — In — }— lim lnA+ lim In A, 

J 2 _i_ ( r \ 2 A^co I 2lT ) A—>oo A-^co 







and the last term can be used to compose (|4.3.7|) . Collecting the above expressions, it is now easy 
to see that 4>q(R) coincides with the one obtained doing the calculation in the other quantization 
scheme, i.e. at a finite temperature. In fact, using the results of Ref. |65| . this quantity can be 
expressed as 

/oo m i 
„5F7=*rr ' (4A ' 7) 

whose expansion in r is given by 

llr °° ( 1 1 \ 

S 2 JR) = — + — ( log — + j E - log 9) + . (4.A.8) 
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Also this result could have been directly obtained computing only the finite part of the 
integral and using the prescription (|4.3.1U|) . 



Appendix 4.B. Elliptic integrals and Jacobi's elliptic functions 



In this appendix we collect the definitions and basic properties of the elliptic integrals and 
functions used in the text. Exhaustive details can be found in |90| . 

The complete elliptic integrals of the first and second kind, respectively, are defined as 



tt/2 



K(k 2 ) 



da 



\/l — k 2 sin 2 a 



E( 



tt/2 

l ) = j (Wl - k 2 sin 2 a . (4.B.1) 



The parameter k, called elliptic modulus, has to be bounded by k 2 < 1. It turns out that the 
elliptic integrals are nothing but specific hyper geometric functions, which can be easily expanded 
for small k: 
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Furthermore, for k 2 — > 1, they admit the following expansion in the so-called complementary 
modulus k' = y/l — k 2 : 
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Note that the complementary elliptic integral of the first kind is defined as 

K'{k 2 ) = K(k' 2 ) . 

The function am(it, k 2 ), depending on the parameter k, and called Jacobi's elliptic amplitude, 
is defined through the first order differential equation 



dam(ii) 
du 



1 — k 2 sin 2 [am(u)] , 



(4.B.2) 



and it is doubly quasi-periodic in the variable u: 

am (u + 2nK + 2imK') = mr + am(ii) . 
The Jacobi's elliptic function sn(u, k 2 ), defined through the equation 



(^) 2 = (l-sn 2 .)(l-^ 2 sn 2 .) , 



(4.B.3) 



is related to the amplitude by snu = sin (amn), and it is doubly periodic: 

sn (u + 4nK + 2imK') = sn(i 



Appendix 4.C. Lame equation 

The second order differential equation 

S^-E-N(N + l)V(u)\f(u) = 0, (4.C.1) 

where E is a real quantity, N is a positive integer and V(u) denotes the Weierstrass function, is 
known under the name of iV-th Lame equation. The function V(u) is a doubly periodic solution 
of the first order equation (see |90j ) 

dV\ 2 

= 4(7>- ei ) (V-e 2 ) (V-e 3 ) , (4.G.2) 

whose characteristic roots ei, e2, e3 uniquely determine the half-periods cj and a;', defined by 

7> (it + 2nw + 2raw') = V(u) . 

The stability equation (|4.3,16|) can be identified with eq. (|4.C,1|) for N = 1, u = x-\-i~K! and 
E = 2 =f-k 2 uj 2 in virtue of the relation between V(u) and the Jacobi elliptic function sn(u, k) 
(see formulas 8.151 and 8.169 of 

k 2 + 1 

k 2 sn 2 {x, k) = V(x + iK!) + — — . (4.C.3) 

Relation ()4.(".3|) holds if the characteristic roots of V{u) are expressed in terms of k 2 as 

2 - k 2 2k 2 -1 1 + k 2 

ei = , e9 = , e3 = , 

3 3 3 

and, as a consequence, the real and imaginary half periods of V(v) are given by the elliptic 
integrals of the first kind 

to = K(fc) , J = iK'(k) . 

All the properties of Weierstrass functions that we will use in the following are specified to the 
case when this identification holds. 

In the case N = 1 the two linearly independent solutions of (|4.C.1|) are given by (see |9T| ) 

f±a{u) = , (4.C.4) 

where a is an auxiliary parameter defined through V(a) = E, and cr(u) and C( u ) are other kinds 
of Weierstrass functions: 

^ --*«>■ ^-««). (^) 

with the properties 

C(n + 2K) = C(u) + 2C(K) , 

a(u + 2K) = - e 2 («+ K )C( K )(7(it) . (4.C.6) 
As a consequence of eq. ([4.C.6)) one obtains the Floquet exponent of f± a (u), defined as 

f(u + 2K) = f(u)e iF ^ , (4.C.7) 



in the form 

F(±a) = ±2*[KC(a)-aC(K)] . (4.C.8) 

The spectrum in the variable E of eq. ()4.C.1|) with N = 1 is divided in allowed/forbidden 
bands depending on whether F(a) is real or complex for the corresponding values of a. We 
have that E < and e2 < E < e\ correspond to allowed bands, while < E < e<i and 
E > ei are forbidden bands. Note that if we exploit the periodicity of V{a) and redefine 
a — > a' = a + 2nco + 2mui' , this only shifts F to F' = F + 2rmr. 

The function ((u) admits a series representation [52] that will be very useful for our purposes 
in Sect. EH 
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where /i = e _7rK //K . The small- A; expansion of this expression gives 
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(note that h ~ (i) 2 + 0(A; 4 )). A similar expression takes place for V(u), by noting that 

Ttf \ dC(u) 
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Finally, in the complementary limit k' = y/l — k 2 — > 0, the Weierstrass functions can be 
expanded as 
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Conclusions and Outlook 



In this thesis, we have described some fruitful applications of the semiclassical methods in the 
study of non-integrable QFT and finite-size effects in two dimensions. 

The mass spectrum of a generic QFT can be explored non-perturbatively by exploiting 
the analytic properties of semiclassical form factors of the local fields between kink states, 
which are available as the Fourier transforms of the classical kinks. With this technique, we 
have specifically studied two non-integrable theories, given by the </> 4 interaction in the Z2 
broken symmetry phase and by the double sine-Gordon model, which describe relevant physical 
phenomena and find many applications in statistical mechanics and condensed matter physics. 
Along similar lines, it would be worthwhile to extend the semiclassical analysis to other non- 
integrable theories which describe the scaling region of interesting statistical systems, like for 
instance the </> 6 potential with three degenerate vacua, associated to the tricritical Ising model 
with both leading and subleading thermal perturbations. Moreover, further attention should be 
dedicated to the false vacuum decay in presence of unstable classical backgrounds, that we have 
studied in the specific case of the double sine-Gordon model. In fact, the semiclassical method 
is capable of providing interesting results about the dynamics of decay processes. 

An opportune generalization of the semiclassical quantization technique has led us to the ana- 
lytical study of finite-size effects in QFT. In fact, once the proper classical solutions are identified 
on the geometry of interest, one can implement on them the relativistically refined Goldstone 
and Jackiw's result and the DHN quantization method in order to estimate, respectively, the 
correlators and the spectrum. In particular, we have studied form factors and correlation func- 
tions for the SG and broken 4 theories on a twisted cylinder, and we have derived the scaling 
functions of the ground (and excited) states in few sectors of the SG model defined both on a 
cylinder and on a strip with Dirichlet boundary conditions. An interesting open problem is the 
systematic investigation of the finite-size spectrum, which consists in finding time-dependent 
classical solutions in the remaining topological sectors of the theory, and in applying to them the 
semiclassical quantization procedure. Furthermore, it would be useful to compare numerically 
the semiclassical results with the ones available from different techniques for the SG model in 
finite volume. This, in fact, would permit to test quantitatively the efficiency of the semiclas- 
sical approximation, as we have successfully done in the infinite volume case. A closely related 
subject is the study of two-dimensional QFT at finite temperature, which are known to repre- 
sent the physical interpretation of the cylindrical geometry in which the compactified variable 
is Euclidean time, instead of space. A more profound understanding of the relation between 
finite-volume and finite-temperature pictures will give the opportunity to gain further insight in 
these both quantization schemes. 

As a final remark, it is worth to emphasize again that the semiclassical method is equally 
suited to describe integrable and non-integrable models, if they admit kink-type solutions. 
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This feature, that we have fully exploited in the analysis of infinite-volume spectra, makes the 
analytical study of scaling functions in non-integrable theories a well stated problem, more 
complicated only from the technical point of view, but conceptually analogous to the procedure 
illustrated for the sine-Gordon model. 
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